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Automatic equipment providing water for intravenous and other parenteral solutions in Jefferson Hospital, Philadelphia. One of 
the two stills is for standby service. Micromax recorder and cell, which record purity of distillate from still in use, are 
shown close-up at left. 


Hospital Has Continuous Proof of Water Purity 


Among many interesting features of the Jefferson Hos- 
pital centralized sterilization room is an automatic device 
for checking and recording the purity of distilled water. 


Water is provided by a ten-gallons-per-hour Barnstead 
automatic still. The effluent flows through a conductivity 
cell which is connected to a Micromax strip-chart recorder. 
This combination measures and records the electrical re- 
sistance of the distillate. A resistance of $2,000 ohms per 
cc. has been adopted by this hospital as satisfactory for 
parenteral use, and since the conductivity cell used with 
the equipment has a constant of 0.1, the resistance shown 
by the recorder becomes 8200 ohms. When the still is 
started and until the resistance is above this figure, the 
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red light is on. If the effluent should drop again to that 
figure, the recorder would flash a red light to summon 
one of the station nurses. 


Principal advantages of this method of analysis are 
that it largely eliminates the human equation, and that it 
provides a graphic record, always available for consulta- 
tion. It checks every drop of water instead of occasional 
grab samples. It wastes no water. Since the purity is 
recorded, the station supervisor can watch not only for 
sudden or temporary changes but for gradual ones due 
to aging of equipment, seasonal differences in raw water, 
etc. 


Catalog EN-95, sent on request, describes the equipment. 
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“IS THE 
STANDARD CELL 
EPPLEY CELL?” 


The far-seeing purchasing agent 
will find it to his advantage to ask 
this question before making a deci- 
_ sion on the purchase of automatic 

temperature control apparatus be- 
cause— 


For nearly a quarter of a century the 
_ leading American manufacturers of 
such apparatus have used Eppley 
Standard Cells in their recorders, py- 
rometers and other instruments utiliz- 
ing the potentiometric 
measurement. 


system of 


Today, more than ever before, these 
_ manufacturers realize that the standard 
cell—the heart of the potentiometric 
system—must be the finest and most 
reliable cell obtainable. 


The technical men responsible for the 
performance of the instruments know 
that the Eppley Standard Cell meas- 
ures up fully to their standards of 
excellence and reliability. 


THE EPPLEY LABORATORY, INC. 


SCIENTIFIC INSTRUMENTS 


| NEWPORT, R.1. 
| U.S.A, 


INCREASE 
EFFICIENCY 


UNCOATED COATED 


SPECTROSCOPES 
MICROSCOPES 
OPTICAL INSTRUMENTS 


By high vacuum coating with suitable fluoride film* 
to minimize surface reflection loss shown above in 
half coated aircraft instrument. 


NATIONAL RESEARCH CORPORATION 
100 BROOKLINE AVE. 
BOSTON, MASS. 


*U. S. Pat. 2,207,656 
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| EXCEPTIONALLY WIDE RANGE—21 to 1140 db 
2 SIMPLIFIED CONTROLS—complete db range cov- 


ered by a single control knob 
. 3 IMPROVED MICROPHONE—special type of great 
sensitivity, highly rugged, unaffected by ordi- 
nary changes in temperature and humidity 
4 PLUG-IN MICROPHONE—normally used directly 
. on folding socket on instrument but can be de- 
tached and used with cable and tripod if desired 


5 SLOW-FAST METER—two speeds—control switch 
TYPE 759-B SOUND-LEVEL METER 


@ Write for Bulletin 
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Simplified Sound Measurements 


With this Accurate, Highly Stable Meter 


ITH THE IMPROVED General Radio Sound-Level Meter the intensity of sound can now be 
measured conveniently, quickly and accurately by anyone. 
throughout industry for all types of noise measurements and studies. 
tained and portable and will measure the intensity of sound from a whisper to a steam whistle. 


FEATURES: | 


GENERAL RADIO COMPANY. Cambridge. Massachusetts 


Please mention JOURNAL OF APPLIED PHYSICS when writing to advertisers 


This instrument is widely used 
It is completely self-con- 


selects either; “fast” position, meter movement 
corresponds to A. S. A. standards; “slow” posi- 
tion, meter heavily damped to measure average 
level fluctuations 
6 INTERNAL CALIBRATION SYSTEM—calibration can 
be checked quickly from built-in calibrator 
7 THREE WEIGHTING NETWORKS—low level, high 
level and substantially flat over-all response 
8 NO COILS—instrument contains no coils or in- 
ductances and accordingly reads accurately in 
the presence of any ordinary magnetic field 


$195.00 


718 for Information 


ea 
ot ® 2 
SS 
| 


WIRE WOUND 


July, 1941 


Never has there been such a con- 
certed development of electrical equip- 
ment requiring resistors for bridge and 
similar circuits where the call is for 
greater accuracy and lower temperature 
co-efficient than heretofore available in 
production quantities. 


IRC Precision Wire Wound Resistors 
were designed for such purposes. They 
are made to a standard tolerance of 
+1% or to as low as 1/10 of 1% on special 
order. They have an extremely low 
temperature co-efficient for stable opera- 
tion over a temperature range from 
-30° C. to 100° C. They are impreg- 
nated against humidity and other atmos- 
pheric conditions. They are available in 
the most complete line of inductive and 
non-inductive types, shapes, sizes and 
terminals on the market today. 


MAXIMUM PROTECTION © CLOSE TOLERANCES * SMALLER SIZES 


This is an actual size view of an IRC Type BW ‘2-watt Insulated Wire Wound 
Resistor —a compact, well-protected unit that can often be used in exacting 
applications where space is at a premium. Other Type BW sizes are 1-watt 
and 2-watts. Tolerances available to as close as 2%. Some sizes and ranges 
can be supplied to even closer tolerances and in matched pairs to 1%. 


The Answer to Many Bridge Circuit 
Resistance Requirements 


Thousands of these units used in a wide 
range of equipment from voltmeter multi- 
pliers for sea coast duty to geophysical 
instruments used in tropical climates 
prove their dependability under the most 
exacting conditions. 


Whatever your resistance require- 
ments, it pays to ask IRC. IRC engineers 
will welcome the opportunity to cooperate. 
They are backed by the world’s largest 
line of fixed and variable resistors. Their 
recommendations are based on many 
years of specialized experience in deal- 
ing with practically every type of resist- 
ance problem. 


Wnite for IRC Resistance Engineering 
Data Bulletin No. IV including both IRC 
Precision and Cement Coated Power 
Wire Wound Resistors. 


INTERNATIONAL RESISTANCE CO. 
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RADIOLOCATORS 


N a radio appeal* on June 18 for volunteers 

for civilian technical workers Lord Beaver- 
brook, British Minister of State, released some 
of the secrets of the “radiolocator,’’ a device 
which, more than anything else, enabled Great 
Britain to ward off the ferocious Nazi air attack 
last fall. According to statements released to the 
press, the radiolocator sends out radio-waves 
which are essentially unaffected by fog, clouds 
or darkness. Any object in the path of the waves, 
such as an enemy plane, reflects back a signal 
which locates the plane and gives instant warn- 
ing of attack. 

A statement by the British Air Ministry in- 
cludes the following: ‘‘The development of the 
radiolocator has created a new military science 
which came into being without the world being 
aware of its existence and which has already 
exerted a profound influence on air, military and 
naval strategy. Today production lines in the 
plants of radio manufacturers throughout Great 
Britain are turning out radiolocators in ever- 
increasing quantities, and their installation 
throughout the country has made it imperative 
to train more men to 
service and maintain 


at least two months of specialized training. 
Thousands of men with a knowledge of radio 
are required to look after the radiolocators now 
being produced. The use of the radiolocator 
in this system has been so successful that inde- 
pendent observers who have seen it in opera- 
tion have expressed the opinion that, but for 
it, the first Battle of Britain might have been 
lost many months ago. But the first battle 
was not lost, and the system of radiolocation 
is now being widely extended. Radio men, ac- 
cepted for this work, receive a course of tech- 
nical instruction lasting approximately eight 
weeks. However much they may already know 
about radio, this course and the experience that 
will follow will extend their knowledge of new 
developments in radio far beyond anything yet 
imagined... .” 

The British themselves need immediately tens 
of thousands of operators and service men for 
radiolocators now coming off the production 
lines. What our own needs will be it is difficult 
to estimate, but we can be sure they will be 
enormous. There seems to be no question but 
that we must train 
many more workers in 


them, although thou- 
sands are already en- 
gaged in this work... . 
They are not simple 
pieces of apparatus. 
Before he is competent 
to maintain one a man 
needs to have a sound 
working knowledge of 
radio to begin with, 
and then he must have 


*For more details of this 
appeal see the New York 
Herald Tribune or The New 
York Times for June 18 and 
19, 1941. 
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electronics and _ radio. 
Funds for the cost of 
instruction will un- 
doubtedly be available 
through legislation 
sponsored by the U.S. 
Office of Education and 
now before Congress. 
Unfortunately, as yet 
no federal provision is 
being made for scholar- 
ships or fellowships to 
attract some of the 
most outstanding grad- 
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and several other short papers and abstracts. uates to advanced 


study in this field. 
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The Physical Aspects of Boundary Lubrication 


By OTTO BEECK 


Shell Development Company, Emeryville, California 


What is Boundary Lubrication? 


F two lubricated surfaces are sliding over 

each other under conditions of high pressure 
and low speed, no truly fluid film of the lubricant 
can exist between the sliding surfaces, and the 
term “‘boundary lubrication’? has been widely 
adopted for this condition. Boundary lubrication 
or non-hydrodynamic lubrication is therefore 
readily distinguished from film’ or 
“hydrodynamic” lubrication, which takes place 
under maintenance of an oil film of sufficient 
thickness so that hydrodynamic properties such 
as viscosity and the relative motion of the two 
surfaces are the sole characterizations necessary. 
Under boundary conditions, phenomena present 
themselves which must be associated with 
physical characteristics of the lubricant other 
than viscosity. They are, in fact, associated with 
chemical as well as physical characteristics of 
both the lubricant and the surface, which makes 
boundary lubrication a very elusive and compli- 
cated scientific problem. Lack of fundamental 
knowledge in this field together with a strong 
realization of its practical importance led some 


time ago to the term “‘oiliness,”’ which simply 
indicates the lubricating value of an oil under 
boundary condition; this term is still widely 
used. Its definition, as adopted later by the 
Society of Automotive Engineers, is as follows: 

“Oiliness is a term signifying differences in 
friction greater than can be accounted for on the 
basis of viscosity when comparing different 
lubricants under identical test conditions.” 

The state of boundary lubrication should, 
however, be defined much more stringently. In 
the author's opinion, the best criterion for this 
state is that the coefficient of friction is independ- 
ent of viscosity and of the sliding velocity. Pure 
boundary lubrication thus defined is extremely 
rare in practice because quasi-hydrodynamic 
effects of the lubricant are difficult to eliminate 
but we shall see that it can be obtained in the 
laboratory. 
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Considerable progress has been made within 
the last 10 years in the understanding of the 
phenomena of boundary lubrication and a good 
deal has been contributed by physicists and 
physical chemists. It is primarily their efforts 
which will be discussed in this summary without 
any claim of completeness. 


The Area of Contact Between Stationary 
and Sliding Surfaces 

It is obviously quite important that before 
speaking of friction between sliding surfaces we 
should have some idea of the actual area of 
contact of such surfaces under various loads, 
since it is this area which determines the actual 
pressure exerted on a protective boundary film 
between the surfaces. Bowden and Tabor! (1939) 
have made an extensive study of this subject. 
They found that the real area of contact of dry 
surfaces is not greatly affected by size, shape 
and even degree of roughness of the surface and 
that it depends mainly on the pressure. This is 
true also for surfaces which are very carefully 
ground and polished. The real area of contact 
was estimated by electrical and by visual means 
and the following results were obtained for steel 
flats of 21 cm’. 


FRACTION OF MACROSCOPIC 


Loap, KG AREA IN CONTACT 
300 1/130 
100 1/700 
20 1/10,000 
3 1/170,000 


This means that even for lightly loaded surfaces 
the local pressure may be sufficiently great to 
cause plastic deformation of the steel. A careful 
investigation of this point through measurements 
of the conductivity of crossed steel or copper 
cylinders leaves no doubt that plastic deforma- 
tion takes place at the local points of contact 
and that the real area of contact between flat 
surfaces is proportional to the normal load, 
whereas on the basis of elastic deformation the 
real area of contact should vary as the two-thirds 


JOURNAL OF APPLIED PHYSICS 


| 
> 


power of the normal load. This, at once, gives 
an explanation of the old law of Amontons 
(1699) which states that the tangential frictional 
force is directly proportional to the normal load, 
the mysterious meaning of which had only 
recently been discussed by Adam? (1938). It 
also explains the well-known other observation 
of Amontons that the friction is independent of 
the macroscopic area of the surface. 

When the real area of contact of two sliding 
surfaces is measured by the electrical conduc- 
tivity methods, rapid fluctuations occur as 
would be expected from an intermittent making 
and breaking of the contacts. However, the 
type and the magnitude of the fluctuations remain 
essentially the same when the surfaces are lubricated 
under boundary conditions with mineral oils or 
other lubricants. 


The Surface Temperature of Sliding Metals 


It is clear from the preceding paragraph that 
when two surfaces slide over each other, the 
actual pressure on the real area of contact may 
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COEFFICIENT OF STATIC FRICTION 


100 200 300 400 500 600 
MOLECULAR WEIGHT 
Fic. 1. Coefficients of static friction showing linear rela- 
tionship with molecular weight for homologous groups. 


exceed several thousand kg/cm? according to 
the pressure at which plastic deformation occurs 
in individual cases. As most of the work done 
against the frictional force is liberated as heat 
in the very small area of real contact, a measure- 
ment of the temperature at the points of contact 
should give information about the actual temper- 
ature of the surface atoms of the rubbing metals. 
Bowden and Ridler*® (1936) have calculated and 
measured such surface temperatures. Calcuation 
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for the simple example of a constantan cylinder 
1 mm in diameter loaded.with 100 g and sliding 
over a mild steel surface with a velocity of 
100 cm per second and having a measured 
value of the coefficient of friction of 0.23 gives 
the following results. 


FRACTION OF SURFACE AREA 
IN RUBBING CONTACT 


CALCULATED 
SURFACE TEMPERATURE 


All of the surface ‘ig > 
1/10 414 
1/100 2372 


The measurements were done by sliding two 
different metals over each other and measuring 
the electromotive force set up by the resulting 
thermocouple.* The temperatures reached de- 
pend upon load, speed, coefficient of friction and 
thermal conductivity of the materials used and 
are in good agreement with the theory. With 
metals of low melting point, the surface temper- 
atures reached correspond to the melting point. 
With less fusible metals surface temperatures 
exceeding 1000°C were observed. Lubricating the 
surfaces under boundary conditions did not prevent 
the surfaces from reaching temperatures of 600°C, 
and it is evident that breakdown or decomposi- 
tion of the boundary film must have occurred. 


Friction Under Boundary Conditions 

The most direct experimental evidence that 
friction of lubricated surfaces under boundary 
conditions is intimately connected with the 
physico-chemical properties of the molecules 
themselves was obtained by Sir William Hardy 
and his co-workers (1922-25) in their classical 
experiments on static friction using as lubricants 
the homologous series of paraffins, alcohols and 
fatty acids. The relationship obtained was 
unexpectedly simple and the friction was found 
to be a function of separate contributions by 
the solid surfaces, the chemical series, and the 
number of carbon atoms in the chain. The 
results were interpreted on the assumption that 
friction took place between monomolecular films 
of lubricant adsorbed on the solid surfaces. 
Figure 1 shows some data obtained by Dr. Givens 

* Blok (1937) (see reference 4) has used the thermocouple 
principle to measure ‘‘flash’’ temperatures on gear teeth 
during meshing and has estimated the highest tempera- 
tures to be of the order of 200-250°C. He also has made a 


general theoretical study of the temperature rise at surfaces 
of actual contact under boundary conditions. 
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the local points of contact, 
and it seems quite impossible 
that under these severe con- 
ditions any monomolecular 
film could stand up unless it 
were self-repairing. Lang- 
muir’ (1934) has described 
some striking experiments to 
illustrate this point. When a 
clean spherical glass rider is 
placed on a clean flat glass 
slide and the slide is tilted, 
the rider will not slide down 
unless the angle is at least 
60°. When a monomolecular 
film of fatty acid is placed on 
the rider, slipping will occur 
at 40° but the rider will only 
slip a few millimeters after 
which the angle must again 
be raised to 60°. When, how- 
ever, the monomolecular film 
is applied to the slide, an 
angle of only 6° to 8° is 
necessary to cause the rider 
to slip off the slide irrespec- 
tive of whether a film also 
had been applied to the rider. 
This is of course due to the 
fact that the rider moves 
continuously to fresh parts 
of the surface, whereas in 
the first experiment, the 
same point of contact is 
maintained with respect to 
the rider and the film applied 
to the rider wears out almost 


Fic. 2. Friction pendulum with photographic record of decrement. 


in these laboratories. Hardy's experiments are 
so convincing that for many years static experi- 


ments have been employed to measure ‘‘oiliness”’ 
or the properties of lubricant under boundary 
condition. The practical value of these data is 
very doubtful, however, as they do not reveal 
any properties of adsorbed films under sliding 
conditions. In fact, we have seen in the preceding 
paragraphs that under conditions of continuous 
sliding very high temperatures are produced at 
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immediately. In another ex- 
periment, the rider moved 
over the slide in recipro- 
cating motion on the same track. If a mono- 
layer of fatty acid was placed on the slide, 
it took 700 strokes to increase the coefficient 
of friction from 0.11 to 0.13 and an addi- 
tional 200 strokes to raise it to 0.24. With a 
film 7 layers thick the coefficient of friction did 
not increase perceptibly even after 4000 strokes. 
Similar experiments were carried out by Dr. 
Givens and the author using a friction pendulum 
pictured in Fig. 2 with either a journal bearing 
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on four balls or a four-ball bearing (Fig. 3).7 
The balls and the journal were covered with 
varying numbers of monolayers of fatty acids 
or related compounds like octadecyl stearate or 
tristearin by the Langmuir-Blodgett method. 
Coefficients of friction could be obtained by 
measuring the decrement of the pendulum using 
the photographic record shown in Fig. 2. In 
agreement with Langmuir it was found that no 
less than seven layers would give a small and 
constant decrement. Seven layers, in fact, 
seemed to wear indefinitely.t 

It was already pointed out that with flooded 
bearings of any type true boundary lubrication 
may be difficult to obtain, and if obtainable the 
coefficient of boundary lubrication should be 
independent of the velocity of the rubbing 
surfaces. Experiments with reciprocating motion 
are not very suitable to investigate this point. 
The first clean-cut case of boundary lubrication 
was obtained by Beare and Bowden’ using a 
flat flywheel with a spherical rider. The fly- 
wheel was set in motion and its deceleration was 
observed by plotting the angular velocity against 
time. With highly polished steel surfaces and 
fluid lubricants such as octyl alcohol, a straight 
line was obtained. Viscous and oily lubricants, 
however, showed always a dependence of their 
coefficients of friction upon velocity, even at 
very low velocities. We have shown at Emeryville 
(Beeck, Givens and Smith’), that even with 
viscous and oily lubricants in the flooded state 
true boundary conditions can be obtained. The 
apparatus used was the Boerlage four-ball top 
shown in Fig. 4. This apparatus in which a single 
ball rotates under variable loads on a support 
formed by three similar balls clamped together in 
an oil cup is particularly adaptable to measure- 
ments of coefficients of friction at low sliding 
velocities of 0 to 1 cm/sec. and high pressures of 
a few thousand kg/cm?* or higher (as calculated 
from the elastic properties of the steel balls). 

A careful analysis of the motion of the rotor 
was made by means of a photoelectric device in 
conjunction with a revolving drum recorder. 


t A friction pendulum was described by Kyropoulos and 
Shobert (reference 5). The four-ball arrangement and the 
photographic recording method are the author's. 

tSee also recent experiments with built-up films by 
Bowden and Leben, reference 6. 
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Instantaneous values for the deceleration were 
obtained by graphical differentiation. The results 
were finally expressed in plotting coefficients of 
friction against velocity. Figure 5 shows some 
typical results. At certain critical and easily 


Fic. 3. Four-ball bearing for friction pendulum. 


reproducible velocities, the curves bend over 
sharply and the coefficient of friction becomes 
constant. Ricinoleic acid, the main constituent 
of castor oil, does not show this effect. Neutral 
or white oils likewise do not show the effect at 
low velocities. It is imparted to these oils by the 
presence of small amounts of such fatty acids 
or other compounds that are readily adsorbed 
in form of well oriented surface films as shown 
by electron diffraction experiments. 

Since built-up films of fatty acids of one or a 
few molecular layers will not produce the effect, 
it has been concluded that the adsorbed layer of 
long chain polar molecules is able to induce a 
wedging effect which allows an early transition 
from the boundary state into a quasi-hydro- 
dynamic state in which lower coefficients of 
friction prevail. Rough surfaces do not show the 
wedging effect as would be expected for two 
reasons. The interlocking of the asperities of 
such surfaces and the extreme pressure and 
temperature conditions thus produced upon 
sliding would make the action of an adsorbed 
layer impossible, and the fact that the degree 
of orientation of an adsorbed layer depends 
greatly on the degree of perfection of the 
surface would make the induction of the wedging 
effect less probable. 
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Fic. 4(a). Four-ball top with photoelectric chronograph. 


We have seen in an earlier paragraph that 
Bowden and collaborators have found rapid 
fluctuation of the surface temperature under 
boundary lubrication conditions. This makes it 
quite probable that with an apparatus capable 
of measuring frictional fluctuations, correlations 
between frictional and temperature fluctuations 
may be found. It is essential for such an appa- 
ratus that the friction measuring device should 
have a high natural frequency and that the 
system which applies the load should have a 
still higher natural frequency. Bowden and 
Leben® constructed such an apparatus and found 
the frictional fluctuation which they called 
“stick-slip.”’ The exact behavior depends on the 
combination of metals, particularly on their 
melting points and, of course, on the lubricant 
if the surfaces are lubricated. In case of dissimilar 
metals the fluctuatious are quite regular. The 
surface of the rider sticks to the moving bottom 
plate (it is, in fact, welded to it) and moves 
with it until, as a result of the increasing force 
which wishes to keep the rider in its original 
position, a sudden and very rapid slip occurs. 
Simultaneously with this slip there is a sudden 
temperature flash which produces the next 
welding, and so on. Surfaces lubricated with 
straight mineral oils behave in the same manner. 
However, long chain fatty acids, for instance, 
may prevent the stick-slip and allow continuous 
sliding to take place. The experiments of 
Bowden and co-workers were carried out at very 
low sliding speeds, and the slider had a low mass. 
It would be wrong to generalize that all boundary 
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friction occurs with the type 
of stick-slip just described. 
Blok* (1940) rightly 
suggested that stick-slip is a 
regime which may or may 
not be entered depending on 
dynamical characteristics of 
the sliding machine parts 
and the relation of static 
and kinetic coefficients of 
friction. In the last analysis 
it is, of course, true that the 
elastic properties of the re- 
straining system may be 
that of the surface irregu- 
larities themselves in which case stick-slip of 
distances of only a few atoms and periods of 
10-* to 10-° second may still occur. 


Boundary Lubrication and W ear 


Although sliding friction and mechanical wear 
are logically closely connected, no significant 
general correlation is found in practice between 
the coefficients of friction and wear with various 
oils with and without addition agents. Under 
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Fic. 4(b). Cross section of four-ball bearing. 


moderate conditions, the coefficient of boundary 
friction varies but little with pressure, tempera- 
ture, and material of both lubricant and surfaces. 
The reason for this lies again in the fact that wear 
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takes place momentarily at isolated spots which 
are small compared to the whole surface, and 
although these isolated spots must exhibit a 
higher coefficient of friction at the moment of 
abrasion, the resulting change for the whole 
lubricated system is too small to affect the 
over-all coefficient of friction unless refinements 
of the type of Bowden’s “‘stick-slip’’ method 
are used. 

Wear is, however, an accumulative effect and 
therefore its measurement is always possible 
even if the rate of wear is very low. Furthermore, 
the engineer is seldom interested in the relatively 
small variations in power loss due to variations 
in the coefficient of friction in machine parts 
which partly or accidentally have to rely on 
boundary lubrication, but he is interested in 
wear. We have, therefore, attempted at Emery- 
ville to investigate wear (Beeck, Givens, Wil- 


COEFFICIENT OF FRICTION 


ANGULAR VELOCITY OF TOP IN R.PS. 


Fic. 5. Coefficients of friction obtained with four-ball 
top. © castor oil; X ricinoleic acid; VY fatty acid fraction 
from castor oil; © 1 percent stearic acid in ricinoleic acid. 


liams® (1940)) under conditions which would 
enable us to link friction and wear measurements 
more closely than has heretofore been possible. 
A motor-driven modification of the four-ball top 
was used as shown in Fig. 6. As wear was 
unmeasurably low at the low speeds at which 
the wedging effect was observed, speeds of several 
hundred revolutions per minute had to be used. 
A successful innovation was the use of steel balls 
which were electroplated with a very thin and 
uniform copper plate, whose wear characteristics 
were found to be similar to those of iron, and 
which were simply run until the track on the 
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Fic. 6. Four-ball wear apparatus. 


upper rotating ball was worn through to the 
steel. This time could be easily measured, and 
the experiments were found to be reproducible 
within 1 percent. Under these conditions long 
chain polar compounds contributed very little 
to wear prevention, but it was found that other 
known wear prevention agents, like tricresyl 
phosphate, showed a wear reduction factor of 
about five. A careful study, especially by Dr. 
Givens, of the chemistry of the whole metal- 
lubricant system led to a complete understanding 
of the mechanism of the action of tricresy] 
phosphate or similarly acting compounds. Their 
action is essentially that of chemically polishing 
the surface and thereby distributing the load 
and lowering the surface temperature. The 
formation of low melting metal phosphide alloys 
plays a major role in the mechanism of this 
action. Here we cannot fail to be reminded of 
the interesting work of Bowden and Hughes 
(1937) who have investigated the mechanism of 
polishing, surface flow, and the formation of the 
Beilby layer. They find that the process of 
polishing is mainly influenced by the relative 
melting point of the polisher and the solid and 
that the relative hardness is comparatively 
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unimportant. Polishing is brought about by 
intense local heating of the surface irregularities 
and subsequent smearing of the softened solid 
over the surface. 

After having gained some understanding of 
the chemical polishing action of phosphorus 
compounds and the action of long chain polar 
compounds, it seems suggestive that chemical 
polishing agents should create just that type of 
surface on which long chain polar compounds 
would act most favorably. This hope has been 
justified beyond expectation. A white oil whose 
wear reduction factor was raised to five by 
addition of 1.5 percent tricresyl phosphate 
showed a wear reduction factor 10 when, in 
addition to the phosphorus compounds, 1 percent 
of oleic acid was added. Factors as high as 17 
have been found with other combinations. This 
self-consistency of results is most gratifying in a 
field where in practice the multiplicity of condi- 
tions may easily obscure the application of fun- 
damental knowledge. 


Surface Finish 


It would indeed be difficult not to say in this 
connection a word about surface finish, although 
the field is so large that it can only be touched 
upon. It is sufficient to say, however, that 
engineers have been well aware of the fact that 
surface finish has an important bearing on wear; 
and the general tendency of the last 10 years 
has been towards better surface finishes of 
lubricated rubbing surfaces, even if such surfaces 
are operating most of the time under hydro- 
dynamic conditions. The importance of this field 
was well demonstrated at a special Summer 
Conference on Friction and Surface Finish at 
M. I. T. June 6-7, 1940 (see proceedings of this 
meeting). It is particularly significant in the light 
of the work on boundary lubrication that the 
mechanical approach from the standpoint of 
surface finish and the physico-chemical approach 
from the standpoint of the addition agents have 
both led to the realization that mechanical wear 
is mainly a function of the quality and main- 
tenance of surface finish. 
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Electron Theory of Thermoelectric Effects 


BY W. V. HOUSTON 


California Institute of Technology, Pasadena, California 


HE development of the theory of solids 

during the past decade has made possible 
an understanding, at least in a qualitative way, 
of the thermoelectric properties of metals.' It is 
now possible to describe, in physical terms, the 
behavior of the electrons that leads to the 
different signs and magnitudes of the Peltier 
effect, the Thomson effect, and the thermoelec- 
tromotive force. At the same time it is possible 
to illustrate the electronic aspect of the two 
types of electromotive force distinguished by 
Bridgman,’ the driving electromotive force and 
the working electromotive force. Since electro- 
motive force itself.is a concept that is frequently 
misunderstood, this illustration may be of value. 


Definitions of the Thermoelectric 
Coefficients 

The thermoelectric effects that we shall con- 
sider are the thermoelectromotive force, the 
Peltier effect, and the Thomson effect. These 
may be defined in connection with the thermo- 
electric circuit shown in Fig. 1. In this figure, 
A and B represent two different metals connected 
together at the ends. The regions around the 
junctions are maintained at temperatures 7, and 
T2, respectively. Let 7, be greater than 7). The 
conductor B is cut so that we may consider the 
open circuit situation; but we may also imagine 
the ends connected by a resistance, or other 
energy-consuming device, through which a cur- 
rent can flow. It will be understood, however, 
that the ends of the conductor on either side of 
the break, and any apparatus connected between 
them, are maintained at the same temperature. 

When the circuit is open there will be a 
difference of potential across the cut in con- 
ductor B. This potential difference is a measure 
of the thermoelectromotive force in the same 
way that the difference of potential between the 
poles of a battery on open circuit is a measure of 
its electromotive force. This electromotive force 
will be designated by 642(71, 72) to indicate 
that it depends upon the two temperatures and 
also upon the two materials. The order of the 
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subscripts indicates that &4(71, 72) is positive 
if a current tends to flow from A to B at the 
warmer junction. It has been emphasized by 
Bridgman that a thermoelectromotive force will 
exist even if the two conductors are merely 
strips cut in different directions from the same 
anisotropic crystal. This fact imposes some im- 
portant restrictions upon the theory. 

If the two parts of B are connected by a 
resistance, a current will flow around the circuit 
in consequence of the thermoelectromotive force. 
This flow of current will be accompanied by an 
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The Thermoelectric Circuit 

Fic. 1. The thermoelectric circuit. A and B are two 
different metals whose junctions are kept at the tem- 
peratures 7, and 7». 
absorption of heat at the warmer junction and a 
liberation of heat at the cooler, superposed on 
the ordinary liberation of the Joule heat. The 
rate of liberation or absorption of this thermo- 
electric heat is proportional to the current. It 
also changes sign with the current, and if a 
current is forced around the circuit in a direction 
opposite to the direction of the thermoelectro- 
motive force, the heat will be absorbed at the 
cooler junction and emitted at the warmer. 
This is a process by which the cooler junction is 
cooled and the warmer junction is warmed by 
the aid of an external agent. The Peltier coeff- 
cient of a junction, P4z(7T), depends upon the 
two metals involved and upon the temperature 
of the junction. It represents the amount of heat 
that must be supplied to the junction to maintain 
its temperature when unit quantity of electricity 
flows from metal A to metal B. 
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In addition to the heat absorbed at the junc- 
tions, there is heat absorption, positive or nega- 
tive, along the lengths of the conductors where 
there exists a temperature gradient. This also 
contributes to the energy necessary to keep the 
current flowing, and is known as the Thomson 
effect. The Thomson coefficient ¢4(7°) represents 
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Fic. 2. The normal distribution of electrons in the absence 
of an electric field and with a uniform temperature. 


the amount of heat absorbed by a conductor of 
metal A at the temperature 7 when unit 
quantity of electricity flows in the direction of 
increasing temperature, and a unit temperature 
difference exists between the ends of the portion 
considered. Both the Thomson heat and the 
Peltier heat are in addition to the Joule heat and 
must be distinguished from it. They are also in 
addition to the thermally conducted heat that 
will always appear in the presence of a tempera- 
ture gradient, and they are distinguished from it 
by being proportional to the current. In the 
theoretical analysis, as well as in experimental 
work, the thermally conducted heat must be 
evaluated and subtracted from the total heat 
energy. 

It is these three phenomena, the thermal 
electromotive force, the Peltier effect, and the 
Thomson effect, that we wish to understand in 
terms of the behavior of the conduction electrons. 
In studying them we must also treat the elec- 
trical conductivity, the contact potential, and 
the thermal conductivity. 

Thermal conductivity is ordinarily divided 
into two parts, one part due to the electrons and 
one part due solely to the ionic interaction. Only 
the electronic conductivity of heat will be of 
importance here, because only it is associated 
with the electric current. 

Statistical Distribution of Electrons in a 
Metal 

The electrons in a metal must be treated 
statistically, since only their temperature is 
given. The electronic situation is described by 
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giving the average number of electrons per unit 
volume in the various states. These states may 
be specified in various ways, but the simplest 
way is to regard them as very similar to the 
energy characteristic states of completely free 
electrons. The states can then be specified by 
three quantum numbers, k., k,, k- that are 
related to the velocities by the equations 


vz=hk,/m*, v,=hk,/m*, v,=hk,/m*. 


The quantity m* is an effective electron mass. 
If the electrons are entirely free m*=m, but the 
effect of the nuclear binding can be roughly 
described by assigning to m* values that are 
somewhat greater or less than m, and _ that 
depend on (k,, k,, k-). The kinetic energy ¢€ of an 
electron in such a state is (A?, 2m*)(k,?+k,7+k-2), 
and the permitted values of k, are integral 
multiples of 2x, L, where L is the length of the 
piece of metal in the x direction. An electron in 
such a state has a definite kinetic energy and a 
definite velocity. It is also assigned a definite 
potential energy that is taken as constant 
throughout a single piece of metal. The specifica- 
tion of the statistical distribution is then the 
specification of the number of electrons having 
certain velocity components, and this carries 
with it the specification of the number with any 
given energy. 

When there is no temperature gradient and no 
electric field, the electrons are distributed over 
the various states so that the average number in 
each state is 


n=2/| exp [(e—e’)/kT ]—1}, (1) 


where & is Boltzmann's constant and T is the 
absolute temperature. ¢’ is a constant energy 
that depends upon the temperature and upon the 
number of electrons per unit volume. For free 
electrons ¢’ decreases when the temperature is 
increased, and this may be expected to be true 
of all good conductors, although the rate of 
decrease will depend upon the nature of the 
material. 

In Fig. 2 the value of ” is plotted as a function 
of k,, when k, and k, are zero. This may be 
regarded as the plot of m along this particular 
diameter of a solid figure in the k,, k,, k. space. 
ach point in such a space represents an electron 
state, and according to the Pauli exclusion 
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principle, not more than two electrons (one with 
each of the two possible directions of spin) can 
occupy such a state. In assigning the electrons, 
they are put into the states of lowest possible 
energy until they are all located. This fills up 
all the available states below a certain maximum 
kinetic energy. For perfectly free electrons, the 
states having the same energy are located on 
the surface of a sphere, and for the absolute zero 
of temperature all of the states within the sphere 
of kinetic energy e’ are filled, and all those of 
greater energy are empty. At higher tempera- 
tures the boundary between the occupied and 
the unoccupied states is less sharply defined, and 
the function n(e) falls less sharply to zero. But 
even at room temperature, the energy spread of 
this surface region is only of the order of 3 
percent of the energy e’, and is entirely too small 
to be shown on the scale of Fig. 2. 

In Fig. 3 the portion of Fig. 2 in the neighbor- 
hood of «=e is very much enlarged in order to 
show the way in which m decreases from its 
maximum value to zero. The value of as a 
function of k,, rather than as a function of €, is 
plotted, although over the very short region 
shown, € and k, are proportional to each other. 


Contact Potential 


One of the things that it is important to 
understand in the discussion of the thermo- 
electric effects, is the relationship between the 
electron distributions in two different metals, 
when these metals are brought into contact. 
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Fic. 3. Enlargement of the electron distribution in the 
neighborhood of «= e’. 


This involves a study of the contact potential 
difference and the understanding of why the 
contact potential difference is not directly con- 
nected with the Peltier effect. In fact a good 
deal of the study of the theory of the thermo- 
electric phenomena consists in understanding 
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why certain things do not produce the effects 
observed. 

Figure 4 illustrates the relationships between 
two different metals when they are brought into 
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Fic. 4. (a) The potential energy of an electron in two 
uncharged pieces of different metals. (b) The potential 
energy of an electron in two pieces of different metals after 
equilibrium has been attained. 


contact. The figures show the energies as a 
function of a single coordinate that is taken 
perpendicular to the surface of contact. In 
region A there is the metal A, and in region B 
there is the.metal B. Region C represents the 
space between them, which may be treated as a 
vacuum. The heavy line represents the potential 
energy of an electron in these three regions, 
respectively. Figure 4(a) shows the situation 
between two uncharged pieces of metal. An 
electron in B has a lower potential energy than 
it has in A, and between them it has a constant 
potential energy, which may be taken as zero. 
The drop in potential energy that accompanies 
the passage of an electron through the surface is 
designated by W, for metal A and Wy, for the 
other. The hatched horizontal lines represent the 
energy of an electron that has a kinetic energy 
equal to €,’ or ex’ in the two metals. The hatching 
on these lines indicates that this is not a sharp 
upper limit for the occupied states; there are 
some electrons with kinetic energies above €’, 
and there are some states whose energies are less 
than e’ that are unoccupied. Nevertheless €’ is a 
perfectly definite kinetic energy that lies some- 
where in the middle of this narrow limiting 
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region. It is, in fact, the kinetic energy part of the 
thermodynamic potential of the electron gas in 
the metal in question. 

Figure 4(b) represents the situation when the 
two metals have been brought into contact and 
separated again, when they are connected at 
some other point, or when equilibrium has been 
set up by some other means. When the two 
metals really form parts of the same statistical 
mechanical system there can be only one value 
of the thermodynamic potential for the whole 
system. This does not mean that there can be 
only one ¢«’ measured as a kinetic energy, but 
that there is an energy, kinetic plus potential, 
that is the same for all parts of the system. Let 
this be &, and let it be measured from a common 
zero of energy. The two hatched lines are now at 
the same level in the two metals. This is brought 
about because a few electrons from metal A have 
gone over to metal B and have raised the poten- 
tial energy in B by enough to bring the tops of 
the two electron distributions together. The 
number of electrons necessary to do this is so 
small that it has no influence whatever on the 
electron densities in the two metals. 

Between the metals there is now an electric 
field that corresponds to the potential difference 
Vpa. In the figure this is drawn as a difference in 
potential energy for electrons. The difference in 
the conventional electrostatic potential will have 
the opposite sign. This potential difference is the 
contact potential difference. If the two metals 
are touched together and then separated, B will 
be negatively and A positively charged, corre- 
sponding to this contact potential difference. 
The difference between the potential energies of 
electrons in the two metals will no longer be 
the same as when they were uncharged, but its 
amount can be determined from the relation- 
ships illustrated in the figure. 

Further consideration of this figure shows that 
the contact potential differences cannot produce 
a current. It is true that the potential energy of 
an electron is less in one metal than in the other, 
but its kinetic energy is correspondingly greater, 
and the equality of the energies ® in the two 
metals is just the condition for statistical equi- 
librium between them. In equilibrium there will 
be no tendency for electrons to flow one way or 
the other. On the other hand, the kinetic energy 
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part of ® is a function of temperature, and the 
temperature variation will be different for the 
two substances. This will produce a temperature 
variation in Vga which might be supposed to 
be connected with the thermoelectric effects. 
But the constancy of ® across the junction is 
the condition for equilibrium at any temperature, 
and there will be no tendency for a current to 
flow one way or the other across the junction. 
For this reason there will be no electromotive 
force that can be associated with the temperature 
dependence of the contact potentials. 


Effect of an Electric Field 


The discussion thus far has referred to the 
distribution of electrons over the various states 
in the absence of an electric field inside the metal 
and in the absence of a temperature gradient. 
The problem of the theory is, however, to find 
the distribution of the electrons under other 
circumstances, when both an electric field and 
a temperature gradient are present. 

In the presence of an electric field along the 
x axis that tends to accelerate the electrons in the 
positive direction, the k, of each electron in- 
creases uniformly and the whole distribution 
moves along the &, axis in the positive direction. 
This uniform motion along the &, axis is the 
representation of a constant acceleration in the 
x direction. Such an acceleration, and such a 
motion along the k, axis, does not continue 
indefinitely. It is opposed by restoring forces 
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Fic. 5. Influence of an electric field on the electron dis- 
tribution. 


due to the collisions of the electrons with the ions 
of the lattice. Eventually there will be an equi- 
librium set up in which the distribution function 
of the electrons will have very much its original 
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shape, but will have its center displaced in the 
direction of positive k,. The situation is indicated 
in Fig. 5 where the light line shows the equi- 
librium distribution in the absence of a field, 
and the heavy line shows the distribution when 
the field is present. 

When the electrons are distributed as shown in 
Fig. 5 there exists a current. The number of 
electrons moving in the positive direction has 
been increased, and the number moving in the 
opposite direction has been decreased. One may 
also consider that an average drift velocity has 
been added to the random motion of each elec- 
tron. This displacement under the influence of an 
electric field is very small. If it is such as to 
correspond to a velocity of 1 cm per sec., it 
represents a tremendous current, some 20,000 
amp. per sq. cm, and such a shift is almost 
infinitesimal! even on the scale of Fig. 5. It repre- 
sents an increase in energy of the fastest electrons 
of about 1.210-' erg or 7.5X1077 electron 
volt. This is very much less than the tempera- 
ture energy k7 at the lowest attainable tempera- 
ture, and to show it at all in Fig. 5, the displace- 
ment has been very much exaggerated. 

The smallness of the change in the distribution 
due to the electric field justifies its treatment as 
a small correction to the equilibrium distribution. 
In Fig. 5 the equilibrium distribution is desig- 
nated by fy and the correction due to the electric 
field by fe=Ee®. Because the shift along the k 
axis is small compared with the distance in which 
the distribution falls to zero, not only the shift 
but also this correction is proportional to E. 
The sum of these two quantities, f=f)+EeQ, 
corresponds to a current, and the magnitude of 
the current is given by the integral 


I= ev.fdr=E f ES. (2) 


dr is the volume element in the space of k,, k,, R:, 
and S is the electrical conductivity. S is a func- 
tion of the temperature and of the metal. The 
integral involving fy vanishes because there is no 
current in the absence of a field. An approximate 
calculation of the current can be made by 
attributing to all of the electrons in fg the speed 
Ye corresponding to the kinetic energy ¢’. For 
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still higher order effects, however, such as the 
thermoelectric effects, it is sometimes necessary 
to take into account the fact that they do not 
all have quite the same velocity. 


Effect of a Temperature Gradient 


To investigate the thermoelectric phenomena 
it is necessary to study the effect of a temperature 
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Fic. 6. Influence of a temperature gradient on the 
electron distribution. The temperature is increasing toward 
the right. 


gradient on the electron distribution. Here again 
the effect is small and can be treated as a small 
correction to the equilibrium distribution. Figure 
6 indicates the nature of the change in this case. 
The right side of the distribution is steeper, and 
the left side is less steep than in the absence of 
the temperature gradient. This can be associated 
qualitatively with the temperature gradient by 
considering that the electrons represented on the 
right are moving toward the right and have 
therefore come from the left. At positions to the 
left of the one in question the temperature is 
lower and the electrons coming from the region 
of lower temperature will reflect that fact. The 
electrons represented by the left part of the dia- 
gram have come from a position of higher tem- 
perature, and reflect that fact in their distribu- 
tion. For the detailed derivation of the shape, 
reference must be made to some treatise on the 
theory of electrons in metals.! 

In this case it is convenient to divide the 
correction to the equilibrium into two parts as 
shown in Fig. 6. The first part fr; has the same 
shape as fz, but in it 2 is multiplied by the 
temperature gradient instead of by the electric 
field. It is also multiplied by the rate of change 
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of ¢’ with respect to the temperature, and there- 
fore fr: is proportional to de’/dx. Combining 
these terms, fri1= — (de’/dx)Q. If we consider at 
first only fri, it is clear that it represents a 
stream of electrons moving toward the high 
temperature end of the conductor. It is, of 
course, impossible to maintain such a state 
unless some means are available to introduce 
electrons at one end and remove them at the 
other. In an ordinary finite wire, the charge will 
accumulate at one end until there is sufficient 
electric field to correct the electron distribution 
in the opposite direction and to stop the current. 
Considering only fr; the current will be 
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e dx 


(3) 


The current will be stopped when E = (1 e)de’. dx. 
Associated with this opposing field there will be 
a potential difference between the two ends of 
the conductor equal to 


2 1 
sv=-f (4) 
e 
The energy corresponding to ¢’ will be the same 
at both ends of the conductor, but the relative 
contributions of kinetic and potential will be 
different. Figure 7 shows the situation in which 
the temperature difference between the two ends 
is correlated with a potential difference. 
Although this first part of the correction to the 
electron distribution will set up a_ potential 
difference between the two ends of a wire, and 
this potential difference will depend upon the 
material because de’ dT is different for different 
materials, such a mechanism will not produce any 
of the ordinary thermoelectric effects. To see 
that this is the case consider Fig. 8. Here it can 
be seen that the potential differences set up in 
the material due to the temperature gradient are 
just balanced by the potential differences that 
are set up at the junctions. The electric field at 
any point is proportional to de’/dx and _ this 
quantity, when integrated around any closed 
circuit will give zero. To the approximation con- 
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sidered here, the quantity ® is constant through- 
out the whole circuit and this constancy guar- 
antees that there will be no electromotive force, 
because it represents equilibrium of the electrons. 

To understand the thermoelectric effects it is 
necessary to consider the next approximation and 
to include the term frez in the distribution func- 
tion. This term is proportional to the temperature 
gradient, dT dx, but it is independent of de’, dx. 
It therefore exists in both metals where there is a 
temperature gradient, but it does not occur at 
the junctions, where ¢’ changes in a region of 
constant temperature. fr2 also contains (e—e’) 7 
as a factor, and this makes it a small quantity 
of higher order than either fg or f71. The complete 
expression for fre is 


and the electric current associated with it is 
1 dT 
| v.frodr= — adr 
T dx 
1 dT 
=—- N, (6) 


where the integral N is a function of the tempera- 
ture and of the nature of the substance involved. 
Both the direction and magnitude of this current 
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Fic. 7. Potential difference set up between the two ends 
of a conductor due to a temperature gradient acting 
through only. 


are very sensitive to the exact value of the 
velocity associated with each electron state and 
also to the exact number of states per unit 
energy range. 

For free electrons, fr2 represents a flow of 
electrons in a direction opposite to that of the 
thermal gradient, and therefore an electric cur- 
rent toward the high temperature end of the 
conductor. If all of the electron states repre- 
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sented a velocity of magnitude v, and the 
number of states per unit energy range were 
constant, fr2 would represent a zero current. 
However, both of these quantities vary con- 
siderably from metal to metal, and thus one can 
account for the wide variety of values for the 
thermoelectric coefficients. 


The Thermal Electromotive Force 


Because of the presence of fr2 and the existence 
of a current associated with it, the potential 
difference between the ends of a conductor in 
which there is a temperature gradient is not that 
given in Eq. (4). The potential difference will be 
either greater or less than this amount, depending 
upon the sign and magnitude of the current 
associated with fre. This means that the quantity 
# will not be a constant throughout the con- 
ductor in spite of its thermodynamic significance. 
This is possible because although # is a constant 
for a system in thermodynamic equilibrium, a 
metal in which there is a temperature gradient is 
not in equilibrium. It is carrying a heat current 
in the effort to set up such an equilibrium. The 
statistical distribution under these circumstances 
is not to be found by investigating the conditions 
for equilibrium, but must be found in other ways. 

The extra potential difference due to the 
presence of fr2 is not compensated by any 
corresponding potential drop across the junction. 
When it is taken into account, the existence of a 
thermoelectromotive force can be understood as 
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Fic. 8. Potential distribution set up in a thermoelectric 
circuit if only fr, were effective. 


resulting from the situation shown in Fig. 9. 
Its magnitude is given by the difference between 
the integrals of the fields in the two metals 
that serve to stop the current flow. 

When all the correction terms are considered, 
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Fic. 9. Potential distribution in the thermoelectric 
circuit when the complete effect of the temperature 
gradient on the electron distribution is considered. 


the total current is given by 


Ide’) 1 dT 
=| — S- N. (7) 


For the steady state on open circuit, when the 
current must vanish, it is necessary that the 
electric field be such as to cancel the other terms. 
This requires a field given by 

Ide’ 1dTN 

k= + , (8) 

edx Tdx 
where both N and S are functions of the tempera- 
ture and of the material. The thermal electro- 
motive force is the integral of the negative of 
this field around the circuit. The first term gives 
no contribution to the integral because ¢’ is a 
single-valued function of position. The second 
term, however, does not vanish. It leads to 
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This represents the integral of the ‘driving 
electromotive force’ of Bridgman. It is located 
in the individual conductors and not in the 
junctions, and it is associated with the fact that 
the effort of the electrons to restore thermal 
equilibrium in the conductors leads them to 
flow along the conductor. This tendency is 
greater in one conductor than the other and so 
produces a resultant electromotive force. For a 
thermal electromotive force to exist in a circuit 
composed of pieces cut in different directions 
from a single crystal it is only necessary that 
N_S be different in different directions. S is 
known to depend upon direction in a crystal and 
it is to be expected that N will also depend on 
direction and in a different way. 
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Driving Electromotive Force 


There has often been a good deal of confusion 
about the significance of the term ‘electromotive 
force.”’ It is properly applied to the line integral, 
around some closed curve, of the vector force 
acting on a unit positive charge, and it is implied 
that such a force is proportional to the charge. 
If this total force, due to all possible sources, is 
designated by Et, then 


(10) 


It is clear that & depends upon the curve about 
which the line integral is taken, since the integral 
between two points is not generally independent 
of the path between them. 

In addition to this use of electromotive force 
as a line integral around a closed curve, it is also 
possible to speak of the electromotive force along 
a specified path between two points, although 
this may be less customary. When this is done, 
the electromotive force between two points close 
tagether, along the shortest path between them, 
becomes equal to the scalar product of Et and the 
distance between the points. This usage may 
lead to a confusion between & and Et, or between 
& and the potential difference between the 
points. 

The total field Et can be divided into two 
parts such that one part is the gradient of a 
potential. Let these two parts be Es and Ea, 
where Es= —gradV. Es is then the field due to 
fixed electric charges, and because it is the 
gradient of a scalar potential, its integral around 
any closed curve will be zero. Then 


= (11) 


The electromotive force about any closed circuit 
is just the integral of Ea around the circuit, and 
the importance of the electromotive force is that 
it separates Ea from Es in this way. The current 
in a circuit is due to Ea only, since Es cannot 
maintain a continuous current. Bridgman* has 
designated this Ea as the source of the driving 
electromotive force of the circuit. 
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Working Electromotive Force 


On the other hand, the electromotive force is 
often regarded as a source of energy, and when 
unit quantity of electricity moves around a 
circuit, the amount of work done on it will be &. 
However, the electrons themselves have no more 
energy after completing the circuit than they had 
at its beginning. Along the path they have trans- 
ferred their energy to the material of the con- 
ductor in various ways, and have transferred to 
it all of the energy gained. The electromotive 
force for the whole circuit, calculated from the 
energy relations, must be the same as the driving 
electromotive force for the whole circuit; but 
there is no reason for believing that they must 
be obtained by integrating the same fields. 
There is no reason for believing that the energy 
given to the electrons is always immediately 
transferred to the material. It may be carried 
along the circuit for a time, and then transferred 
to the material of the conductor at a distant 
point. 

Bridgman defines a ‘‘working electromotive 
force’ in terms of the energy given to a unit 


quantity of electricity when it moves around the 


circuit. Let Ew be a vector of the dimensions of 
force per unit charge and such that Ew-ds is the 
energy transferred from the material body to the 
electrons when one unit of charge moves a 
distance ds. The assumption of such a vector 
implies that the rate of absorption of energy is 
proportional to the current and will be positive 
or negative depending upon the direction of the 
current. In addition to this energy exchange 
there is the Joule heat. The Joule heat represents 
energy transferred from the electrons to the 
material of the circuit to the amount of RJ’, 
where R is the total resistance of the circuit. 
If there are energy transformations proportional 
to other powers of the current they are neglected. 
The total rate of energy transfer to the electrons 
in the whole circuit is then 


This is set equal to zero because a group of 
electrons after completing the circuit have the 
same energy as at the beginning. 
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From Eq. (12) it follows that the integral of 
Ea around the circuit is the same as the integral 
of Ew. The total electromotive force of one kind 
is the same as the total electromotive force of 
the other kind. It does not follow, however, that 
Ea and Ew are the same. In fact, in the thermo- 
electric case under consideration they are quite 
different, as was established by Bridgman from 
thermodynamic arguments. To evaluate Ew on 
the basis of the statistical theory, one must 
consider the flow of energy in the electron gas. 

The heat current, or energy current, carried by 
the electrons is given by 


Woe f vee fet fraddr 


1 de’ 1 dT 
-(z-- K, (13) 


e dx T dx 


where 


L=N+¢S/e. 


The form of this expression shows that an energy 
current will exist in the presence of an electric 
field, a temperature gradient, a change of e’, or 
any combination of these three in which they do 
not neutralize each other. By means of Eq. (7) 
it is possible to eliminate the electric field, and 
to obtain an expression for the heat current in 
terms of the electric current and the tempera- 
ture gradient. The result is 


L 
+ —~—K (14) 
S Tdx\S | 


This gives the heat current entirely in terms of 
fre when there is no temperature gradient, but in 
terms of fr; and fre when no electric current 
exists. 

The quantity that is of importance is not the 
flow of heat energy, but the accumulation of this 
energy at various places. Since the electron 
distribution is based upon the establishment of 
equilibrium due to interaction with the metallic 
ions, a change in the rate of heat flow represents 
an absorption of energy from the metal, or by 
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the metal from the electrons. Let Q be the heat 
energy that accumulates per unit time per unit 
length of the circuit. This reduces to 


Q EI Wo dx 


d dT 
=——]T ( ) + (: ), (15) 
S dx\TS dx\ dx 


where « is the thermal conductivity. The quantity 
Q is made up of two parts. The first part is the 
energy put into the electron gas by the electro- 
static field E, and the second is the energy that 
is left at a given point because of the decrease in 
the energy stream. Strictly, this represents the 
accumulation of energy in the electron gas, but 
because the gas is in such intimate contact with 
the metal ions, and because it is always in 
thermal equilibrium with the metal ions, Q is 
the total deposition of energy at a point in the 
conductor. 

The energy Q is expressed as a sum of three 
terms in the last form of Eq. (15). The first 
term is the ordinary Joule heat, and the last 
term is the accumulation of heat that would be 
attributed to a changing thermal conductivity. 
Neither of these is included in the idea of the 
working electromotive force. The working elec- 
tromotive force is based on the integration of the 
coefficient of J around the circuit. 


d N 
Ey = —T ( ). (16) 
dx\TS 


E,, is different from zero in a homogeneous 
metal in which there exists a temperature 
gradient, because of the presence of 1/7° in 
the derivative, and also because both N and S 
are functions of the temperature. It is also differ- 
ent from zero at the junction between two metals 
because N and S depend upon the metal. The 
energy absorbed by the electron stream under the 
first condition gives rise to the Thomson effect. 
From its definition the Thomson coefficient is 
given by 


dT adTj1N 

dx dx|(TS adTX\S 
Under the second set of conditions the heat 
absorbed gives rise to the Peltier effect. From the 
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definition of the Peltier coefficient it follows 
that 


BdsN 
Pap(T) -{ E,dx = ( 
A A dx = 


A given electric current flowing in metal A 
carries with it a different thermal current than 
the same electric current in metal B. The differ- 
ence must be supplied at the junction and 
constitutes the Peltier heat. On the other hand, 
the thermal current accompanying a_ given 
electric current in a single metal depends upon 
the temperature, so that heat must be added or 
subtracted as the current flows toward regions of 
higher temperature. The exact values of these 
differences are very sensitive to the detailed 
properties of the electron states in the metals, 
and to the exact form of fre. 


The kelvin Relations 


From the forms of Eqs. (17) and (18) it 
follows directly that 


d Pap O4—CR 
( )- ; (19) 
T 


This is one of the Kelvin relations that was 
originally based on a thermodynamic treatment 
of the phenomena. It also follows from Eqs. (9), 
(17) and (18) that 
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which is the other Kelvin relation. It must not 
be concluded from the fact that these relations 
come out of the statistical treatment, that the 
statistical treatment provides a justification for 
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the assumption of reversibility used in the 
thermodynamics. It is merely shown that fo the 
approximation of the statistical treatment the 
Kelvin relations are valid. The whole statistical 
method is an approximation method. The elec- 
tron distribution in the presence of the field and 
the temperature gradient is not obtained exactly, 
but is made up of the distribution without the 
field plus some small corrections. It is gratifying, 
although not at all surprising, that this treatment 
gives results similar to those obtained by Kelvin 
and Bridgman. 

In his discussion of the thermoelectric phe- 
nomena, Bridgman worked out expressions for 
the quantities EF, and Ey by thermodynamic 
reasoning. According to him, 


E,,=o0dT dx (21) 


in the body of a homogeneous substance and the 
integral of E,, through a junction is the Peltier 
heat. This is identical with the result given in 
Eqs. (16) and (17) above. This identity is really 
only a question of the definition of E,,. Bridgman 
also works out the driving electromotive force as 
the integral of E,, and he gives for Ey, the ex- 
pression 


dT( 


dx 


This can be identified with the expression in Eq. 
(8) if ¢, T and (NTS) are taken as approaching 
zero as 7 approaches zero. This is to be expected 
both from Nernst’s third law of thermodynamics 
and from the significance of the integrals N 
and S. It is also necessary that f(7’)(d7T dx) be 
identified with —(1 e)de’ dx. This is not strictly 
a function of the temperature only; but it is a 
function of position in the circuit only, and gives 
a zero integral around the circuit. However, it 
has a value at the junctions as well as in the 
region of a temperature gradient in the homo- 
geneous metal. 

It can be seen from this discussion of the 
electron theory of the thermoelectric effects that 
the quantitative computation of the values of the 
thermoelectric coefficients is not a simple matter. 
It requires a very high precision calculation of 
the energy states for the electrons in a metal, 
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and it requires more precision in this calculation 
than is usually possible at the present time. 
Nevertheless it is possible to understand the 


qualitative features of the phenomena, 


variety of magnitudes and signs and the tem- 
perature dependence, in terms of a picture that 
is identical with that used for the other properties 
of the metals. 
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Frictional Phenomena. I 


This series of articles will be concerned with the present knowledge of viscous and 
related frictional processes, to be considered from a logical, unified standpoint. It will 
consist of three main parts covering, respectively, gases, liquids and solids. Theories, 
experimental data, and technical applications in the fields of mechanics, electricity, 


and chemical technology will be discussed. 


BY ANDREW GEMANT 
The Detroit Edison Company, Detroit, Michigan 


Chapter I. Frictional Processes in General 


N this introductory chapter a general and 

rather formal theory of viscous processes is 
presented that is valid, with certain restrictions, 
for all three states of matter. In particular, the 
fundamental differential equation that can be 
used for the computation of any particular case 
is derived, and is then illustrated in the flow of 
a gas or liquid between two parallel walls. 


1. Definition 


By frictional phenomena are understood those 
surface processes that are involved in any rela- 
tive motion of two bodies in contact with each 
other. This definition can be generalized by 
extending it to continuous media, any two 
neighboring elements of which being considered 
as the two bodies mentioned before. Frictional 
processes will then occur within any continuous 
medium if there exists in it a finite velocity 
gradient. These statements, and some more 
detailed ones to follow in the next section, are 
of a phenomenological nature. An explanation 
of the underlying physical processes for all three 
states of matter will be given in the three 
respective chapters. 

Considering the fact that nearly any motion 
in nature can be described by relative velocities 
of neighboring elements as just mentioned, the 
great generality of the conception of friction 
will easily be understood. In fact, only astro- 
nomical motions, on the one hand, such as the 
rotation of the earth, for instance, and motions 
of atomic constituents such as the stationary 
rotation of electrons around nuclei, on the other, 
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can be reckoned as performances without fric- 
tional accompaniment. But it is almost impos- 
sible to devise a motion in the intermediate 
range devoid of any frictional process. 

If this be true for motions with which pure 
physics is concerned, the more will it hold for 
technical appliances of any kind. Thus, frictional 
processes will have a considerable bearing in the 
realm of the technical sciences, especially in 
applied mechanics, hydraulics, and electricity. 


2. General Characteristics 


The phenomena of friction show a great 
variety, as will be shown in detail in this mono- 
graph, yet it is desirable to point out here in the 
beginning the essential features that can be 
traced in any special type. In order to find this 
common feature, let us consider a rather typical 
case of motion involving friction. A solid body, 
the shape of which is not essential, should be 
assumed to float on water and to be dragged 
along with constant speed on the surface of a 
large water container. 

The details of the velocity distribution of the 
water in the neighborhood of the floating body 
should not be considered. It is known that the 
water is carried along with the body at a speed 
that rapidly diminishes with increasing distance 
from the solid, and that, at a distance of several 
times the average dimension of the body, the 
water will be practically at rest. In any case, 
the water in the neighborhood will have a lesser 
speed than that of the body, and according to 
our definition given above, some frictional 
process will be operative between the two. 
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The chief result of this process will be a 
certain force exerted upon the body which is 
called the frictional force. The direction of this 
force is opposed to that of the speed of the body 
relative to the water. This is a general rule and 
will hold for any frictional process. Viewing 
viscous forces from this standpoint reveals them 
as forces that have to be overcome by the 
external forces in order to maintain a given 
motion. 

According to the principle of relativity, the 
body can be treated as being at rest and the bulk 
of the water as moving in a direction opposite 
to that in which the body moved. Viewing the 
process from this angle, it follows that the body 
must exert a frictional force upon the water. As 
this force must be opposed to the relative 
velocity of the water, it must act in the direction 
of the motion of the body. The magnitude of 
the two forces is, of course, the same. This is 
just a special application of the general law of 
action and reaction. Looking at the viscous 
forces from this angle they are revealed as forces 
that result from a transmission of any motion, 
maintained by external forces, to other parts of 
a body. Both conceptions are instructive for the 
understanding of viscous processes. 

Although the direction of frictional forces is 
determined as stated above, nothing of similar 
generality can be said with regard to the function 
expressing its magnitude. As a matter of fact, 
several of the different mathematical possibilities 
are encountered in nature. The most usual case 
is that the force is a function of the relative 
velocity itself. In this group again the more 
common type of relation is the simple propor- 
tionality between force and velocity. This was 
the case in our above example. However, any 
deviation from this linear law is realized in 
nature, as is the case with certain colloidal 
solutions, for instance. 

The force might also be entirely independent 
of the characteristics of the motion, being just 
a constant. This is the case with some types of 
external solid friction. 

Again, it is possible that the force, although 
entirely independent of the relative velocity, is 
a function of the displacement, reckoned from a 
certain position of equilibrium. Thus the internal 
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frictional force, occurring during vibration in 
solids, is often a linear function of the amplitude 
of vibration. But, whatever the function might 
be, the sign of the force acting on a part is 
always opposite to that of the velocity of the 
part relative to its neighbor. 

There is another general law which is a direct 
consequence of the first one. Let us return to 
our concrete example. In moving the body we 
have to perform a certain work which for the 
time unit will be 


A uf, 


if w=velocity, f=frictional force. This work is, 
according to our previous rule, always positive. 
Whatever the type of motion, A, is positive in 
each of its elementary phases, and_ steadily 
accumulates. There is no motion imaginable, 
not even a purely vibrational one, moving the 
body to and fro through a position of equi- 
librium, in which even a part of the time integral 
JS A;dt should be zero. The reason for this is 
that f remains always opposed to u, whatever 
its magnitude. 

The work A, is totally converted into heat 
energy. This must be the case, since in a vibra- 
tory motion the whole system is exactly the 
same at the end as it was in the beginning. 
An accumulation of any other kind of energy is 
not possible, thus the work lost on the system 
can only turn into heat. In other words all 
frictional phenomena are of the dissipative type. 

In mechanics, dissipative processes are treated 
in contradistinction to conservative ones in 
which latter there is an interconversion of work 
and kinetic or potential energy. 

In thermodynamics frictional processes are 
called irreversible, as distinguished re- 
versible processes, such as the conversion of the 
internal energy of a system into external work 
and latent heat during an infinitely slow iso- 
thermal process. 

3. Deformation Velocities 

After this general review we proceed in 
formulating the equations that allow frictional 
processes to be treated mathematically. 

As just pointed out, the expression repre- 
senting the frictional force as a function of the 
characteristics of the motion is not universal. 
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This fact, unfortunately, prevents a generally 
valid theory from being formulated. There is, 
however, a way out of this difficulty, in con- 
sidering that the most frequent relation is linearity 
between force and relative velocity. It is possible 
to formulate the equations under this assump- 
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Fic. 1. Neighboring elements in continuous media. 


tion, and these equations will be presented in 
the following. 

One must, however, bear in mind that they 
are not generally valid, but mainly restricted 
to gases and liquids. For solids they cannot be 
applied without further comment. The reason 
for this is that molecular mechanisms underlying 
solid friction are entirely different from those 
encountered with gases and liquids. 

Purely formally it is possible—and sometimes 
even useful—to apply these equations to solids 
as well, but those quantities that for gases and 
liquids represent characteristic physical con- 
stants will have only the significance of an 
auxiliary quantity, which is not a constant but 
might be any function of any other variables. 
We are referring here to the constant that will 
be defined later as the so-called viscosity constant 
of gases and liquids and which has a real physical 
sense. If introduced into cases dealing with solid 
friction, it is just a mathematical expression, 
and is some times referred to as the “equivalent 
viscosity.” 

We shall come back to this point more in 
detail in the third part of our monograph, 
dealing with solids, but it is important to 
mention it now in order that the scope of the 
following presentation may be clearly understood. 

It is not our purpose to present a complete 
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theory that is available in handbooks or text- 
books; some essential points only should be 
given and explained, since they will be needed 
when dealing with special problems in the 
following chapters. 

We consider a continuous medium (Fig. 1) 
and fix our attention on two neighboring volume 
elements (1, 2) possessing a velocity relative to 
each other, thus giving rise to frictional forces, 
as stated in Section 1. In other words, we 
imagine a certain surface element o in the 
medium, separating the two volumes, and shall 
define the direction of the frictional force as 
being exerted from 1 to 2 along the surface o. 
The surface might have any orientation with 
regard to the arbitrary coordinate axes x, y, 2, 
but it is essentially sufficient to consider surfaces 
perpendicular to the three axes. If the behavior 
of the forces on these particular surfaces is 
known, then the construction of the forces on 
any other surface element is merely a problem 
of vector and tensor analysis. 

Moreover, as we are mainly interested in 
isotropic bodies in which all directions are 
equivalent, it will be sufficient to consider one 
special surface element only, since the other two 
are perfectly analogous. Thus, we choose the 
surface perpendicular to the x axis and will put 
down expressions for the force acting from 
volume 1 to 2 along the surface unit of o. Instead 
of the force, which is a vector, its three compo- 
nents should be considered. The x component 
is designated by X’,. The prime should indicate 
that the force in question is a frictional one. 
The index x indicates that the surface in question 
is normal to the x axis, and that the force is 
transferred in the direction of the positive x axis; 
that is, from element 1 to element 2. In the same 
way Y’, and 2’, designate the y and z components 
executed upon a plane perpendicular to x. 

All these forces should be linear functions of 
the relative velocities of neighboring particles. 
The velocity components of a particle with 
regard to the coordinate axes should be u, v, w. 
The relative velocities can then be expressed as 
the appropriate functions of the gradients of 
u, v, and w with regard to the coordinate axes. 

If there is a flow in the y direction, for instance, 
then the velocity of element 2 in Fig. 1 relative 
to 1 is given. by (dv/dx)dx, if the distance 
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between the centers of 1 and 2 is dx. Or, if u 
is taken into consideration, the relative velocity 
of 2 is (du/dx)dx. 

There are altogether nine such gradient 
components, describing the velocity distribution 
within the medium, but, as the detailed theory 
shows, there are only six quantities essential for 
calculating the frictional forces. The six quanti- 
ties are: 

xX, Ox; 


02; 
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dx Oy az dz Ox 


ov ou dw 


the expressions x,, x,, etc., being nothing else 
but shorthand writings for the respective gradient 
components. Thus, the nine components occur 
only in the six combinations given by Eq. (1). 
Owing to the symmetry of the expressions, it is 
just as well to write y, for x,, 2, for y. and x, for z,. 

The six expressions (1)—to be exact, the 
quantities: x,, Vy, termed 
the components of the tensor of the deformation 
velocity. The name originates from the theory of 
elastic media, and it will be quite useful to 
illustrate the correctness of this designation. 
Let us consider the term du/dy. A plane ABCD 
of the body (Fig. 2) situated in the xy plane of 
the coordinate system should assume the shape 
A’B'CD after a certain flow in the x direction 
has taken place. The displacement of any point 
should be denoted by a, being a function of y 
as well as of the time, ¢. 

Now the velocity of each point is u=da/dt, 
and the gradient of this velocity as considered 
above is (0/dy)(da/dt). 

Since differentiation with regard to y and ¢ is 
independent, their sequence can be altered, and 
the gradient written as (0/dt)(da/dy). The 
gradient in the bracket is a measure of the 
deformation which in this special case is a 
simple shear. The whole expression is then the 
change of shear with time, or the deformation 
velocity, also called rate of shear in this case. 
Hence, it is shown that the conceptions of 
velocity gradient and deformation velocity are 
equivalent. This result is of importance in con- 
sidering elasto-viscous materials. 

According to our fundamental assumptions, 


the force components X’,,.. are linear functions 
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of the deformation velocity components x,... 
and the question is how many proportionality 
constants will be required. A consideration of 


y 


A B 
A 


Cc D x 


Fic. 2. Diagram of deformation. 


energy relations in conjunction with require- 
ments arising from the isotropy of the body 
shows that generally two constants are sufficient 
to describe the frictional behavior of an isotropic 
medium. Whether both are necessary is another 
question, and the generally adopted view is that 
there must be some simple relation between the 
two, leaving only one characteristic constant 
for an isotropic medium. 

The linear equations to be obtained by means 
of the energy relations just mentioned are the 
following : 


— — +9, +22) 


(2) 


Z' — 122, 
the two constants being 7 and & As mentioned 
before, similar equations are valid for the planes 
normal to the y and z axes. Although the de- 
duction of (2) is omitted here, a kind of inter- 
pretation in two special cases will be given in 
the next section. 

1. Stress Components 

relations 
between » and & should be put forward. The 
sum (x,+¥,+2-) Means, as is known from vector 
analysis, the total expansion of the volume unit 
per second, as is designated by divU, if U is 
the velocity vector at any point. 


Before explaining Eq. (2), the 


Further, one has to remember that in liquids 
and gases there is always present the hydrostatic 
pressure p which acts in the same direction as 
X’,, namely, normal to any surface. Thus, one 
can express the total surface forces as the sum 
of the frictional forces and the hydrostatic 
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pressure, and designate them as X,, Y,, and Z,. 
One has then, 


X,=p—2nx,—£ divU 


3 


For incompressible liquids evidently 
divU =0 
and, from the first Eq. (3): 
p=3(X.+ Y,+2Z,), 


the equation for Y, and Z, being similar to that 
for X,. For incompressible liquids the hydro- 
static pressure proves to be the mean value of 
the three normal force components. 

This rather plausible consequence should now 
be assumed to hold for compressible liquids 
(i.e., gases) as well, an assumption that is in 
keeping with experimental observations. Carry- 
ing out the addition, one sees that the condition 
necessary for this assumption to hold is: 


— 


Whereas, therefore, the first of Eqs. (3) for 
liquids assumes the final form 


p—2nx,, (4) 
that, for gases is 
X,= p—2nx.+3n divU. (5) 


The only constant remaining in the equation 
is now 7; it is characteristic of the medium in 
question and depends only on such variables as 
pressure and temperature. It is called the 
viscosity constant, and will be considered more 
in detail in several of the following chapters. 

As to the interpretation of Eq. (3), let us first 
consider the case wherein of the nine velocity 
gradient components only one, du/dy, is finite 
and positive; all the others are zero. Eq. (3) 
then reduces to 


(6) 


The other force components are zero. Figure 3(a) 
shows at the, bottom the distribution of u with 
regard to the x—y axes, and above, a volume 
element. The particular type of motion, as 
assumed, is maintained by certain external 
forces, acting at the boundaries of the body, and 
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the frictional forces must be such that they 
transmit the particular motion, in this case a 
pure shear, into the interior of the body. This 
is done exactly by the forces parallel to the 
x axis, as given by Eq. (6) (X,). 

Simultaneously, however, there are also forces 
parallel to the y axis, as required by Eq. (6), 
and shown in the figure. The significance of 
this latter is the following. 

X, is in reality a couple, trying to rotate any 
volume element in a clockwise direction. Y,, on 
the other hand, produces a couple in an anti- 
clockwise direction. In every case of a stationary 
flow those two opposite couples will exactly 
compensate each other; the resulting torque is 
zero and, should the volume element rotate, 
which generally is the case, then its momentum 
of rotation must remain constant. 

In the second case let us consider du/dx 
finite and positive, all other components vanish- 
ing. Then Eq. (3), with the value introduced 
for —, becomes 


4 du 2 du 
X,=p-—-n—; (7) 
3 dx 3 Ox 
the other components being zero. 

Figure 3(b) illustrates this case. The velocity 
gradient is shown at the bottom. Above is a 
volume element with the frictional forces 
(X’,...). Since the motion as assumed leads to 
an expansion of the medium, there must be a 
negative pressure component in the x direction 
as required by the first Eq. (7). In addition, 
there are positive pressure components in the y 
and z direction, the magnitude of which is such 
that the sum of these three friction components 
should be zero, or the sum of the total compo- 
nents should be 3p. 

The consideration of these two special cases 
should suffice as an interpretation of Eq. (3), 
since all the others are different combinations 
of these two elementary ones. 

5. The Fundamental Differential Equation 

In dealing with experimental cases involving 
friction another step is necessary. Up to now 
we have formulated the frictional forces as 
acting on surface elements. This is, indeed, 
correct and in keeping with the physical picture 
according to which the force is transmitted from 
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one body to another across their common 
boundary, or from one volume element to 
another in a continuous medium. If, however, 
the velocities of these volume elements have to 
be calculated—and this is one of the most 
frequent problems—then these frictional surface 
forces have to be transformed into volume forces. 
This is a frequent operation in mathematical 
physics, and we will consider how it is done in 
our special problem. 

A cube of volume one, with edges parallel to 
the coordinate axes, should be considered within 
the medium in question. There are generally 
four types of forces acting upon it. 

(a) Let us begin with the frictional ones. The 
six forces acting upon its six surfaces will 
evidently be equivalent to the forces acting on 
the volume. In summing them up, it will be 
useful—following the procedure adopted in the 
previous sections—to separate the vectors in 
question into their three components. Only the 
x components should be dealt with. 

The x components of the frictional forces 
acting on the yz planes are X’, acting on the 
one surface and — X’, acting on the other. Were 
X’, constant throughout the medium, the two 
would just compensate each other. Only in 
case there is a finite gradient of this force, will 
the resultant differ from zero. One of the forces 


a b 
y 
x x 
z 


Fic. 3. Diagram of stresses: (a) du/dy finite; (b) du/dx 
finite. 


is then X’,, the other —(X’,+0X’,/dx), since 
the cube has been taken as of volume one, and 
their sum is —dX’,/dx, acting in the positive 
x direction. 

The same situation holds for the X’, and X’, 
components. The actual force acting on the 
volume is always given by their respective 
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differential quotients. Thus, the total frictional 
force is 
ax’, aX’, ax’, 


Ox oy Oz 


(b) In the same way the hydrostatic pressure 
has to be taken into consideration. Only in 
case of a pressure gradient is a force exerted 
upon that volume element; its magnitude is 
—dp/dx. The forces under (a) and (b) can be 
dealt with together, by replacing X’ with X. 

(c) There might be any external force acting 
on the volume elements, such as gravitation or 
electric forces in case of free electric charges, 
the x component of which should be designated 
as F,. 

(d) For nonstationary processes the velocity 
is not constant with time. The inertia force has 
then to be considered too, having an x component 
—pdu/dt, if p is the density. 

Adding the forces reviewed under (a) to (d), 
the following equation of motion valid for each 
volume element is obtained: 


Ou OX, OX, OX, 
F,—p -( + + )=0 (8) 


Ox oy Oz 


and exactly analogous ones for the y and ¢ 
components. 


6. Energy Relations 


At the end of this introductory chapter the 
main energy relation controlling frictional pro- 
cesses should be given, such as can be derived 
on the basis of the above theory. The same 
limitations hold with regard to solid friction, 
as indicated in Section 3. 

The derivation of the equation will not be 
given; we again content ourselves with having 
it explained. 

In the last section external forces F have 
been introfuced which, together with the hydro- 
static pressure p, can be considered as_ the 
driving forces. They will perform a_ certain 
external work per second, the amount of which 
should be A. This work will be used up in three 
different ways. 

First, it might be transferred into the kinetic 
energy L of the system. This transformation 
will take place during transitory, or initial 
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processes, and the increase per second of the 
kinetic energy should be denoted by AL. This 
quantity is evidently zero for stationary motions. 

Second, it might be stored as the elastic energy 
U of the system, playing a role only in gases 
possessing a finite compressibility and in solids 
with finite elasticity (also elasto-viscous liquids 
see the corresponding chapter). Its increase per 
second is AU. 

Third, a certain amount of heat is developed 
due to friction which, as has been shown, is a 
dissipative process. The heat produced per 
second and per volume unit of the medium 
should be W. 

Then we have the equation 


A=AL+A u+f Wad2, (9) 


the last integral being taken over the total 
volume. 

Our special problem is now to express W as a 
function of the force and deformation velocity 
components, introduced into our general theory. 
This function is the following : 


—(X' yyy 224+X' x, 
+ (10) 


Instead of a rigorous proof which is outside 
the scope of this monograph, an interpretation 
should be given similar to that of Eq. (3) in 
Section 4. 

Take @u/dx as finite, all other velocity 
gradient components vanishing, and consider 
again a cube of unit volume as in Section 5. 
Equation (10) then reduces to 


W= —X’,du/dx, (10a) 


where X’, acting on the two yz planes is given 
by Eq. (2). Since the motion of the cube in 
question is relative to the surroundings, the left 
of the two yz planes (see Fig. 3(b)) should be 
considered for the moment as fixed in space and 
only the right one moving. Its velocity relative 
to the space and in the positive x direction is 
then du/dx. X’, is the force acting upon the 
volume unit on this right surface, being directed 
toward the negative x direction. Its magnitude 
in the positive x direction is then —X’,. The 
work done per second by this force is then, 
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according to elementary mechanics, — X’,du/dx 
(force times path). Since, however, this work 
must be converted entirely into heat (Section 2), 
we have exactly the above equation. 

Again, let us take only du/dy as different 
from zero. Equation (10) becomes 


W=—X’',du/dy. (10b) 


X’,, acting on the x—z planes, is given by (2). 
The lower x—z plane (Fig. 3(a)) should now be 
fixed, the upper moving with a velocity du/dy 
in the positive x axis. The force acting on this 
plane and in this direction is —X’,, its work 
per second — X’,du/dy, as Eq. (10b) requires. 

In summing up over all gradients in the same 
way, one obtains the complete form given by 
Eq. (10). 

It is, of course, possible to substitute Eq. (2) 
into (10), and obtain expressions, which sdme- 
times are more advantageous. Thus, (10b), for 
instance, becomes: 


W=n(du/dy)*, (10c) 


showing that, whatever the sign of the gradient, 
the developed heat is positive. 

It is important to note that, although in some 
cases the actual external forces are restricted to 
certain boundaries of the medium in question, 
yet heat is being developed in each volume 
element of the medium possessing a_ finite 
velocity gradient. Certain cases will be dealt 
with in the next chapters, illustrating this 
behavior. 

7. Flow Between Parallel Walls 

A special case should be considered here, in 
order to give an application of the equations 
deduced. Two parallel walls in a distance of 2h 
(Fig. 4) should form a channel in which a flow 
of the velocity « in the direction x takes place. 
The driving pressure difference per unit length 
should be p. The height of the channel (y 
direction) is unessential. All data below will 
refer to a height of unity. 

Generally we have to assume that u is a 
function of z. Thus, writing down the components 
of the deformation velocities (1), we have only 


x,=du/dz, 


all others being zero. 
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As to the stress components (3), we have 
according to Eqs. (6): 


du 


dz 


X,.=Y,=Z,=p. 


Now we apply the differential equations (8) 
to a volume element of thickness dz as shown in 
the Fig. 4. There are no volume forces F present, 
and in the steady state in which we are interested, 
the inertia forces disappear. Consequently, we 
have for the x component 

d*u dp 


n—_——=0 
dz* dx 


whereas the other two components yield 
Op 
=—=0 


dy Oz 


indicating that a pressure difference acting in 
the x direction is the only one. Denoting —dp/dx 
by p as mentioned above, and integrating, we 
have 

du/dz= — pz/n, (11) 


the integration constant being zero, since, for 
reasons of symmetry, du/dz is zero for s=0. 
For the next integration the boundary condition 
at the wall has to be settled; this question will 
be dealt with in the next chapter. Let us here 
assume the most frequent case, namely, that 
there is no slip at the wall: «=0 for z=h. Then 


u=(p/2n)(h?—=). (11a) 


The velocity distribution along the channel 
thickness is a parabolic function with its maxi- 
mum at the center plane. 

The total volume V' flowing through the unit 
height is given by 
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+h 
v= f udz= (11b) 
3n 


being proportional to the third power of the 
channel thickness. 


2 | 
x 


Fic. 4. Diagrams of flow between parallel walls. 


The heat produced is given by Eq. (10). 
Using (10c) for the volume element dz, we have 


p?2* 
dW=nxdz= dz (11c) 
n 
and 
w= dW =2p*h* /3n, (11d) 


which, as can be seen, yields the relation 


W=pV. (12) 


In other words, the heat developed is equal to 

pressure gradient times flow. The resistance of "2a 
the channel, R, per unit length can be defined 

by means of 


R=p/V, (13) 


thus, giving in this case 


R= 3n/2h'*. (11le) 


The simple case discussed is most instructive 
in clarifying the theory of viscosity. Similar 
cases, which are more important from the 
standpoint of experimental devices, will be 
treated in some of the following chapters. 
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The Theory of the Plastic Properties of Solids.” [IV 
BY FREDERICK SEITZ 
Randal Morgan Laboratory, University of Pennsylvania, Philadelphia, Pennsylvania 
AND 
T. A. READ 


Westinghouse Research Laboratories, East Pittsburgh, Pennsylvania 


Part B. Polyerystals' 
INTRODUCTORY DISCUSSION 


HE properties of polycrystals are influenced 
by two separate factors—first by the 
intrinsic properties of the single crystal con- 
stituents or grains, and second by the restrictions 
neighboring grains exert on one another. If it 
were not for the second influence, all properties 
of polycrystals could be derived by taking some 
simple average of the properties of corresponding 
single crystals for various orientations. This 
simple procedure is possible in a few simple 
cases, of course. For example, the elastic moduli 
and electrical resistivity of pure polycrystals 
may be approximated closely? by taking the 
appropriate averages for single crystals. On the 
other hand, many properties of polycrystals, 
such as shear strength and internal friction, are 
affected by the discontinuities in structure that 
occur in a way that cannot be explained simply 
by treating the system as a set of disoriented 
isolated single crystals. We shall refer to the 
additional factor as grain boundary influence and 
shall attempt to unify present knowledge of this 
factor in the next section. It should be borne in 
mind that by grain boundary influence we mean 
not only the effect arising from the atoms in the 
transition region between two grains, but also 
the effect that neighboring grains exert on one 
another. 
On the whole, the amount of purely scientific 
curiosity that has gone into the experimental 
* The previous installments of this series appeared in 
the February, March, and June, 1941, issues of this journal. 
' The analytical conditions used in mathematical treat- 
ments of practical problems in plasticity will not be sur- 
veyed here. They are discussed in the book by A. Nadai, 
Plasticity (McGraw-Hill Book Company, New York, 1931), 
and that by M. Gensamer, Strength of Metals Under Com- 
rev Stresses (American Society for Metals, Cleveland, 


?See, for example, E. Schmid and W. Boas, Kristall- 
plastizitat (Springer, Berlin, 1936), Section 81. 
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study of polycrystals is far less than that 
expended in the investigation of single crystals. 
This fact will be evident at many points in the 
following pages, for only in a few cases is it 
possible to draw conclusions from available 
experimental work comparable with those that 
may be drawn from the basic experiments on 
single crystals. 


7. GRAIN BOUNDARY INFLUENCE 


a. The grain boundaries as barriers for 
thermal and electrical flow 

Since different grains presumably grow from 
different nuclei, it follows that grain boundaries 
will naturally be the regions where insoluble 
impurities will tend to aggregate. For this 
reason, the grain boundaries may act as barriers 
for thermal and electrical conductivity, particu- 
larly in more impure materials. We shall see 
later that this fact has an influence on the 
internal friction of polycrystals. 
b. The nature of the transition layer be- 

tween grains 

Even if the material of which polycrystals are 
made is very pure, or if all impurities are highly 
soluble, we may expect an abnormal arrangement 
of atoms in the immediate vicinity of the 
boundary between grains. That is, we may 
expect a transition layer of atoms which occupy 
positions resulting from a compromise between 
the forces of the atoms in both grains. These 
atoms will not be so tightly bound as the atoms 
in the interior of the grains and, as a result, 
may be expected to be more mobile at a given 
temperature than interior atoms. This does not 
necessarily mean that slip will occur more easily 
along grain boundaries than within grains for, 
as we have seen in the previous sections, slip is 
ordinarily determined by the ease of formation 
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and motion of dislocations rather than by the 
ease with which individual atoms may move. 
At most we may conclude that the diffusion of a 
given atom would take place more easily*® at 
the grain boundaries than inside of grains. In 
fact, it is possible that, as a result of irregularities 
occurring in the intergranular regions, it would 
be difficult for dislocations to move through 
them. 

The width of the transition region between 
grains in very pure crystals is not known from 
experimental work. We know from other sources, 
however, that the forces between atoms in solids 
are of short range, extending with appreciable 
intensity only over a few atom distances. For 


this reason, it seems safe'to conclude that the 


width of the intergranular region is of the order 
of five This 
conclusion is not valid, of course, in materials 
containing a high percentage of insoluble im- 
purities, for in such cases a large fraction of this 
material may be localized at the grain boundaries. 


interatomic distances at most. 


POLYCRYSTAL 
A 


Fic. 51. Slip in a single crystal of zine bounded by a 
polyerystal. (After Miller.) The heavy contour is the 
shape of the specimen after slip, the dotted line the shape 
before. The line AD marks the boundary between the single 
crystal and polycrystal. Slip has not occurred in the region 
OAD neighboring this boundary. The light lines in this 
region designate the position of the latent slip planes, and 
show the orientation of the basal plane in the entire single 
crystal before slip. Simple slip has occurred in the regions 
OBA and in the part of the specimen lying to the left of 
BD, but the latter region has been bent relative to the 
plane BD. The light lines show the slip bands in this 
region. It may be noted that slip occurs by bending in the 
region OCD, which contains planes previously intersecting 


AD. 


From a study of the behavior of very pure 
tin at temperatures near its melting point, 
Chalmers' has concluded that the grain boundary 
material has a slightly lower melting point than 
the bulk material. This conclusion was drawn 


3 It does not necessarily follow that bulk diffusion occurs 
more rapidly along grain boundaries, for the total amount 
of intergranular material is probably very small in a pure 
polycrystal. 

4B. Chalmers, Proc. Rov. Soc. A175, 100 (1940). 
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from the fact that grains separate along their 
boundaries at temperatures somewhat below the 
melting point. The difference between the 
separation temperature and the true melting 
temperature for any pair of grains turns out to 
be independent of the relative orientation of the 
two crystals and of the amount of impurity, 
providing it does not exceed 0.02 percent. For 
high purity tin the measured temperature 
difference was 0.14°C. These results are appar- 
ently in good accord with the conclusion drawn 
previously that the transition layer of atoms is 
thermodynamically less than 
atoms. 


stable interior 

There is considerable evidence that the bond 
at grain boundaries is very strong at tempera- 
tures not too near the melting point, in spite of 
this lower thermodynamic stability. For example, 
fracture occurs most commonly® through grains 
rather than at: their boundaries in rupture tests 
well below the melting point. On the other hand, 
intercrystalline fracture is common near the 
melting point. At first sight, observations of 
this kind seem to support the view that the 
intergranular material actually is stronger than 
the bulk material at low 
weaker at 


temperatures and 
high temperatures. However, an 
alternative explanation of these facts is as 
follows. At low temperatures the grain boundary 
material is somewhat weaker than the grains 
and rupture starts at grain boundaries. The 
difference in strength is not so great, however, 
that a crack will automatically follow grain 
boundary surfaces regardless of their inclination 
relative to the plane of greatest tensile stress. 
In fact, once started, a crack will occur in the 
plane of greatest tensile stress even if this plane 
cuts through grains. At high temperatures it is 
possible that the greater relative instability of 
grain boundaries increases and a crack will 
follow such boundaries. Indeed, Chalmer’s 
experiments indicate that the grain boundaries 
are extremely weak just below the melting point 
of the bulk material. 

The observation that low temperature fracture 
is principally transcrystalline rather than inter- 


5 See, for example, the discussion in the book by Z. 
Jefferies and R.S. Archer, The Science of Metals (McGraw- 
Hill Book Co., Inc., New York, 1924), p. 67. 
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granular was first explained by Beilby® on the 
assumption that each grain in a polycrystal is 
surrounded by a thick coating of brittle amor- 
phous material having much greater strength 
than the interior of the grains at temperatures 
not too near the melting point. Apparently this 
notion must be almost completely discarded, for 
investigation’? of the electrical behavior of 
evaporated films of many common metals indi- 
cates that crystallization can occur at tempera- 
tures near 100°K. Thus amorphous coatings 
many atom layers thick would be highly unstable 
under ordinary conditions. 

More direct evidence against a thick transition 
layer hypothesis has been given by Chalmers® 
from a study of the critical shearing stresses in 
bicrystals of tin. He investigated cases in which 
the axis of tension was symmetrically disposed 
relative to the two grains in the specimen, so 
that the stresses were the same in corresponding 
crystallographic planes in each grain. Results 
showed that the critical shearing stress increased 
continuously as the relative orientations were 
changed, starting from the case in which they 
Thus he 
concluded that there is no essential discontinuity 
in the nature of the grain boundary in passing 


were the same (1.e., a single crystal). 


from a single crystal to a polyerystal, as would 
follow if the thick transition layer theories were 
correct. The origin of the increase in shearing 
stress with increasing disorientation will be 
discussed below. 


Grain boundaries as the seat of stress 
magnification and dislocations 

If block boundaries in single crystals can 
contain weak spots for stress magnification, as 
experiments on slip and rupture indicate, we 
may expect grain boundaries to have similar 
properties. Thus we may expect grain boundaries 
to be the source of dislocations. Excellent 
evidence for this may be derived from experi- 
ments on internal friction of the type discussed 
in Section 2 Part d. If single crystals and course- 
grained polycrystals of very pure copper are 

®* A survey of the theories of intercrystalline bonding is 
given by E. H. Bucknall, Metals Industry 311, 369, 396 
(1929). See also B. Chalmers, Proc. Roy. Soc. A162, 120 
(1937). 

7R. Suhrmann and G. Barth, Physik. Zeits. 36, 841 


(1935); Zeits. f. Physik 103, 133 (1936); R. Suhrmann and 
W. Berndt, Zeits. f. Physik 115, 17 (1940). 


540 


carefully annealed, both have very low decre- 
ments, of the order of 10 

Now we saw in Section 2 that, in the case of 
single crystals, this type of internal friction is 
very sensitive to mechanical treatment, presum- 
ably because the number of dissipating centers 
(dislocations) is increased as a result. It is found 
that the coarse-grained polycrystals are much 
more sensitive. For example, the internal friction 
of an ordinary sized single crystal of copper is 
scarcely affected when the specimen is dropped 
on a wooden surface from a height of one inch, 
whereas the internal friction of an otherwise 
identical polyerystal is raised by a factor of 
about ten. This result indicates that grain 
boundaries are a ready source of dislocations 
and implies, in turn, that slip nuclei are easily 
formed there. Whether this is due to the fact 
that there are larger numbers of weak spots at 
grain boundaries or to other causes cannot be 
said at present. 


d. Slip interference at grain boundaries 
Neighboring grains exert a strong restricting 
influence on the amount and kind of slip that 
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Fic. 52. Comparison of the stress-strain curves for 
single crystal and polycrystalline specimens of zinc. (After 


Elam.) 


may occur within the grains. A_ significant 
experimental investigation of this effect has been 
carried out by Miller’ on specimens of zinc 
consisting of a large single crystal bounded by a 
polyerystal (Fig. 51). Measurements were made 
on specimens extended at 180°C since slip occurs 
smoothly at this temperature. Miller found 
that those slip planes of the single crystal that 
intersect the polycrystalline region cannot oper- 
ate as freely as those which do not intersect. 
Thus in the typical case illustrated by the figure, 


*R. F. Miller, Trans. A. 1. M. E. 111, 135 (1934). 
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stresses sufficient to cause slip in the uninhibited 
regions of the crystal were not sufficient to cause 
slip in the region OAD immediately bordering 
the polycrystal. Although slip occurred in the 
region OBA which does not border the poly- 
crystal, the deformation was highly restricted 
by the presence of the unstrained region in that 
rotation of the slip planes was not permitted. 
Similarly, slip occurred in the region OCD but 
was limited to the amount allowed by bending 
on the plane BD. Miller states that the total 
slip in this region is somewhat less than that 
occurring in the unrestricted part of the single 
crystal. 

The fact that slip can occur in region OCD 
makes it seem surprising, at first sight, that slip 
is absent in OAD, for if the shearing stress there 
had the same value as in the rest of the specimen, 
we might expect some bending, as in OCD. A 
reasonable explanation of this fact is as follows. 
Initially, the shearing stress is uniform through- 
out all slip planes (exclusive of microscopic 
fluctuations), but after an undetectably small 
strain has occurred in OAD, the stress becomes 
redistributed in such a way that the shearing 
stress is lower in this region. This evidently 
requires that the polycrystalline region exert a 
transverse tensile stress on the material to the 
left of the boundary AD, thereby reducing the 
shear in the slip planes. A condition sufficient 
for this is that slip does not occur in the inter- 
granular surfaces along AD. Naturally, the 
transverse tension across AD will diminish if 
the polycrystalline material becomes deformed, 
but it is clear in any case that the slip occurring 
in OAD and in the polycrystalline region must 
be closely correlated. 

Now if a single crystal grain is favorably 
oriented for slip and is entirely surrounded by 
less favorably oriented grains, as occurs in the 
interior of a polycrystalline specimen, stresses 
similar to those occurring in Miller’s cases will 
be exerted across the boundaries of this grain 
and slip will be permissible only when the stress 
is sufficient to allow all grains to deform. Natu- 
rally, the more numerous and more randomly 
oriented the planes of easy slip are, the more 
probable it will be that one of the easy slip 
planes of an arbitrarily chosen grain is favorably 
oriented for deformation. Thus the stresses 
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Fic. 53. Same as Fig. 52 for aluminum. In this case the 
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needed to produce a given deformation should 
be smaller for cubic crystals, which have many 
easy slip planes, than for hexagonal crystals, 
such as zinc, which have only one. Figures 52 
and 53 bear out this conclusion by showing that 
the difference between the stress-strain curves 
for single and polycrystalline specimens is much 
greater in zinc’ than in aluminum,!® which is a 
face-centered cubic crystal possessing four equiv- 
alent planes of easy slip. 

It is natural to ask at this point whether or not 
the actual stress-strain curves for polycrystals 
can be explained quantitatively by regarding 
them as a cluster of single crystals subject to a 
simple constraint, such as that there shall be no 
relative motion at grain boundaries. This prob- 
lem has not yet been subject to careful analysis 
so that it is not possible at present to say that 
the constraining condition is a simple one. We 
shall return to this topic briefly in the next 
section. 

One of the facts indicating that the strength 
of polycrystalline specimens cannot be derived 
with the use of a simple constraining condition 
is that the Brinell hardness," the form of the 
stress-strain curves, and the breaking strengths 
of polycrystals seem to vary with grain size™ 
in a range in which the grain size is much smaller 
than the dimensions of the specimen and of the 
ball used in the test. For example, Fig. 54 shows 
the variation of the hardness of brass with grain 
size. Results of this kind are very surprising if 
taken at their face value, for if the width of the 
grain boundary region is as small as we have 
supposed, we should expect the mechanical 
properties of polycrystals to be independent of 
grain size as long as the size is small compared 
with the dimensions of the specimen and of the 
device used in measuring the hardness. A possible 
explanation of this discrepancy lies in the fact 

® Taken from C. F. Elam’s book, The Distortion of Metal 
Crystals (Oxford University Press, 1935), p. 53. 

'’°R. Karnop and G. Sachs, Zeits. f. Physik 41, 116 
(1927). 

" The Brinell hardness is determined by measuring the 
size of the indentation produced when a 1-mm ball is 
pressed into the specimen with a specified load. (See A. S. 
M. Handbook (1939), p. 112 et pete 

® Investigations bearing on this point have been carried 
out by the following: G. Masing and M. Polanyi, Zeits. f. 
Physik 28, 169 (1924); W. H. Basset and C. H. Davis, 


Trans. A. I. M. E. 60, 428 (1919); Wood, Phil. Mag. 10, 
1073 (1930). See reference 9, p. 52. 
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that specimens must be treated in different 
manners in order to obtain different grain sizes 
and it is possible that this treatment has a 
profound effect upon the mechanical properties, 
particularly if the materials are not perfectly 
pure. Good evidence for this viewpoint has been 
given by Corson® from a study of very pure 
copper. He found that the tensile strength and 
hardness of coarse-grained ingots of very pure 
copper are practically the same as the corre- 
sponding properties of fine-grained specimens 
and concluded that the weakness ordinarily 
observed in large-grained materials is to be 
attributed to the influence of gaseous impurities. 
His results suggest that the gaseous impurities 
congregate at the grain boundaries and, for 
reasons not yet completely understood, weaken 
the bond between grains. In_ coarse-grained 
specimens the total amount of grain boundary 
area is small compared with that in fine-grained 
materials so that the impurity is more concen- 
trated and has a larger effect. 

It should be noted at this point that the con- 
straints provided by the intergranular material 
possibly are not the only factors contributing to 
variations of stress within polycrystalline ma- 
terials on a scale of dimensions of the order of 
grain size. In this connection, Barrett and 
Levenson have shown that the deformation 
within individual grains in compressed poly- 
crystalline aluminum and iron is not uniform. 
In particular, the relative orientation of different 
parts of the grains changes progressively as 
compression proceeds. For example, in the case 
of aluminum the spread of orientation ranged 
from 7° to 10° for 10 percent compression, from 
15° to 25° for 30 percent compression and from 
35° to 45° for 60 percent compression. This 
Variation in orientation was accompanied by 
the appearance of narrow bands on either side 
of which the orientation was different. These 
deformation bands were also observed when single 
crystals of aluminum were compressed by 50 
percent between carefully lubricated plates. In 
these cases the difference in orientation on either 
side of the band was of the order of 3°. The 
origin of the bands has not vet been given a 


148M. G. Corson, Trans. A. I. M. E. 128, 398 (1938). 
4 C.S,. Barrett and L. H. Levenson, Trans. A. I. M. E. 
137, 112 (1940). 
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satisfactory explanation. It seems possible, how- 
ever, that they are regions such as the line BD 
in Fig. 51 at which bending takes place. If so, 
their appearance indicates that the stresses are 
not uniform throughout individual grains. This 
would not be surprising if the bands were ob- 
served only in the polycrystalline specimens for 
the non-uniformity could then be interpreted as 
being related to the influence of different faces 
of the grain. However, the observation of bands 
in single crystals indicates that the grains them- 
selves are inhomogeneous on a scale of distance 
much larger than that of the spacing between 
slip bands. 


8. Stipe IN POLYCRYSTALS 


It is a natural scientific desire to hope to base 
a general discussion of the behavior of poly- 
crystalline solids in various deforming processes, 
such as tension, compression, rolling and torsion 
on a few primary empirical laws. Present knowl- 
edge of the basic principles of slip in grains was 
outlined in the early sections of this series of 
articles. Although this knowledge cannot be said 
to be complete, inasmuch as the role of such 
factors as deformation bands is not yet under- 
stood, we at least have a qualitative formulation 
of the laws of plastic flow in single crystals. 
Unfortunately, as we saw in the preceding 
section, present knowledge of grain-boundary 
influence is not sufficient to complement the 
work on single crystals and provide us with a 
good foundation for treating polycrystals. In 
spite of this, we shall outline briefly the principal 
experimental facts concerning slip in polycrystals. 
a. Crystalline orientation” 

When polycrystalline specimens in which the 
grains are randomly oriented are subject to 
uniform but directional deformation several im- 
portant and probably closely related changes 
occur. In the first place, the orientations of the 
grains become altered. This reorientation fre- 
quently occurs in such a way that the resultant 
distribution of grains is no longer random, but 
exhibits preferred arrangements. The degree of 
preferred orientation is never nearly as sharp as 
in a single crystal and usually depends both 
upon the material of which the specimen is made 


% This topic is discussed in reference 9, Chapter V. 
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Fic. 54. Variation of Brinell hardness of brass with grain 
size. (After Elam.) 


and upon the method used to produce the 
deformation. For example, when face-centered 
cubic metals are placed in tension, there is a 
tendency for both the (111) and (100) directions 
to become oriented parallel to the direction of 
tension. In aluminum the existing evidence indi- 
cates that most of the grains become oriented so 
that the (111) direction is along the axis of 
tension, but both orientations occur in other 
metals having the same structure. Similarly, in 
compression the (110) direction of face-centered 
cubic crystals tends to become aligned normal to 
the plane of compression. In some cases the final 
orientations of the grains in polycrystals appear 
to be similar to those occurring in single crystals 
that have been subject to the same deformation, 
but more often they are not. This fact indicates 
once again the importance of the grain boundary 
influence in determining the actual stresses 
within the grains. For example, in tension tests 
with single crystals of aluminum, the (110) direc- 
tion tends to become parallel to the axis of 
tension, whereas the (111) direction tends to 
become aligned in this direction in polycrystals. 

The reorientations occurring during rolling are 
of particular practical importance since much 
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useful metal is fabricated by this procedure. 
The stresses occurring in this process are more 
‘complicated than in simple tension or com- 
pression tests, because the surface of the rolled 
specimen is subject to normal stresses, and to 
shearing stresses both in the direction of rolling 
and at right angles to this direction. Although 
strong preferred orientations ordinarily are pro- 
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Fic. 55. The upper continuous curve is the computed 
stress-strain relation found by Taylor for aluminum using 
single crystal data. The crosses and plus signs represent 
experimental values. The lower curve is a representative 
stress-strain curve of a single crystal of the same material. 


duced by rolling, the results are greatly de- 
pendent upon the conditions of rolling and 
cannot be summarized in terms of a few concise 
laws which could be conveniently discussed here. 
Apparently rolled orientations can sometimes be 
duplicated in compression tests in which the 
specimen is constrained so as to spread in only 
one direction in the plane of compression. 

One of the interesting effects of extensive 
deformation on polycrystalline materials is the 
development of fiber structure. This structure 
results from the elongation of individual grains 
in the direction of the deformation, which gives 
them the appearance of fibers of woven fabrics. 
During the elongation neighboring grains appear 
to maintain close adhesion at their boundaries, 
for the rupture strength of the materials is not 
reduced except in extreme cases of working. 
This fact indicates once again the great strength 
of the intergranular bond. Moreover, it shows 
that such elongation takes place under great 
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restrictions. It can be shown that the develop- 
ment of fiber structure would not be possible in 
hexagonal crystals if slip in the basal plane were 
the only possible means of producing deforma- 
tion. Actually twinning seems to provide the 
additional degree of freedom in such crystals. 
It is also possible that inner-crystalline deforma- 
tion of the type observed by Barrett and 
Levenson is an important means of deformation 
in many materials. 
b. Work hardening 

Polycrystals harden with cold work, as may 
be seen from the stress-strain curves shown in 
Figs. 52 and 53. Undoubtedly work hardening of 
the type observed in single crystals is a very 
significant factor in polycrystals for the latter, 
like the former, may be softened by heating. 
It is possible, however, that additional factors 
related to the presence of grain boundaries 
enter into the hardening of polycrystals. For 
example, there are indications'® that complete 
recovery in polycrystals requires recrystallization 
of the specimen. In contrast with this, we have 
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Fic. 56. Creep curve for polycrystalline cadmium. (After 
Andrade and Chalmers.) 


seen in Section 2 Part g that single crystals can 
recover completely far below their recrystalliza- 
tion temperature. It is not unlikely, of course, 
that the temperature accompanying recrystal- 
lization and not the change in structure is the 
important factor in determining softening. Since 
polycrystals recrystallize at far lower tempera- 
tures than single crystals, it is plausible to expect 
that polycrystals happen to recrystallize in the 
temperature range effective for softening, which 
is the same for single crystals and polycrystals. 

A highly interesting treatment of the stress- 
strain curve for polycrystalline aluminum has 


'® See reference 9, p. 161. 
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been given by Taylor." Briefly, he computed the 
stress-strain curve on the basis of the following 
assumptions : 

(a) The strain is uniform in all grains of the 
material. This evidently contains implicitly the 
assumption that no slip occurs at grain bound- 
aries, but it is obviously far more stringent. 

(b) The deformation takes place in each grain 
by slip on the octahedral planes, as he and Elam 
had found for single crystals, and in such a way 
that the energy required is a minimum. The 
second part of this assumption is essentially a 
statement of the principle of least work. 

On the basis of these assumptions, Taylor was 
able to compute the stress-strain curve for a 
polyecrystal of randomly oriented grains with the 
use of curves for single crystals. The upper con- 
tinuous curve of Fig. 55 is the theoretically 
derived one and is accompanied by experimental 
points. It may be seen that the computed curve 
agrees substantially with these points. 

Unfortunately, the study carried on by Barrett 
and Levenson on the orientation of grains 
in “homogeneously”? compressed polycrystalline 
specimens of aluminum indicates that assump- 
tion (a) is not correct and thus detracts from the 
significance of the agreement found by Taylor. 
Moreover, Barrett and Levenson found that the 
average change in orientation of the grains 
during compression was frequently different from 
that predicted on the basis of Taylor's treatment 
and thus opens assumption (b) to question. The 
constructive implications of Barrett and Leven- 
son’s work have not yet been determined, but it 
seems definite that Taylor's theory must be 
modified in some essential respects, at least in 
the case of compression. 


ec. [Impurity hardening 


The hardening influence of soluble impurities 
in polyerystals seems to be qualitatively the 
same as in single crystals and indicates that the 
hardening is basically intragranular. 


9. TWINNING 


Although twinning has been extensively in- 
vestigated in polycrystals, none of this work 
seems to have been carried out with the view of 


7G. 1. Taylor, J. Inst. Metals 62, 307 (1938). 
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Fic. 57. Dependence of the creep rate upon stress for a 
bicrystal of tin. (After Chalmers.) 


finding differences between the behavior of poly- 
crystals and single crystals. 


10. CREEP" 


Just as in the case of single crystals, the 
elongation versus time curves for polyerystals 
may be divided into a transient and a steady- 
state part. Such a curve is shown in Fig. 56 for 
polycrystalline cadmium!® and the possibility of 
a division of this type is clearly indicated. The 
dependence of the steady-state creep rate upon 
both temperature and stress has been investi- 
gated by a number of workers’ and seems to be 


'’ For helpful discussions of this subject we are indebted 
to Drs. W. Kauzmann of the Westinghouse Research 
Laboratories and A. Lawson of the University of Penn- 
svlvania. 

'"E. N. da C. Andrade and B. Chalmers, Proc. Roy. 
Soc. A138, 348 (1932). 

20 An analysis of available experimental data is given by 
W. Kauzmann, Trans. A. I. M. E. (February meeting 
1941), Tech. Pub. No. 1301, 
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governed by a function of the form 


R=A exp (—(e'—be)/kT), (1) 


where A and ¢’ are constants, o is the applied 
stress and } is dependent upon temperature, but 
not upon stress. It is to be emphasized that this 
law is derived from a study of results corre- 
sponding to the practical range of stress, which 
is of the order of magnitude of 1 kg/mm?. This 
is from ten to a hundred times larger than the 
range of stress employed in the creep experiments 
on single crystals of tin described in Section 3. 

Figure 57 shows the creep rate as a function 
of stress for a bicrystal of tin, as determined by 
Chalmers.* This result indicates that Eq. (1) 
fails in this case for stresses of the order of 100 
g mm’. Figure 57 should be compared with 
Fig. 44 for single crystals, for it indicates that 
there is a distinct difference in the creep behavior 
of single crystals and polycrystals at low stresses. 
In particular, there seems to be a limiting stress 
for creep in the case of the bicrystal whereas 
there is none in the case of single crystals. 

Extension versus time curves for polycrystalline 
specimens of tin stressed in the range of Fig. 57 
are shown™ in Fig. 58. It may be seen that the 
steady-state rate is practically zero for stresses 
below 200 ¢ ‘mm? and becomes finite in the range 
from 200 to 300 g/mm*. A similar effect has been 
observed by Chalmers™ in polycrystalline lead 
but it has not been studied in other metals. 

Values of the constants appearing in Eq. (1) 
as determined by curve fitting of experimental 
results are given in Table IX. The values of } 
are found to vary with temperature in the 
manner B exp (a7) in which B and a are con- 
stants; the range of b shown in the table gives the 
variations over a several hundred degree tem- 
perature range. 

On general grounds we might expect two 
sources of creep in polycrystalline materials: 
First, a creep analogous to that occurring in 
single crystals and arising from relative motion 
within grains and second, creep at grain bound- 
aries. Slip corresponding to the second type of 
motion does not seem to be observed near room 
temperature in any of the metals that have been 


21 B. Chalmers, Proc. Roy. Soc. A156, 427 (1936). 
2 B. Chalmers, J. Inst. Metals 61, 103 (1937). 
23 B. Chalmers, Proc. Phys. Soc. 47, 352 (1935). 
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studied up to the present time. However, if it 
did exist we should expect to observe it just in 
the range in which creep ordinarily is measured. 
If this type of creep did occur, it presumably 
could be separated from inner-granular creep by 
investigating the dependence of creep upon 
grain size, for we should expect it to play a 
larger role in fine-grained specimens having a 
comparatively large amount of grain-boundary 
surface than in coarse-grained specimens. Ac- 
tually there does not seem to be any conclusive 
evidence indicating that there are two contribu- 
tions to creep. It will be seen in Section 12 that 
there is indirect evidence that intergranular flow 
occurs in zine in the range from room tempera- 
ture to 100°C; however, in lieu of more definite 
evidence, we shall discuss creep on the basis of 
dislocation theory. 

We saw in Section 3 that two simple mecha- 
nisms of creep may be considered in discussing 
creep on the basis of dislocations. First, we may 
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EXTENSION 200 G/MM2 
RECOVERY 300 G 
RECOVERY 200 
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Fic. 58. Upper two curves: extension versus time curves 
for specimens of tin stressed at 300 g/mm? and 200 g/mm?. 
The lower two curves give the subsequent contraction 
when the stress is released. (See Section 13.) (After 
Chalmers.) 


consider creep as caused by the relatively slow 
migration of dislocations, in which case 
dS/di= PrN, (2) 


where N is the density of moving dislocation 
lines, v is their average velocity of motion and \ 
is the slip distance associated with the passage 
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of one dislocation. In. the steady state, N is a 
temperature-dependent constant so that the tem- 
perature dependence of creep rate is determined 
by the dependence of both v and N. 

It was also pointed out in Section 3 that there 
is evidence that the activation energy for motion 
of dislocations is small in a weil-annealed speci- 
men, for which N is small. The same conclusion 
cannot be expected to hold during a practical 
creep test, since N may then be much larger. 
Thus it is difficult to analyze the constants 
appearing in Eq. (1) in terms of the quantities 
appearing in (2). However, we shall attempt such 
an analysis in the following way. 

We shall assume that the rate at which dis- 
locations are generated is given by an equation 


TABLE LX. Empirically determined creep data for metals. 
(After Kauzmann.) (The values of ¢' are given in ev, those of 
A in units of sec.~' for specimens of unit length. b 1s given in 
units of A’, 


LOGw A b 
Pb 0.17 -8 108 
0.1—0.35 —Tto —5 104 
Su 0.23 —4 5.108—5 « 104 


Brass (60 Cu 40 Zn) 0.6 —4 10?9—5x« 10 
Steel 1.2 —3 104 


of the form 
R,= (N, exp (€,—d,o) kT), (3) 


where .V, is the number of generating centers, 
and that the rate of annihilation is 


R,=(N_ ra) exp RT), (4) 


where N is the number of dislocations present. 
Then in the equilibrium state 


N=N,(ta/t,) exp (—(Ae—Aao)/kT), (5) 


where 
Ae=€,—€., Ad=a,—d,. (6) 


In addition, we shall assume that 


kT kT (7) 
Ta 


where d is the distance between neighboring 
equilibrium positions for the dislocation, ra is 
the relaxation time for passage of a dislocation 
from one equilibrium position to another, and 
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TABLE X. Tensile strengths of annealed pure polycrystals. 
(The values are given in units of kg/mm?.) 


TENSILI TENSILI 
STRENGTH STRENGTH 
METAL Oss. Catc. | METAL Oss. CaLc. 

Al 9.2 372 Pb 1.1 261 
Au 12.0 782 Pd 1 
Cs 6.4 Pt 125 
Cu 22.6 885 Sn 1.4 357 
Mg 9.2 Zn 13 580 
Ni 32 


the exponential coefficients give, respectively, 
the influence of temperature and stress on the 
probability that the dislocation will jump. For 
small stresses, the second function presumably 
should be a hyperbolic sine of the same argument, 
but we shall consider the practical range of 
stresses. Substituting (5) and (7) into (2), we 
obtain 


dS hd Ta 
=—N,— exp (—(e’—ba)/kT), (8) 
dt Td Tg 


where 


=€,—€aten, b=a,—Aatay. (9) 


Now tr, and +r, presumably are of the same 
order of magnitude, since they are associated 
with inverse processes. Moreover, \d is of the 
order of magnitude 10~-" cm?, so that we obtain 
by comparison of (1) and (8) 


(N,/ra)10-8=A. (10) 


We are at loss for a precise estimate of N,, but 
it seems reasonable to assume that this quantity 
is of the same order of magnitude as the number 
of blocks per unit cross section, namely 10%. 
The values of ry required to account for the 
observed values of A vary between 10 and 10™4 
sec. for all of the cases listed in Table LX. The 
larger of these values of ry agrees in order of 
magnitude with that obtained in Section 3 from 
a similar analysis of the rate of creep of single 
crystals of tin on the basis of Eq. (2). We note 
again that if Eq. (2) is correct, it seems to imply 
that the natural relaxation time for motion of 
dislocations is very long compared with an atomic 
oscillational period. 

A second possible simple mechanism for creep 
is based on the notion that dislocations move 
very rapidly a distance L from the point at 
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which they are generated and then become stuck. 
The corresponding equation was discussed in 
Section 3 and is 


dS dt. (11) 


As we have seen in Section 3, the purely transient 
creep observed in single crystals of tin, zinc and 
cadmium at low stresses can be explained more 
suitably on the basis of Eq. (2) than on the 
basis of (11), but it is possible that Eq. (11) is 
valid in the practical range of creep. 
Substituting into (10) from Eq. (3), we obtain 


dS N, 
=\L— exp (—(€,—a,0) kT). (12) 
dt Te 


Thus in the present case 
A=10*LN,/r,. 


If L is of the order of magnitude 10~' cm, as is 
assumed in Taylor's theory of hardening, we 
obtain 

A=10-"N,/r,. 


In the cases listed in Table IX, N, +7, varies 
between 10° and 10*. If our previous estimate of 
N,, namely 10%, is valid, we obtain the result 
that r, varies between 10~' and 10* sec., which 
seems to be much too large to be reasonable. 
In other words, the parameter values derived 
from the use of Eq. (2) seem to be more reason- 
able than those derived from (11). It must be 
admitted, however, that the evidence in favor of 
Eq. (2) is far from satisfying. At most, we may 
say that the theory of dislocations seems to 
provide a rough semi-quantitative means of 
correlating creep data. 

Comparing Eqs. (1) and (6), we see that the 
quantity 6 in (1) is composite if Eq. (2) is 
correct. It is interesting to note that the experi- 
mental values of b are so large that be is of the 
order of 1 ev for stresses near 1 kg, mm?*. This 
implies that either a, or a, in Eq. (9) may be 
large. In the first case, it would be concluded 
that there is very large stress magnification near 
the slip nuclei and in the second that dislocations 
have lengths of the order of 10-° cm. It is not 
easy, however, to understand the strong de- 
pendence of the experimental values of ) upon 
temperature. Kauzmann*® has suggested that 
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this dependence is related to an increase with 
temperature of the lengths of either slip nuclei 
or dislocations. Additional evidence for such a 
variation possibly is furnished by the difference 
between the values of + we have obtained above 
from comparisons of Eqs. (2) and (11) with Eq. 
(1) and the value we might expect to obtain, 
namely 10-" sec. According to Kauzmann’s sug- 
gestion the increase in the length of dislocations 
with temperature gives rise to an increase of the 
activation energy with temperature. If the creep 
data were analyzed on the basis of Eq. (1), 
taking into account a temperature dependence 
of e’, a value of A larger by a factor of 10'° could 
be obtained. In this case r would come out to be 
of.the order of magnitude of 10~'® sec. 


11. RupTURE IN POLYCRYSTALS 


It was pointed out in Section 7 that rupture 
is transcrystalline except at temperatures just 
below the melting point. This indicates that 
grain boundaries are not sources of great weak- 
ness at sufficiently low temperatures and we 
might expect to correlate the breaking strengths 
of polycrystals with those of single crystals. 

Such a correlation is made difficult for several 
reasons. In the first place, stresses within grains 
are not the same as the applied stress because of 
grain boundary influence. In the second place, 
specimens of the same material prepared in 
different ways will extend differently during a 
rupture test and hence will usually have different 
orientations when fracture occurs. Fortunately, 
the second factor may be eliminated to an 
appreciable extent by comparing specimens which 
already possess the maximum degree of preferred 
orientation of the type that would be produced 
during the rupture test. 

Now we saw in Section 5 that each crystallo- 
graphic plane of a single crystal appears to 
possess a characteristic tension stress at which it 
will rupture. Moreover, evidence was cited which 
indicates that this stress is independent of 
previous deformation. Reasoning from this we 
should expect the rupture stress of a polycrystal- 
line metal to be independent of deformation as 
long as the deformation does not produce widely 
different orientations. This expectation seems to 
be borne out in the case of copper. Figure 59 
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shows™ stress-strain curves for specimens of 
copper wire prepared from the same material. 
Prior to the test the wires were drawn to different 
diameters, during which process comparable 
degrees of preferred orientation presumably were 
attained. The specimens were then subjected to 
tension stress and the strain was measured in 
terms of the contraction in cross section. The 
» specimens began to neck at the stress corre- 
sponding to point 2, and fracture occurred at 
the points 3. It may be seen that the breaking 
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Fic. 59. Stress-strain curves for variously worked speci- 
mens of polycrystalline copper. It may be seen that the 
rupture stresses are closely the same in each of the cases. 
The strain is measured in terms of the contraction in cross- 
sectional area. 


stresses are nearly the same in all four cases, 
which indicates that the previous strain had 
little effect on the breaking stress. 


J. vy. Mollendorf and J. Czochralski, Zeits. Verein. 
deuts. Ing. 54, 931 (1913). 
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Just as in the case of slip, there seems to be 
evidence that the grain size has an important 
influence on the breaking strength even when the 
crystal size is small compared with the dimen- 
sions of the specimen. As was pointed out in 
Section 7, this result is inexplicable unless it is 
assumed that the different treatment required 
to produce various grain sizes results in different 
distribution of impurities at grain boundaries. 
This interpretation of the effect is supported by a: 
the experiments of Corson, discussed in Section 7. he 

comparison — of 


experimental — breaking 
strengths and those computed from Polanyi's 
equation with the use of experimental values of 
the surface tension is given in Table X. It may 
be seen that the experimental values are smaller 
than the theoretical ones by a factor of the order 
of a hundred or more, just as in the case of 
single crystals. There is a slight indication that . 
polycrystals are somewhat stronger than single 

crystals, but this effect may be related to the 

redistribution of stresses resulting from the in- . 
fluence of grain boundary restrictions. 


12. INTERNAL FRICTION 


The principal source of internal friction in non- 
ferromagnetic single crystals seems to be inti- 
mately associated with plasticity, as we have 
seen on Part d of Section 2. The same source is 
sometimes important in polycrystals, but, as 
was first pointed out by Zener,” is often com- 
pletely masked by another. The second source 
has been extensively studied by Zener and his 
co-workers and arises from 
effect in the following way. 


the thermoelastic 


When a crystal is stressed suddenly its tem- 


perature changes, heat being generated or ab- 
sorbed, depending on the sign of the dilatation. 
The change in temperature is proportional to the 
stress and reverses its sign if the sign of the 


stress is reversed. Since the stresses in a poly- 
crystal are not uniform but vary from grain to 
grain, it follows that the temperature resulting 
from the thermoelastic effect will vary from 
point to point. Now if the stresses vary slowly 
compared with the time required for heat to flow 

%C. Zener, Phys. Rev. 52, 230 (1937); 53, 90 (1938); 
W. Otis and R. Nuckolls, Phys. Rev. 53, 100 (1938); R. H. 


Randall, F. C. Rose and C. Zener, Phys. Rev. 56, 343 
(1939); C. Zener, Proc. Phys. Soc. 52, 152 (1940). 
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between regions of different temperature, the 
temperature will remain constant and the trans- 
fer of mechanical energy into heat, and vice 
versa, takes place reversibly. In this case the net 
amount of heat generated during a complete 
cycle is zero and there is no contribution to the 
internal friction. Similarly, if the variations in 
stress take place so rapidly that no heat flows 
during a stress cycle, the process is said to occur 
adiabatically since each region of the material 
behaves as if it were thermally isolated. It 
follows that the net conversion of elastic energy 
into heat during a complete cycle also is zero in 
this case. However, if the frequency or the 
thermal conductivity are such that the process is 
neither isothermal nor adiabatic, the conversion 
will not be reversible and there will be a finite 
contribution to the internal friction from the 
thermoelastic effect. 

Zener has pointed out that if the average 
applied stress is uniform throughout the speci- 
men, the critical distances for thermal diffusion 
in a polyerystal should be the dimensions of the 
grain, for stresses should vary relatively abruptly 
from one grain to another. Reasoning from this 
on the basis of dimensional analysis, he suggested 
that the internal friction should depend on the 
dimensionless quantity vd*/D. Here d is the 
grain diameter, v is the frequency and D is the 
thermal diffusion constant defined by the equa- 
tion 
thermal conductivity 


D= -, 
(specific heat) (density) 


In other words, for a given material the rela- 
tive value of the internal friction, should be a 
universal function g(vd? D). Figure 60 shows 
schematically the manner in which this function 
should depend on its argument. It has a peak 
near the value unity and decreases to zero on 
either side. This type of dependence of internal 
friction has been verified in brass. 

If the average stress is not uniform throughout 
a specimen, as during the transverse vibration of 
a reed, there will be an additional contribution 
to the thermoelastic internal friction resulting 
from the variations in temperature across the 
specimen. This contribution is determined by a 
function similar to that shown in Fig. 60, the 
argument being vL?/D in this case, where L is 
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the distance between the points of maximum 
temperature difference, which is the thickness of 
the specimen in the case of the vibrating reed. 
The thermoelastic internal friction was _ first 
discovered and studied systematically in inhomo- 
geneously stressed systems of this type. 

The factors contributing to the absolute value 
of the thermoelastic internal friction have been 
investigated by Zener, but will not be discussed 
here since they would lead us too far afield. 
This work shows that in general we should expect 
the internal friction arising from stress variations 
between grains to be smaller the more nearly 
elastically isotropic the material is. 

We might expect to find peaks of the type 
shown in Fig. 60 in single crystals as a result of 
stress variations either from domain to domain 
or between larger regions. The first of these 
cases corresponds to an effective grain size of 
10~* cm and would produce a peak in the mega- 
cycle range of frequency. This type of internal 
friction actually seems to be negligible in Read's 
experiments. 

It should be mentioned that the value of D 
will be dependent on the ease with which heat 
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Fic. 60. Relative value of the internal friction as a 
function of the dimensionless variable yd?/D. (After Zener.) 
The function has its peak in the region where the argument 
is unity and decreases on either side. 


flows across grain boundaries and is not neces- 
sarily an average of the values for single crystals. 
Actually, the influence of grain boundaries has 
not been studied. 

The separation of the internal friction of 
plastic origin from the thermoelastic contribution 
can be carried out in several ways. For example, 
if measurements are carried out for frequencies 
extending into the adiabatic or isothermal regions, 
the residual value of plastic origin will appear as 
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an additive constant on a function of the form 
of Fig. 60. The range of frequency needed for 
this evidently will depend upon d and D. Thus 
in the case of three specimens of brass of different 
grain size,?® a residual decrement of 0.15.1075 
was determined by extrapolating to the region 
of low frequencies. 

Barnes and Zener®’ carried out a_ similar 
separation at fixed frequency in zinc (impurities 
0.5 percent). Measurements were made on speci- 
mens of four grain sizes for various temperatures 
in the range from 0°C to 100°C. The values of 
vd? D at the frequency used were 0.008, 0.008, 
0.02 and 0.8 for the four specimens. In the first 
three cases the internal friction varied with 
temperature in the manner A exp (—e k7), the 
constant e being the same in all cases and the 
constant A varying inversely as the square of 
the grain size. In the fourth case, in which the 
thermoelastic internal friction should be near its 
maximum value, the internal friction was ex- 
pressible in the form 


A=C+A,exp (—€ (1) 


in which C was independent of temperature and 
equal to 44-10~°, « was the same as in the other 
three cases, and the ratio of A, to the values of A 
occurring in the other cases varied inversely as 
the grain sizes. The investigators conclude that 
the temperature-dependent term had its origin 
in the plasticity, whereas the constant in (1) 
corresponds to the thermoelastic dissipation. 
They also conclude from the dependence of the 
coefficient of the Boltzmann factor upon grain 
size that the source of the temperature-dependent 
internal friction lies on the grain boundaries, 
rather than in the interior of the grains. 


13. SECONDARY PLASTIC EFFECTS 


In addition to the primary plastic effects, slip, 
creep, twinning and rupture, there are a number 
of interesting secondary effects which, like 
fatigue, undoubtedly find their origin in the 
primary effects. Some of these are observed only 
in polycrystals; others are observed in single 
crystals as well. We shall discuss both here. 


26°C. Zener and H. Randall, Trans. A. I. M. E. 137, 41 
(1940). 


27 A. H. Barnes and C. H. Zener, Phys. Rev. 58, 87 
(1940). 
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a. Creep recovery 

If a specimen of polyerystalline material is 
extended by creep, it is found to contract slowly 
with continually decreasing rate when the applied 
stresses are removed. The lower two curves in 
Fig. 58 show this effect for two specimens of 
polycrystalline tin. In these cases the ordinate 
gives the decrease in length as a function of 
time after the specimens which yielded the upper 
two curves were unloaded. The rate of contrac- 
tion is rapid at first and then decreases to zero. 
This effect did not occur in the single crystals 
studied by Chalmers, which indicates that it is 
primarily a polycrystalline one. The effect is also 
observed after practical creep tests. 

This type of creep recovery apparently can be 
given a satisfactory explanation on the basis of 
a principle first proposed by Masing.** He pointed 
out that the type of extension occurring in 
different grains of a polycrystal during any 
deformation will be different because of differ- 
ences in orientation. Some grains will deform 
almost entirely as a result of slip whereas others 
will deform elastically. For small deformations, 
such as those dealt with in Chalmers’ experi- 
ments, we may expect a continuous range 
between grains which are deformed almost en- 
tirely elastically and those which have deformed 
almost entirely as a result of plastic flow. Con- 
sequently the actual stress in the interior of the 
metal will vary from a maximum value in those 
grains which have undergone. the least plastic 
deformation to a minimum value in those which 
have undergone the most. Thus when the ex- 
ternally applied load is removed and the average 
stress reduced to zero, some of the grains will be 
stressed in tension and some in compression, the 
latter being those which had previously tlowed. 
If the temperature is sufficiently high, these 
stresses will cause creep and the specimen will 
be observed to contract in the manner of Fig. 58. 
b. The Bauschinger effect 

If the previous interpretation of creep recovery 
is correct, we might expect that jusce after a 
polycrystalline specimen has been unloaded, 
following an elongation or compression, there are 
forces present in some of the grains which would 


28 (G. Masing, Wiss. Siemens Konzern 3, 231 (1924); 4, 
74, 244 (1925); 5, 135, 142 (1926). 
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aid in producing a deformation of opposite sign. 
Thus we might expect the critical stress required 
to produce plastic compression following a small 
elongation to be lower than the critical com- 
pression stress in the virgin state. An effect of 
this kind was first discovered in polycrystals by 
Bauschinger*® and was explained qualitatively in 
the basis of Masing’s principle, discussed in 
part a of this section. Figure 61 shows*’ the 
Bauschinger effect in fine-grained brass. Figure 
61(a) is the stress-strain curve obtained during 
an extension of a virgin specimen of the material 
by a few percent. The critical shearing stress is 
about 10 kg mm*. Essentially the same value 
would have been obtained if the virgin specimen 
were compressed rather than extended. Figure 
61(b) shows the stress-strain curve for com- 
pression following the extension. It may be seen 
that the critical shearing stress for this deforma- 
tion is very small for the range of compression 
corresponding to the original extension. Figure 
61(c) shows the type of stress-strain curve ob- 
tained if the specimen is further extended (with- 
out compression) after the original extension, 
and exhibits the usual effects of work hardening. 
The last figure shows the stress-strain curve 
obtained in compression after annealing the 
originally extended specimen at 150° for several 
hours. The Bauschinger effect has disappeared, 


2” J]. Bauschinger, Ziviling. 27, 289 (1881). 
8° G. Sachs and H. Shoji, Zeits. f. Physik 45, 776 (1927). 
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Fic. 61. Illustration of the 
Bauschinger effect in fine-grained 
polycrystalline brass. (a) shows the 
extension versus stress curve for the 
virgin material. (b) is the stress- 
extension diagram obtained after 
previous extension. (c) is the same 
curve obtained after previous com- 
pression. It may be seen that the 
elastic range is very small, illustrat- 
ing the Bauschinger effect. (d) is 
the stress-strain curve similar to (c) 
obtained when the specimen is an- 
nealed at relatively low tempera- 
tures after compression. The hard- 
ness remains but the Bauschinger 
effect is gone. 


but the hardness remains, showing that the 
factors contributing to the Bauschinger effect are 
not immediately connected with the hardening. 

A quantitative development of Masing’s theory 
of the Bauschinger effect was carried out by 
Heyn,*' who showed that curves of the type of 
Fig. 61(b) could be obtained for polycrystals in 
a reasonable manner on the basis of the theory. 

It is unfortunate for the simple theory that a 
well-defined Bauschinger effect has been ob- 
served by Sachs and Shoji*® in single crystals of 
brass. A typical example of their results is shown 
in Fig. 62. In this the ordinate is stress and the 
abscissa is strain, the positive and negative 
directions corresponding, respectively, to tension 
and compression. The curve obtained in a stress- 
strain cycle starting with a virgin specimen is 
somewhat reminiscent of a magnetic hysteresis 
curve. The segment D, of the entire curve 
represents the stress-strain curve obtained in the 
first extension and rises sharply in the elastic 
region, as in Fig. 61(a). When the tensile stress 
was relieved, the elastic contraction to the 
point A occurred. Following this, the specimen 
was compressed and the curve D2 was obtained. 
It is evident that the decrease from the point A 
is much more gradual than the original rise from 
the origin, showing the Bauschinger effect. If this 
part of the curve were inverted about the point 


Heyn, Festband Katser-Wilhelm Gesellschaft 
(1921), 131. 
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A, it would bear the same relation to D,; that 
Fig. 61(b) bears to Fig. 61(a). The main differ- 
ence between the polycrystalline single 
crystal specimens appears to lie in the fact that 
the latter are somewhat softer, as may be seen 
by comparing the ordinates in Figs. 61 and 62. 
The curve Ds; represents another stress-strain 
curve for extension following the compression Do». 

A summary of the experimental facts obtained 
from this tvpe of investigation on single crystals 
of brass is as follows: 

1. The primary stress-strain curves obtained 
in tension or compression for virgin materials are 
nearly identical. 

2. The Bauschinger effect is observed in the 
opposite deformation following the primary one, 
regardless of the sign of the primary one. We 
shall call the deformation following the primary 
the secondary one. The secondary curves are 
very nearly identical if the primary extensions or 
compressions are identical. 

3. Secondary curves, and curves such as D; 
obtained by later deformation, are very nearly 
alike as long as the over-all extension and com- 
pression remains fixed in each cycle. The prin- 
cipal difference lies in the fact that a gradual 
work hardening takes place; that is, the curves 
rise to larger absolute values of the abscissa. 

4. If at any period in a sequence of cycles, a 
deformation larger than the preceding is carried 
out, the Bauschinger effect will be more clearly 
marked in the deformation immediately fol- 
lowing. 

5. Annealing for two hours in the range from 
250°C to 400°C removes the Bauschinger effect 
for the next deformation, but does not appre- 
ciably affect the work hardening. 

6. When correction is made for the slightly 
greater hardness of polycrystals, the Bauschinger 
effect in single crystals of brass may be said to be 
fully as well defined as in polycrystals. 

At first sight, the existence of a well-developed 
Bauschinger effect in single crystals seems to be 
in contradiction with Masing’s theory. Actually, 
there are two ways of interpreting the results: 

(a) In the first place it is possible that the 
single crystals of brass on which the measure- 
ments were made actually are far less homo- 
geneous than the best single crystals of mon- 
atomic substances. Thus it is possible that the 
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F1G. 62. Stress-strain diagram for single crystal of brass 
obtained during a cyclical deformation process. Negative 


strains correspond to compressions, positive to clongations. 
D, corresponds to the virgin material. 


deformations are not homogeneous and_ that 
parts of the crystal undergo almost purely elastic 
deformations. This possibility could be tested by 
more extensive investigations of single crystals 
of other materials. Apparently, pure monatomic 
crystals are not well suited for tests of this 
type because of their great softness. It should be 
added that Sachs and Shoji made measurements 
on single crystals of aluminum-copper alloys, but 
found the effect smaller than in brass and did 
not extend the investigation. At the present 
time, this interpretation seems to be the most 
plausible one. 

(b) In the second place, it is possible that 
there are two contributing causes for the Bausch- 
inger effect in  polyerystalline materials—one 
being the effect suggested by Masing, and the 
other an effect that occurs even in single crystals. 
Since Heyn’'s work indicates that Masing’s theory 
is sufficient to explain at least the magnitude of 
the effect in polycrystals, we must conclude that 
if two co-existing effects occur they are of com- 
parable importance. 
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Fic. 63. The elastic after effect in a compressed single 
crystal of brass. 


It has recently been suggested by Masing 
that the theory of dislocations may contain 
features usable for an explanation of the Bausch- 
inger effect in single crystals. In particular, he 
has proposed the possibility that dislocations 
which have moved into the crystal as a result 
of the action of a given stress move backward 
under a lower reverse stress than is needed to 
make them move forward. According to Masing, 
this unidirectional hardening is a consequence of 
the fact that the newly formed dislocations will 
encounter dislocations of opposite sign, generated 
on the opposite sides of the particular block in 
which they are moving, if they continue to pro- 
ceed in the direction in which they were started. 
They will not encounter these dislocations if their 
direction of motion is reversed. 

This picture seems difficult to justify on the 
basis of the theory of dislocations developed here 
for we should either expect the positive and 
negative dislocations in the same block to attract 
and lower the shearing stress for motion in the 
same direction, or we should expect them to 
become intermixed in the manner suggested by 
Taylor and form a lattice that possesses equal 
rigidity in either direction. The following, how- 
ever, is an alternate possibility. In discussing 
resoftening of work-hardened crystals, it was 
pointed out that newly formed dislocations near 
surfaces are subject to a strong attractive 
“image” force which drops off as the first power 
of the distance from the surface. It is readily 
seen that this force operates in just such a 
direction as to contribute to the Bauschinger 
effect. Evidently a quantitative estimate of the 
magnitude of this contribution is necessary be- 
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fore we can draw a definite conclusion as to 
whether the mechanism is a likely one. 

The fact that the Bauschinger effect can be 
completely removed by an annealing process 
that does not produce complete resoftening seems 
to indicate strongly that dislocations are not 
primarily responsible for the effect. 


c. The elastic after effect 

Another interesting secondary effect observed 
commonly in polycrystals and in some single 
crystals, such as the brass specimens used by 
Sachs and Shoji, is the elastic after effect illus- 
trated in Fig. 63. It is found that the stress- 
strain curve obtained when a specimen is un- 
loaded and then reloaded, following an extension, 
is not a single-valued function but has the form 
shown in the figure. 

This effect may readily be explained in a 
qualitative manner on the basis of Masing’s 
principle of inhomogenous strain. We may postu- 
late that when the specimen is first unloaded the 
contraction is entirely elastic. As soon as the 
elastic forces in the regions that have undergone 
the most extensive plastic deformation are re- 
leased, the remaining regions will exert on them 
stresses of such sign as to reverse the original 
deformation. These stresses increase as the un- 
loading proceeds and cause a small amount of 
irreversible plastic flow, analogous to the flow 
occurring during creep recovery (Fig. 58). In fact, 
it seems likely that such creep recovery actually 
would have been observed on the single crystal 
specimens of brass used in obtaining Fig. 63. 

Evidently, the questions that arise concerning 
Masing’s theory and the existence of the elastic 
after effect in single crystals of brass are basically 
similar to those discussed in connection with the 
Bauschinger effect in single crystals. 

If we grant that all three of the secondary 
plastic effects presented in this section have 
similar origin, we may conclude from the fact 
that single crystals of tin do not show appreciable 
creep recovery that they would not show a 
Bauschinger effect or elastic after effect com- 
parable with that observed in brass. From this 
viewpoint, then, it follows that the brass crystals 
on which Sachs and Shoji made their observa- 
tions were not as homogeneous as the single 
crystals of tin studied by Chalmers. 
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University of Michigan Summer Session 


In addition to a number of degree courses in physics, 
the University of Michigan Summer Session, 1941, is 
presenting a Symposium in Theoretical Physics extending 
from June 30 to August 22. Lectures in this symposium 
will be given by Professor Wolfgang Pauli, Institute for 
Advanced Study, Princeton University; Assistant Pro- 
fessor Frederick Seitz, University of Pennsylvania; Dr. 
Julian Schwinger, University of California; and Assistant 
Professor Victor F. Weisskopf, University of Rochester. 


* 


Dayton C. Miller Research Fellowship 


In memory of Dr. Dayton C. Miller, who was con- 
nected with Case School of Applied Science for more than 
fifty years prior to his death last February, the Dayton C. 
Miller Research Fellowship in Physics has been established, 
according to an announcement received from Robert S. 
Shankland, present Head of the Physics Department at 
Case. “‘We feel that this fellowship will do much to pre- 
serve the traditions of research in pure science that Dr. 
Miller established here,’’ Dr. Shankland states. 

Dr. Miller’s will made provision for a fellowship award 
to be available at a later date, but funds set aside by his 
friends and former students have made possible its es- 
tablishment this year. The first recipient of the award, to 
be made to a graduate assistant in physics during his 
second year of work towards an advanced degree, is P. J. 
Rice. Mr. Rice plans to continue his research studies on 
the electron microscope. 


* 


Symposium on Light and Lighting 


“Light and Lighting” will be the topic of a symposium 
to be sponsored by the Ohio Section of the American 
Physical Society at the October, 1941, meeting at Cleve- 
land, Ohio, according to a preliminary announcement 
made by Richard H. Howe, Secretary of the Section. 


* 


Members of Scientific Mission Return 


Dr. Warren Weaver, Director of the Natural Sciences of 
the Rockefeller Foundation and Edward J. Poitras, of the 
California Institute of Technology, returned to New York 
by clipper on April 28, having spent five weeks in England 
as members of the scientific mission appointed by President 
Roosevelt from the membership of the National Defense 
Committee. Dr. Weaver, Mr. Poitras, Professor Louis B. 
Slichter, of the Massachusetts Institute of Technology, and 
Dean John T. Tate, of the University of Minnesota, were 
together in London throughout the nine-hour bombing 
raid of April 16 and 17—the most severe London has 
experienced. 

Science 
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Flight Similitude Chamber Aids Meteorologists 


Capable of simulating temperature and pressure condi- 
tions of altitudes up to thirty miles, a flight similitude 
chamber is being used at the University of Chicago to 
train radiosonde operators, urgently needed in the Govern- 
ment’s long-range weather forecasting program. Dr. 
Michael Ference, Jr., instructor in physics and a member of 
the staff of the University of Chicago's Institute of Meteor- 
ology, designed the “stratosphere box,” which is one of 
three in the United States. 

The University of Chicago’s Institute of Meteorology is 
headed by Dr. Carl G. A. Rossby. Dr. Rossby is largely 
responsible for the development of five-day forecasts 
through his work with “‘three-dimensional” airflow charts. 
Thirty army, navy, air corps and weather bureau students 
are being trained at the Institute for government service. 


* 


Officers Elected for the Ohio A. P. S. Section 


The fifth meeting of the Ohio Section of the American 
Physical Society was held in joint session with the Ohio 
Academy of Science at Case School of Applied Science and 
Western Reserve University, Cleveland, Ohio, Saturday, 
May 10. Officers of the Ohio Section elected for the year 
1941-1942 were: W. E. Forsythe, Chairman, General 
Electric Company, Nela Park; Paul B. Taylor, Vice 
Chairman, Aircraft Radio Laboratory, Wright Field; 
Richard H. Howe, Secretary-Treasurer, Denison University. 


* 


Radio Amateur Wins Award for Radio School 


Marshall H. Ensor, high school teacher and radio 
amateur of Olathe, Kansas, received the fourth annual 
William S. Paley Amateur Radio Award recently for his 
outstanding service to the nation in training thousands in 
the radio code and radio fundamentals over a ten-year 
period, during which he broadcast regular lessons in those 
subjects from his own transmitter to thousands of “air 
students.”’ Paul W. Kesten, Vice President of the Columbia 
Broadcasting System, awarded the trophy on behalf of its 
donor, William S. Paley, CBS President. 


* 


British Scientists to Work in America 


Sir Lawrence Bragg, Cavendish Professor of Experi- 
mental Physics at the University of Cambridge, and Dr. 
Charles Darwin, Director of the British National Physical 
Laboratory, recently left England for Canada. Sir Lawrence 
will work in Ottawa at the National Research Buildings 
in cooperation with Professor Fowler as liaison officer 
with Canadian scientific men. Dr. Darwin will work in 
Washington. 

Dr. Horace S. Uhler, Professor of Physics at Yale 
University is retiring from active service after being con- 
nected with the University for thirty-four years, according 
to a recent note in Science. 
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Louis E. Levy Medal Awarded 


“Combined Stress Experiments on a Nickel-Chrome 
Molybdenum Steel,”” a paper published in the August, 
1940, issue of the Journal of the Franklin Institute, has 
won the Louis E. Levy Medal for its authors, John M. 
Lessels and Charles W. MacGregor. The co-recipients of 
the annual award are both Associate Professors in the 
Department of Mechanical Engineering at the Massachu- 
setts Institute of Technology. 

The Levy Medal, presented‘annually by the Committee 
on Science and the Arts of the Franklin Institute, is 
awarded “to the author of a paper of especial merit, 
published in the Journal of the Franklin Institute, preference 
being given to one describing the author's experimental 
and theoretical researches in a subject of fundamental 
importance.” 


* 


Bulletin of Mathematical Biophysics 


The September, 1941, issue of the Bulletin of Mathe- 
matical Biophysics, edited by N. Rashevsky and published 
by the University of Chicago Press, will contain the follow- 
ing articles: 


Electrical Charges and Potentials in Cells Resulting from Metabolism 
of Electrolytes, Ronert R. WILLIAMSON 


Note on the Mathematical Biophysics of Temporal Sequences of 
Stimuli, N. RASHEVSKY 

A Note on the Nature of Correlations Between Different Character- 
istics of Organisms, N. RASHEVSKY 


Mathematical Biophysics of the Galvanic Skin Response, Ciype H. 
CoomBs 


A Theory of Steady-State Activity in Nerve-Fiber Networks: II. 
The Simple Circuit, ALston S. HOUSEHOLDER 


* 


National Academy of Sciences Elects 


Among those elected to membership by the National 
Academy of Sciences during its April meeting held in 
Washington, are the following physicists: Charles Christian 
Lauritsen, Professor of Physics, California Institute of 
Technology, Pasadena; Alfred Lee Loomis, Director, The 
Loomis Laboratories, Tuxedo Park, New York; J. Robert 
Oppenheimer, Professor of Physics, University of Cali- 
fornia, and California Institute of Technology; and Joseph 
Slepian, Research Engineer, Westinghouse Electric and 
Manufacturing Company, Pittsburgh, Pennsylvania. 

This spring meeting of the Academy was also the occasion 
of the presentation of the Henry Draper Medal to Dr. 
Robert Williams Wood, Research Professor of Experi- 
mental Physics at the Johns Hopkins University, ‘in 
recognition of his contributions to astronomical physics.”’ 


* 


Five Research Fellows Appointed 


Five new Westinghouse Research Fellows have received 
appointments enabling them to carry on their studies at the 
Westinghouse Research Laboratories in East Pittsburgh, 
it was announced recently by Dr. E. U. Condon, Associate 
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Director of the Laboratories. The appointees are: Dr. 
Daniel Alpert, Stanford University; Dr. John W. Coltman, 
University of Illinois; Dr. Thomas W. Dakin, Harvard 
University; Dr. Theodore Holstein, New York University; 
and Dr. James S. Koehler, University of Pennsylvania. 

Comprising the fourth group selected under the Westing- 
house Research Fellowship plan inaugurated in 1938, the 
appointees will carry out original research on scientific 
problems of their own choosing in an attempt to discover 
new fundamental facts in the realm of pure science. 


Calendar of Meetings 

July 

23-25 American Society of Civil Engineers, San Diego, California 

August 

27-29 American Institute of Electrical Engineers, Yellowstone Na- 
tional Park 

September 

2-6 American Mathematical Society, Chicago, Illinois 

8-12 American Chemical Society, Atlantic City, New Jersey 

8-10 American Astronomical Society, Williams Bay, Wisconsin 

22-26 American Association for the Advancement of Science, Chicago, 
Illinois 

25-26 Society of Automotive Engineers, Millwaukee, Wisconsin 

October 

6-10 National Safety Council, Chicago, Ilinois 

12-15 American Society of Mechanical Engineers, Louisville, Ken- 
tucky 

13-16 Society of Motion Picture Engineers, New York, New York 

15-17 American Society of Civil Engineers, Chicago, Illinois 

20-24 American Society for Metals, Philadelphia, Pennsylvania 

20-24 American Welding Society, Philadelphia, Pennsylvania 


New Book 


Industrial Research Laboratories of the United 
States 
Seventh Edition. Pp. 372. National Research Council, 


Washington, December, 1940 (National 
Council Bulletin No. 104). Price $3.50. 


The revised edition of this directory contains informa- 


Research 


tion concerning 2264 industrial concerns which maintain 
one or more research laboratories. This list, with certain 
subtractions, includes 469 more concerns than were listed 
in the sixth edition which was issued.in December, 1938. 

The directory gives, in most cases, the following items 
for the laboratories listed: address of the laboratory; name 
of the president of the corporation; name of the research 
director; classification and number of members of the 
research staff; research activities; and cooperative research 
and laboratory publications when the company sponsors 
these enterprises. 

There are several indexes to the directory which cover 
the geographical distribution of laboratories, the personnel, 
and the subject of the research activities. The personnel 


index includes the names of about 6900 research workers. 
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Temperature Correction Methods in Calorimetry 


ALLEN KING AND HORACE GROVER 
Rensselaer Polytechnic Institute, Troy, New York 


(Received March 4, 1941) 


There are two fundamentally different methods of 
correcting calorimetric data for heat losses to surroundings. 
The first method, applicable only to experiments in which 
the supply of heat to the calorimeter can be neglected after 
a certain time, uses the afterperiod cooling curve, together 
with a similar curve from data taken before the heat supply 
started, to evaluate appropriate correction terms. This 
problem is discussed from a general point of view and 
correction formulas previously reported by Regnault and 
Pfaundler, Schultes and Niibel, Roth, and others are 
shown to be variations of a general scheme. A second 
method is needed in such experiments as the measurement 
of specific heat when it is not safe to neglect the heat 
supplied by the hot body even after a very long time. In 
this case, there is no afterperiod cooling curve and addi- 
tional assumptions concerning the heat exchange between 
the calorimeter and the hot body within it must afford the 


LTHOUGH extraneous thermal exchanges 
in calorimetric experiments may be con- 
trolled, usually they cannot be neglected. For 
this reason, numerous methods have been de- 
vised to obtain an approximate correction from 
actual measurements.' Some of these methods 
are erroneous; others are unnecessarily elaborate. 
Often methods suitable for one type of experi- 
ment have been applied to another type for which 
they are inadequate. We shall consider here the 
general problem of correction for thermal ex- 
change and discuss the conditions appropriate to 
the several cases that occur in practice. 
Calorimeter data usually consist of readings of 
temperature as a function of time.* At the outset, 
the time rate of change of these calorimetric 
temperatures will be assumed to consist of two 
terms. One, the “heating function,’’ depends in 
general upon the temperature difference between 
the source of the heat to be measured and the 
calorimeter fluid; and the other, the ‘“‘cooling 


1A discussion of adiabatic methods in calorimetry is 
outside the scope of this paper. 

2In the following account it is assumed that thermo- 
metric lags have been allowed for. 
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needed information. The assumptions originally suggested 
by Rowland are re-examined. It is shown that previous 
investigators (Rowland and Hoare) deduced inadequate 
correction methods from these assumptions. A_ better 
method is developed and a simple approximation of it is 
suggested. Application of the new method to a set of data 
gives results different from those obtained by the Rowland- 
Hoare methods, from those obtained by the Regnault- 
Pfaundler method (which, despite previous use, should not 
be applied to such data), and from those obtained by the 
method of Dickinson (which is recommended by the 
ASTM standards). The differences are significant in pre- 
cision work. A brief critical review of graphical methods of 
obtaining the temperature correction concludes the paper. 
This review includes a graphical analog of the new method 
mentioned above. 


function,’’ depends upon that between the calo- 
rimeter and the surrounding jacket. The source 
of heat may be a hot solid as in specific heat 
measurements, an ignited bomb as in combustion 
experiments, or simply a hot wire as in certain 
electrical equivalent of heat experiments. The 
calorimeter ‘‘fluid”’ frequently is water, but may 
even be copper. Following the usual practice, 
we shall suppose that the temperatures of the 
source and the fluid are uniform throughout the 
respective bodies. 

In symbols the above assumption is repre- 
sented by the expression 


d0/dt= (1) 


in which © is the temperature of the source of 
heat, @ is the temperature of the calorimeter 
fluid, and 6; is the temperature of the jacket. 
The two functions are expanded into power 
series in terms of their respective arguments, so 
that 

F(O—6) =h+b(O—6)+--- (2) 
and 

(3) 
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Quadratic and higher terms may be added if 
necessary; but with proper experimental pre- 
cautions these terms usually are negligible. 
Values of the “heating constant” 4 and the 
“heating coefficient’’ b obviously depend upon 
the nature of the particular experiment to which 
Eq. (2) is applied. The term —c(@—6;) in Eq. (3) 
is the analytical expression of Newton's law of 
cooling in which c is the so-called “cooling 
coefficient’ of the calorimeter. Heat exchange 
may take place by evaporation, radiation, con- 
duction, and convection. Evaporation is made 
negligible and the other losses are made small by 
proper construction of the calorimeter.* Under 
these circumstances Newton’s law is appropriate 
for the small temperature ranges encountered in 
precision measurements. The discussion below 
will be restricted to cases where c and 6; are 
constant throughout an experiment. The term g 
in Eq. (3) depends, for the most part, upon the 
rate of development of heat due to stirring. 
Since this rate is supposed constant, it is con- 
venient to define an effective jacket temperature 
6, so that 
f(0—0;) = (4) 


In calorimetric experiments, the source may 
supply heat to the calorimeter for a limited 
period of time as in cases of electrical heating or 
else for an unlimited period of time as in measure- 
ments of specific heats of solids. These two 
types of experiments in general require different 
methods for correcting the temperature rise of 
the calorimeter fluid. 

In the first type of experiment, temperature 
readings are made at equal time intervals before 
the addition of heat (foreperiod), during the 
addition of heat (heating period), and afterwards 
(afterperiod). The change in temperature during 
these periods should not display discontinuities 
or marked unevenness. To illustrate these re- 
marks, Fig. 1 is drawn for a case of electrical 
heating. Although numerous methods for making 
the correction in such experiments have been 
suggested, a perusal of Part I of this article will 

*W. P. White, The Modern Calorimeter (The Chemical 
Catalogue Co., New York, 1928). H. C. Dickinson, Sci. 
Pap. Bur. Stand. 11, 189 (1914). Dickinson ascribes one- 
fifth of the heat exchange to radiation and four-fifths to 
conduction and convection (p. 195). Actually convection 


is made small at the expense of increased conduction but 
such that the total correction term is as small as possible. 
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Fic. 1. A typical temperature-time graph of calorimetric 
data with a limited heating period. 


reveal that all of these are modifications of a 
single general scheme. 

The fundamental problem discussed in Part II 
is that in which heat energy is being liberated 
continually to the calorimeter fluid over an 
unlimited heating period. Data for such an 
experiment consist of temperature readings at 
equal time intervals before the unknown quantity 
of heat is supplied to the calorimeter (fore- 
period) and during the time that the source 
supplies heat (heating period). The heating 
period has been considered to end when the 
temperature-time curve begins to obey, within 
experimental error, a pure exponential relation. 
Our analysis of this problem indicates that this 
assumption is not valid. Even though the 
amount of heat transferred from source to fluid 
per unit time continually diminishes, at the same 
time the rate of loss of heat to the surroundings 
diminishes as well; so that the first effect seldom 
becomes negligible with respect to the second. 
It is for this reason that the exponential char- 
acter of the heating period after sufficient length 
of time does not imply pure Newtonian cooling. 
The straight line portion of the @ vs. ¢ curve 
beyond the maximum as illustrated in Figs. 2 
and 3 is not a cooling curve. For precision 
work, therefore, the method of Part I should not 
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be applied to specific heat and bomb calorimetric 
data. .Rowland originally recognized this diff- 
culty and suggested a method based upon some 
knowledge of the heat function F(Q@—86). An 
examination of this and similar methods reveals 
fundamental errors. Consequently, in Part IT we 
develop a method which we believe to be 
basically better. 

Equation (2) offers a classification of these 
correction methods. In Part I experiments are 
considered for which F(©—@) is zero during fore- 
period and afterperiod but is not otherwise 
specified. Part II deals with two cases in which 
F(©—6@) is specified; in the first, 4 is constant 
and & is zero; in the second, h is zero and 6 is 
constant. Finally several graphical methods of 
making the correction are discussed critically in 
Part IIT. 

Because a large number of symbols is necessary 
for the discussion of these correction methods, 
a table has been appended to facilitate the 
reading of the paper. The complete meaning of a 
given symbol must be obtained from the text, 
since the definition in the table is necessarily 
limited. 


TABLE OF SYMBOLS 


te" time during foreperiod, heating period, and 
afterperiod 

to, tm time at beginning, and at end of heating period 

| ae time when calorimeter temperature is a maxi- 
mum 

tint time when calorimeter temperature-time curve 


has a point of inflection 

a’, 8, 0" calorimeter temperature in foreperiod, heating 
period, and afterperiod 

9a, “‘average’’ calorimeter temperature each 


period 

Ao, Om calorimeter temperature at beginning and at 
end of heating period 

6; jacket temperature 

0; effective temperature of surroundings (jacket 


temperature plus a correction term due to 
rate of stirring) 

Sox maximum calorimeter temperature 

Ging temperature at point of inflection of calorimeter 
temperature-time curve 

o* ideal calorimeter temperature that would exist 
if there were no heat exchange with sur- 
roundings 

i) temperature of source of heat 

e* ideal temperature of source in absence of ex- 
change with surroundings 

@ o=0,* temperature of source at beginning of heating 
period 
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v’, v” slope of foreperiod, and afterperiod curve 

a slope of calorimeter temperature-time curve at 
the point of inflection 

a cooling coefficient for the source of heat 

b heating coefficient for calorimeter 

c cooling coefficient of calorimeter 

g rate of increase of calorimeter temperature due 
to stirring 

h heating constant of calorimeter 

p abbreviation for }(a+4h+c) 

€ abbreviation for p—(p?—ac)! 


Part I 


Limited heating period with heating function 
unspecified 


Measurements of latent heats of fusion or 
vaporization or of the electrical equivalent of 
heat are examples of the class of experiments in 
which heat is supplied to the calorimeter during 
a limited time. If the time rate of change of 
calorimeter temperature cannot be formulated 
explicitly for such experiments, then the methods 
discussed in this section are used. 

During the foreperiod of such an experiment, 
1+1 successive temperature readings 6), 41, 

-, are made at the times fo’, «++, ty’. At 
the time ¢,/=f , the quantity of heat to be 
measured is allowed to begin to flow into the 
calorimeter. During the ensuing heating period, 
m—1 readings of temperature 6;, 02, +++, On—1 
are obtained at times f), fo, «++, fm—1. Finally, at 
time t,,=to'’ the addition of heat is stopped and 
during an afterperiod n+1 readings of tempera- 
ture 0)", ---, are taken at times 

-+, f,’". Usually the time intervals between suc- 
cessive readings are kept equal in length through- 
out the whole experiment. 

Clearly the observed rise in temperature of the 
calorimeter is just the difference in temperatures 
at the beginning and end of the heating period, 
namely @,,—4. During this rise, however, heat 
exchanges with the surroundings occur so that a 
correction term C must be added to get the 
change in temperature A@ of the calorimeter due 
to the unknown quantity of heat alone. Sym- 
bolically, 


AO =O», — (5) 


We wish to examine methods of evaluating C 
from the experimental data. Of course, the 
accuracy of the value for A@ depends not only 


559 


- 
j 
| 
| 
) 
| 
| 


upon how well C can be estimated, but also upon 
the precision with which @, and 6) are found 
from the data. 

Since the heating function need not be specified 
explicitly let F(Q—0@)=F(t) so that Eq. (1) 
becomes 


d6/dt = F(t) —c(@—4,). (6) 


None of the heat being measured is allowed to 
enter the calorimeter during foreperiod and after- 
period ; hence during these periods F(t)=0 and, 
upon integrating Eq. (6), 


0’ = 0+ exp[ —c(t’—t,’) ] (7) 
and 
= O.+ (00 —&) exp[ J. 


For the heating period, on the other hand, 
F(t) #0 so that 


0+ f Fndt— (8) 
to to 


For t=t,, the first integral on the right of Eq. (8) 
is seen to be just A@ and the second integral is 
the correction term 


tm 
c= c(0—0,)dt. (9) 
to 


An average temperature @, for the heating 
period is now defined such that 


0.4 = (tm — to) Adt. (10) 
to 


Suppose the times to be reckoned in units 
equal to the time interval used; then ¢,,—f9=m 
and 

C=c(0,—O)m. (11) 


Difficulties in finding c, @., and 6 from the 
observations arise mainly because the tempera- 
ture readings are made at finite time intervals. 

It should be possible to get c and 6, from either 
foreperiod or afterperiod data when analyzed by 
Eqs. (7). To get both ¢ and @ from the data of 
one period requires, however, sufficient accuracy 
to detect variations of d@/dt with temperature 
within this period. Usually these variations are 
so small that the graphs representing foreperiod 
and afterperiod data may well be fitted with 
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straight lines. In such cases both sets of data are 
used to evaluate and 

Evidently the slopes v’ and v”’ of these straight 
lines may be considered to equal the slopes of 
the exponentials at temperatures 6,’ and @,”’ 
within foreperiod and afterperiod, respectively. 


Thus 


y= dt').= —c( Ba (12) 
and 
= ‘dt’’) = —c( 2” 


where @,’ and are appropriate ‘‘average”’ 
temperatures. 

If v’, v’’, and can be obtained from the 
experimental data then the constants c and 4, 
may be computed from Eqs. (12); so that 


C= (v’ —v"’). — (13) 
and 


= (v' —v""0,')/(v' —v"’). (14) 


Finally these values of c and 6 are substituted 
in Eq. (11) to get 


C= — mv’ — (v' —v"’) (m0, — 82’) (15) 
= — mv" —(v' —m8,) — 


Many of the methods proposed in the literature 
agree in Eqs. (15) but differ in the manner of 
computing v’, 01’, Oa, Ao, and On. 

In view of the assumed linearity of the fore- 
period data it is clear that by fitting a straight 
line to these data both v’ and 8,’ may be deter- 
mined. Application of the method of least squares 
yields the result that 


r=l 
(2r—1) 6," (16) 


r=0 


and 
6,/=(1+1)— 6,’, (17) 


r=0 
where 6,’ represents the observed temperature 
at time f¢,’ in the foreperiod. Similar relations 
may be found for v” and @,’’. Equations equiva- 
lent to these are used by Roth.‘ 

A second method that has been used to obtain 
v’ and @,’ may be understood by supposing the 
end points of the intervals within the foreperiod 
to be joined by straight lines. It is now assumed 


4W. A. Roth, Thermochemie (Sammlung-Géschen, 
Berlin and Leipzig, 1932), p. 26. Note that Roth applies 
his method to bomb calorimetric experiments for which the 
heating period is not limited. 
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that v’ is equal to the arithmetic mean of the 
slopes of these lines and that 86,’ is equal to the 
arithmetic mean of the temperatures at their 
midpoints. Evidently this method may not yield 
an adequate value for v’ since the value obtained 
depends entirely on the observed temperatures 
at the ends of the foreperiod. Of course, v’’ and 
6,’ may be determined in a similar manner. 
Schultes and Niibel® propose this method and 
claim that the value of 6,’ so determined is 
superior to that given by Eq. (17). Regnault and 
Pfaundler® find v’ in the above manner and use 
Eq. (17) for 6,’. 

Once suitable values for v’, v’’, and 
have been found, then 6) and 6, may be com- 
puted from the relations 


= + 3lv’, (18) 
Om = — nv (19) 


as suggested by Roth and by Schultes-Niibel. 

Three methods have been proposed for com- 
puting @,. First, the mean temperature in the 
heating period is satisfactory if the data in this 
period obey a linear relationship. Thus 


r=m 


6.=(m+1)"' > (20) 


r=0 


Second, Regnault-Pfaundler, Roth, and Schultes- 
Niibel all assume that @, is the average value of 
the temperatures at the midpoints of the straight 
lines between successive observed points in the 
heating period. This averaging leads to a slightly 
different expression, namely 


r=m—]) 


6.=m"| |. (21) 
r=1 


Finally, White’ suggests that Simpson's rule may 
be used to get @,. For instance, if m is even 


(09 
+++) (22) 


It must be remembered that 6, is defined by 
Eq. (10) in terms of the area under the @ vs. ¢ 
curve for the heating period. For an approxi- 
mately linear curve, the simple expressions in 


5 W. Schultes and R. Niibel, Die Warme 58, 15 (1935); 
Brennstoff-Chemie 15, 466 (1934). 

6 L. Pfaundler, Ann. d. Physik 205, 102 (1866). Griffiths’ 
arithmetical method is equivalent to this one. See Glaze- 
brook’s Dictionary of Applied Physics, Vol. 1, p. 61. 

7W. P. White, J. Am. Chem. Soc. 48, 1146 (1926). 
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Eqs. (20) and (21) are quite adequate; for some 
experiments more laborious ways of estimating 
the area may be needed. 

It is obvious that simplifying approximations 
in the calculation of any one of the seven con- 
stants can often be made for a given set of 
experiments. Although any such simplifications 
can be adopted only for particular problems, 
some suggested approximations are of sufficiently 
general interest to deserve mention here. 

Schultes and Niibel® observe that if v’ and v” 
are taken to correspond to 6) and @,, instead of 
6,’ and @,"’ and if @, is given by Eq. (21) then 

C= — mv’ —(v"’—v') {3 > (0,—00)}. 
= (23) 

Langbein® notes that if the rise in temperature 
in the heating period takes place so rapidly that 
the elevation occurs during the first interval’ and 
if 0.’ and then 


C= — mv’ — }(v’—v”). (24) 


On the other hand, Moser!® observes that Eq. 
(15) may be written 


C= —G(v'—v"’), (25) 


where G is nearly constant for a_ particular 
apparatus and a given class of experiments. 
Actual values for G (from bomb calorimetric ex- 
periments) varied from 1.0 to 1.5. 

White’ points out that if 6,’’— 96, is small it is 
convenient to write Eq. (15) 


C= — mv’ —mk(0,"’ —0,). (26) 


Then, if k may be considered constant for a set 
of experiments, a single preliminary determina- 
tion of k eliminates the necessity of taking fore- 
period data for each run. 

Graphical analogs of some of the foregoing 
analytical methods are discussed in Part III. 


Part II 
Heating function specified 


Sometimes the unknown quantity of heat is 
supplied in such a manner that the time rate of 


8 F. Stohmann, C. Kleber and H. Langbein, J. prakt. 
Chemie 147, 503 (1889). 

® Strictly this is impossible when the heating period is un- 
limited. Langbein applies his method to the class of experi- 
ments discussed in Part IT. 

10H. Moser, Physik. Zeits. 37, 529 (1936). 
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change of calorimeter temperature can be ap- 
proximated by a simple analytical function. This 
fact allows another approach to the problem of 
calculating the temperature correction. 

To illustrate this approach, we shall consider 
the case in which the heating function F(O—@) 
is constant over a limited time (the heating 
period) and zero otherwise. This situation may 
occur, for example, in measuring the electrical 
equivalent of heat. Eq. (1) then becomes 


(27) 
Its solution may be written 
—c(t—t,) ], (28) 


where the point (/,, 8.) lies in the heating period. 
The constants ¢ and @& can be evaluated from 
foreperiod and afterperiod data as described in 
Part I. 

The expansion of the exponential in Eq. (28) 
about a properly chosen value of ct, through the 
linear term is sufficiently accurate for many 
experiments. Expand, insert the end values of 
the heating period, and eliminate ¢, between the 
resulting two equations to get 


h(t» = ty) = Ao +c(0,— )(tm — to). (29) 


The temperature 6* which the calorimeter and 
contents would have had in the absence of 
stirring and of heat exchange with the surround- 
ings is the solution of Eq. (27) with the correction 
term omitted, namely 


6* — =h(t—to), (30) 


where 6*=6)*=6) when (=fy. Finally, the de- 
sired corrected temperature rise (0,,*—69*) is 
obtained from Eqs. (29) and (30). Thus 


— = On, — Oo 94) (tm—to). (31) 


Equation (31) has the same form as the relation 
obtained from Eqs. (5) and (11) in Part I. 
Generally the term / in Eq. (27) predominates 
over the small correction term so that the @ vs. ¢ 
curve is nearly linear over the heating period. 
It is then sufficiently’accurate to let 0,= ) 
so that 


= = 6,,, 69+ ty) (09+ é.. 20, (32) 


562 


Jaeger and von Steinwehr" and Hoare” obtain 
this relation. Of course, better approximations 
may be found. 

In his famous investigation on the mechanical 
equivalent of heat, Rowland develops a method 
in which the nature of the heat exchange between 
a hot body immersed in the calorimeter fluid and 
the fluid is assumed to be known. Hoare” gives 
a more mathematical treatment of this problem. 
The Rowland and the Hoare treatments, how- 
ever, are based on a fallacious assumption and 
lead to erroneous results as will be shown later." 
Because of this fact and because the Regnault- 
Pfaundler treatment is improper for experiments 
on specific heats, bomb calorimetry, and the like, 
we believe an extended discussion of the Rowland 
problem is justified. 

If results good to one-tenth of one percent or 
better are desired, then great care must be 
exercised in determining the corrected tempera- 
ture rise. In such cases not only are the Regnault- 
Pfaundler and Rowland-Hoare treatments inade- 
quate but also Dickinson’s method,'® which is 
recommended by the ASTM standards, becomes 
seriously questionable. Of course, if errors of 
one-half of one percent or worse are acceptable 
then almost any of these methods may be 
applied. 

We make three fundamental assumptions. 
First, at the initial time fy the charge immersed 
in the calorimeter is at a uniform elevated tem- 
perature and the calorimeter fluid is at a uni- 
form lower temperature.'® As time goes on, the 
charge temperature decreases and the calorimeter 
temperature increases, but always in such a way 
that both temperatures are uniform over the 
entire bodies in question. In practice one may 
conceive of a thin conducting shell enclosing the 
charge-fluid boundary and having a heat capacity 
small compared with that of the charge and that 
of the fluid. Across this shell there is a tempera- 
ture gradient proportional to the temperature 

uW. Jaeger and H. von Steinwehr, Ann. d. Physik 64, 
305 (1921), see p. 334. 

2 F, E. Hoare, Phil. Mag. 29, 52 (1940), see p. 57. 

19H. A. Rowland, Physical Papers (The Johns Hopkins 
Press, Baltimore, 1902), p. 400. The problem presented 
here will be termed the Rowland problem. 

4 See page 565. 

1° H. C. Dickinson, reference 3, pp. 230-234. Also see 
p. 567 of this article. 


-16 Of course, the charge temperature may be lower than 
that of the calorimeter. This situation does not often occur. 


JOURNAL OF APPLIED PHYSICS 


4 
| 
| 


— 


J 
8 
8 
20h 
6, 
s 10 20 2s 
TIME —> 


Fic. 2. Plot of the data for an unlimited heating period 
(see text). Note that an afterperiod does not exist for 
such data. In this figure 6, should read 6,. 


difference between charge and calorimeter fluid. 
Dickinson"? and White’ imply that this assump- 
tion is not wholly incompatible with the true 
situation. More likely models are readily visual- 
ized but these require too complicated analyses. 
At any rate, other modifications are introduced 
due to the act of immersing the charge or the 
process of igniting the bomb. 

Second, the charge temperature © is considered 
to change at a time rate which is proportional to 
the temperature difference between charge and 
calorimeter; so that 


d@®/dt= —a(O—8), (33) 


in which a, the proportionality constant, is a 
“cooling coefficient”’ for the charge. Actually the 
rate is undoubtedly a more complicated function 
of (9—8@); the right-hand member of Eq. (33) 
may be taken as the second term in a power series 
expansion in which the first term and terms 
higher than the second are negligible. Note that 
d® dt does not depend explicitly upon 4. 


17H. C. Dickinson, reference 3, p. 203. 
Is \W. P. White, reference 3, p. 88. 
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Third, the time rate of change of the calo- 
rimeter temperature is supposed to be given by 


d0/dt=b(O— 0) —c(@—94,). (34) 


Here we have expanded the function F(@—@) 
into a power series and neglected the first and 
the higher terms. That the first terms in the two 
series contained in Eqs. (33) and (34) are zero 
seems plausible when we remember that after a 
long time both /dt and d6/dt should approach 
zero. Whether the terms higher than the second 
are always negligible may be seriously questioned. 
But in view of the “idealness”’ of the first assump- 
tion, we may accept the second and_ third 
assumptions as adequate. 

Evidently, if no extraneous heat exchanges 
occur, then from the above assumptions with the 
correction term omitted, the temperature of the 
calorimeter will rise exponentially as shown by 
the upper curve (the 6* curve) in Fig. 2. On the 
other hand, if exchanges do occur as in the case 
of observed data, then (for @)<6;<@max) the 
temperature @ rises at first more rapidly than 6* 
and later more slowly than 6*; it reaches a 
maximum value @y.x and finally falls off slowly 
toward @,. This is shown by the lower curve (the 
6 curve) in Fig. 2. If such observed data be fitted 
by the analytical expression of a @ vs. ¢ curve, 
constants of the corresponding 6@* curve are 
thereby evaluated so that the desired corrected 
temperature rise (6,* — @9*) may be calculated. 

Equations (33) and (34) yield the standard 
differential equation 


(35) 


where 2p=a+6+c and g=ac. The solution may 
be written 


where «= p—(p*—q)!, p?>q,'* and A and B are 
constants of integration. 

Clearly the straight line portion of the 6-curve 
beyond the maximum is not a cooling curve. Since 
2p>e, after a sufficient length of time, the third 
term in Eq. (36) is negligible with respect to 
the other terms and the @ vs. t curve then depends 
upon e-*' rather than upon e-*. Figure 3 shows 


1° Since a, b, and ¢ are positive, solutions of Eq. (35) in 
which p?=q are inadmissible. 
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how much the slope in this region may differ 
from the slope of a cooling curve. From Eqs. 
(33) and (34) it follows that 


= (dé, curve 


=(2p/a)(d0/dt)o curve=(atb+c)a/a (37) 


at the point of inflection which is usually in the 
straight line portion of the actual curve. Evi- 
dently the two slopes are practically equal only 
if (6+c¢)<a which is rarely satisfied.2° Hence it 
is improper to use data after the maximum tem- 
perature has been reached as afterperiod data 
from which to obtain c and 6 in the manner 
discussed in Part I. 

To obtain the corrected temperature rise, the 
equation resulting from Eqs. (33) and (34) with 
c=0 is solved. Thus 


= A*+ (38) 


where A* and B* are constants of integration. 
If, at and then Eq. (34) 
with c=0 and Eq. (38) yield” 


— = b( Oo — 80), (a+b). (39) 


In general it is difficult to use Eq. (39) in this 
form. However, it is easy to develop an expression 
for the corrected temperature rise in terms of 
and Using Eqs. (34), (39), 
and the expressions resulting from Eq. (36) 
evaluated at ‘=0 and again at t=fmax: 


— = (09 — &) 


‘+(2p—e)(at+b)y '(Onax— exp(€tmax)- (40) 
Additional useful relations are 
xe 7=ye ‘= — a(tmax — tk) (Omax — 9%) (41) 
and 
lint 2(tmax — tk), (42) 


where — te), y=(2p—e)(tmax a is 
the slope of the @ curve at the time fjn¢ where it 


20 Apparently a and > are inversely proportional to the 
heat capacities of charge and calorimeter, respectively. 
From the well-known equation of mixtures there then 
results 

a/b=(@Qo— ©.,*)/(0..* — 40), 


where ©.,* =@,,* is the final temperature of the mixture. 
Equation (39) follows at once. Evidently a/b usually is 
not very large. 

21 An equation involving @ instead of (@max—) is 
readily obtained. Other forms of possible use in certain 
cases may be suggested by examination of the various 
relations given in the text. 
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Fic. 3. Upper part of the 6 curve in Fig. 2. Note that 
the straight line portion beyond the maximum is not a 
cooling curve. A cooling curve which is appropriate to the 
given data (c=0.006 int~') has been drawn through the 
point of inflection. 


has its point of inflection, and ¢, is the time at 
which 4. 

For many sets of data wherein values of 6 have 
been observed over all three periods, the quanti- 
ties 90, Omax, tmax, @ and v’ may be obtained di- 
rectly. Although readings of @ and ¢ in the 
neighborhood of the point of inflection may be 
fairly precise, the point itself is difficult to 
locate since the curve is so nearly linear in this 
region. Values of 6 over the principal rise in the 
heating period are often not precise and may be 
affected by lags in making the readings as well as 
other factors mentioned above. Hence applica- 
tion of the theory to such data is difficult and 
sometimes ambiguous. By trial and error, values 
of 6 and ¢ are found which satisfy not only 
eq. (41) but also Eq. (36) over a large portion 
of the heating period and the “‘afterperiod.” If an 
unambiguous value of @, can be so determined ; 
then ¢ may be found from foreperiod. data by 
K-q. (12) of Part I, ¢ and 2p from Eqs. (41), and 
finally the corrected temperature rise 
Eq. (40). 


from 
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In another type of experiment” the quantities 
69=6, to, a, and sometimes rough values of 
Omax ANd tmax are found. Readings of sufficient 
number are taken.so that the exponential char- 


acter of the ‘‘afterperiod”’ is clearly evident. On 
the other hand, no data are obtained over the 
principal rise of the heating period. In this case, 
values of (@— 6) over the ‘‘afterperiod” are fitted 
with a pure exponential since the third term on 
the right-hand side of Eq. (36) is negligible in 
this region. Thus e and A are determined. By 
trial and error values of @:ax and tmax are found 
so that both Eqs. (36) and (41) are satisfied and 
thus (2p—e) is obtained. If c is found from 
foreperiod data, then Eq. (41) gives the corrected 
temperature rise. 

A better fit for data of the first type can be 
made if c has been determined for a given 
apparatus from preliminary experiments. Then 
6, can be found from foreperiod data. If readings 
during the heating period are sufficiently accu- 
rate, (, may now be read directly from the 6 
curve. Otherwise a plot of In(@—6,) for the 
“‘afterperiod”’ is used to determine e€; in this case, 
Eq. (41) then gives &. On the other hand, if ¢ 
is small and the ratio of the heat capacities of 
charge and calorimeter can be estimated, then 
Eq. (37) gives the approximate slope of the 
cooling curve in the neighborhood of the inflec- 
tion point.” Both ¢ and 6, can now be computed 
by the method of Part I. ° 

Even if c and 6 are known, it is not, in general, 
correct to use the method of Part I to compute 
the desired 6* curve from the observed @ curve. 
Consideration of Eqs. (33) and (34) shows that 


t 
= -{ (43) 
0 

Since ©* and © are not equal over the entire 
heating period, it is incorrect to neglect the 
second term on the right-hand side of Eq. (43). 

It is convenient to mention at this point the 
approximations in the treatments of Rowland 
and Hoare. The latter deduces an expression for 


2 F. M. Jaeger and E. Rosenbohm, Proc. K. Akad. Amst. 
30.8, 965 (1927). 

23 An examination of Eqs. (37) and (41) shows that an 
error in @ has a greater influence in the first equation. A 
correction of the observed average ‘‘afterperiod’’ slope 
can be found from the relation 
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6,* by neglecting in Eq. (43) the term involving 
©* and © and by replacing @ under the integral 
by 6*.24 Since, for small values of c, the effect of 
the second approximation is not important, 
Hoare’s equation usually gives too high a value 
for 6,*. However, instead of using the equation 
directly (which is difficult in practice), Hoare 
develops a graphical solution similar to that 
originally proposed by Rowland. This graphi- 
cal solution involves additional approximations 
which sometimes override the approximations in 
the derivation of the equation, so that the net 
result is too low a value for 6,*. 

For many experiments where (c—e) is very 
small compared with (a+b) the method based 
upon Eq. (36) may be simplified. The initial 
temperature 4) is supposed to be equal to 4; 
that is, with the stirrer in operation, the calo- 
rimeter and fluid are allowed to come to constant 
temperature before the charge is dropped in or 
before the bomb is ignited. The first term in 
K-q. (40) is then equal to zero. 

Over the “‘afterperiod” the last term in Eq. 
(36) is negligible to the approximation we are 
considering ; hence the temperatures at the inflec- 
tion point (8jn¢) and at fmax(8,) will be given by 


Bint —O0=A exp(— éfins) (44) 
and 


0,—0,=A exp(—étmax), (45) 
respectively. Remembering Eq. (42) there results 
exp max) = (0, — Oo) (Bins — 9). (46) 


Now if (c—«)<(a+b), then Eqs. (40) and (46) 
yield 


(6, 90) (Omax = Ao) = Ao). (47) 


If the “‘afterperiod”’ is linear, eq. (47) becomes 


6, *_ = (Omax 


x (atmax) (Bing — Ao) } . (48) 


A simple graphical construction for this method 
is given in Part IIT. 

It is, of course, difficult to give an estimate of 
the amount of error involved in any proposed 
method of evaluating the corrected temperature 
change, since this error depends upon the con- 

* The first of these approximations is equivalent to 
replacing the assumption expressed in Eq. (33) by 


d@/dt= 


ine 
zt 
} 
ang 
- 
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stants of the particular apparatus used. In our 
search through the literature we failed to dis- 
cover published data of sufficient extent to allow 
an unambiguous determination of the constants 
in Eq. (36) and so to afford a numerical com- 
parison of the result of this equation with the 
results from other proposed methods. Conse- 
quently, we shall illustrate the errors involved by 
a set of artificial data computed by assuming 
reasonable values of the constants in Eq. (36). 

The curves in Figs. 2 and 3 are plots of Eqs. 
(36) and (38) with the chosen values 6)=0°, 
time intervals, «=0.005 int~', 2p=0.505 
int~', and c=0.006 int“. Other quantities in the 
above discussion then have the values: a =0.4167 
int-', b=0.0823 int“, a= —0.004510 deg./int, 
6,=1.7082°, 2.6532°, Oings=2.6193°, tmax 
= 11.3034 int, and ¢;,;= 20.6068 int. For this data 
the corrected temperature rise is 0,.* = 2.6931°. 

The “corrected infinity temperature’ 6; cal- 
culated by the method of Part I (assuming the 
“afterperiod” data to give a true cooling curve) 
should be less than 6,*, since in this calculation 
the effect of the transfer of heat from the charge 
to the calorimeter fluid during the “‘afterperiod”’ 
is neglected. The corresponding temperature 04 
from Hoare’s equation should be too high; on the 
other hand, the value @g-4” computed by the 
graphical method of Rowland and of Hoare may 
well be too low. A similar analysis shows that the 
value 6p calculated by Dickinson’s method may 
also be expected to differ from @,*. The results 
of such errors in the case of our artificial data 
are indicated in Table I. The approximate 
method suggested on page 565 should not be 
used for the above data. Application of this 
method to similar data (for which 6)= 6) gave a 
result only 0.1 percent too high. 

In conclusion, it must be noted again that the 
differences among the results of various correc- 


TABLE I. Restlts of applying various correction methods to 
illustrative data. 


“CORRECTED 


INFINITY 
TEMPERATURES” ©) ERROR 
2.6931 0 
0; 2.6859 0.27 low 
On 2.7000 0.26 high 
Orn 2.6884 0.17 low 


Op 2.6979 0.18 high 


566 


tion methods are significant only in precision 
work where a final error of not more than 0.1 
percent is desired. 


Part III 
Graphical methods 


Graphical methods for determining the cor- 
rected rise in temperature in calorimetric experi- 
ments are designed to avoid undue and tedious 
calculations. To question the relative accuracy 
of the algebraic and graphical methods in terms 
of the size of the graph paper needed appears 
trivial to us. Usually, the total corrected tempera- 
ture rise is determined directly by graphical 
methods so that large errors are likely to creep in. 
Occasionally, only the correction term is found 
by graphical methods; in such cases, the evalua- 
tion of the uncorrected temperature rise remains 
an important problem. 

In the original Rumford®*> method, the value of 
the initial temperature 6 is chosen on the basis of 
preliminary experiments so that = 
This method has been modified in such a way 
that 4) is chosen so as to make the areas A; and 
A, in Fig. 4 equal. The corrected temperature is 
now supposed to be (@max— 40). Of course, it is 
not always true that an observed 6,,. is the final 
temperature of the heating period; then it is 
difficult to apply this method. Furthermore, the 
time required to find an appropriate value of 4 
may be greater than that to calculate a correction 
by methods suitable for use with an arbitrary 
value of 4. 

Regnault and Pfaundier® suggest a semi- 
graphical method*® based upon their analytical 
treatment (see Part I). A straight line is drawn 
through the points (v’, and (v”’, on graph 
paper with slope-temperature coordinates. Values 
of v are read off this curve at temperatures which 
are equal to the arithmetic means for the time 
intervals during the heating period. The sum of 
these values gives the desired correction. Evi- 
dently the criticisms given in Part I apply here. 

Jaeger and von Steinwehr*’? determine the 
constants c and @, from the fore- and afterperiod 
8 The Complete Works of Count Rumford (Estes and 
Lauriant, Boston, 1870-1875), Vol. II, p. 376. 

*6 For a similar construction see P. E. Klopsteg, J. Opt. 
Soc. Am. 13, 589 (1926). 


27 W. Jaeger and H. von Steinwehr, Zeits. f. physik. 
Chemie 53, 153 (1905). 


JOURNAL OF APPLIED PHYSICS 


. 
| 
| 
4 
A 


TEMPERATURE 


TIME —> 


Fic. 4. Illustration of the Rumford and Jaeger-von 
Steinwehr graphical methods. 


data as described in Part I. Areas abc and cefgd 
(Fig. 4) are measured with a planimeter or by 
counting squares. The correction term C is ob- 
tained by multiplying the difference between 
these areas by c. If the rise ac is very rapid, so 
that abc is small, they suggest that area acdgh 
may be measured and subtracted from the 
rectangular area abgh. It is difficult for us to 
see that a more accurate evaluation of the 
required area is gotten in this manner. They also 
show that the position of point f is immaterial 
as long as it is located on the linear portion of the 
afterperiod curve. 

In many methods a time ¢, is assumed to 
exist such that the difference in readings on the 
extrapolated foreperiod and afterperiod curves at 
t, is equal to the desired corrected temperature 
rise. For that matter, once t, has been found C 
may be calculated by means of the relation 


C= —v'(t,—to) (tm— tz). (49) 
We wish to emphasize that a tacit assumption is 
made in all these methods: namely, that the rate 
of heat exchange during the foreperiod and the 
afterperiod depends only upon the cooling con- 
stant c and the temperature difference (@—6,). 


VOLUME 12, JULY, 1941 


These methods, therefore, are not to be applied 
to experiments in which the heating period is 
unlimited. 

Ferry® proposes that ¢, should be taken at 6; 
on the heating period part of the curve; Eucken?® 
suggests that 6, should be used instead. Needless 
to say, in general, neither of these methods is 
correct. Roth*® uses a value for ¢, corresponding 
to the mean temperature 4$(69+8,). Smith*! 
erects a perpendicular at such a time that the 
area between the extrapolated afterperiod curve, 
the latter part of the heating period curve, and 
this line is equal to the area between the extrapo- 
lated foreperiod curve, the earlier part of the 
heating period curve, and this same line. The 
methods of Roth and Smith are nearly equivalent 
to the Regnault-Pfaundler method® and ought 
not be used for experiments employing the 
method of mixtures. 

In Dickinson's method? the time ¢, is obtained 
in the following manner. On the plot of the 
calorimetric data lines are drawn at 69 and 6,, 
parallel to the time axis. Then a perpendicular 
to these lines is erected at such a place that the 
area between the @,, line, the latter part of the 
heating curve, and the perpendicular is equal to 
the area between the 49 line, the earlier part of 
the heating curve, and the perpendicular. That 
this procedure yields a proper value for f, may 
be readily shown. 

For equality of the two areas evidently 


tr 
(6 f (0, —O)dt=0. (50) 
to ty 
But 


(A, f f (0—0,)dt 
to to to 


and (51) 


tm tm tm 
f f (0—0,)dt. 
t t t 


r z 


28 Glazebrook’s Dictionary of A pplied Physics, Vol.1, p. 62. 

29 A. Eucken, ‘‘Energie und Warmeinheit,’’ Zandbuch der 
Experimental physik, Vol. 8, p. 50. 

8° Reference 4, p. 30. 

31E, S. Ferry, A Handbook of Physics Measurements 
(John Wiley and Sons, New York, 1929), Vol. 2, p. 77. 

% Roth’s method is identical with that of Regnault- 
Pfaundler if @4'=4, and Smith's 
method agrees with the Regnault-Pfaundler formulation 
if 6.’ and @,"" are replaced by the temperatures at times 
4(t;+1,) on the extrapolated foreperiod curve and }(t,+tm) 
on the extrapolated afterperiod curve, respectively. 
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Fic. 5. A simple graphical analog of the analytical 
method described in Part II. It is to be used for bomb 
calorimetric and similar experiments. 


Hence the correction term is just 


C= (0@—0,)dt= — —6,)(t,—to) 


to 


—€(Omn — (tm —t,). (52) 


The constants ¢ and 6 are found as in Part I. 

Dickinson attempts to apply his method to 
bomb calorimetric experiments and derives an 
analytical expression for this purpose. Essentially 
he supposes that fp =0 and that 


J f (53) 
0 tr 


Insert Eq. (38) with A* and B* evaluated and 
integrate. There results that t,=(a+)~'. Now 
solve Eq. (38) for @,* to get (0,*—69*) =0.632 
X (6..* — Dickinson replaces 6,.* with 6,, and 
considers the other temperatures to be on the 
6-curve; the value of ¢, corresponding to 6, is 
then used in Eq. (53). For his apparatus Dickin- 
son finds that 0.60 gives better results than does 
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0.632. It is this method which is recommended 
by the ASTM as a standard for finding the 
temperature correction® in bomb calorimetry. 

None of the above methods is strictly appli- 
cable to specific heat or bomb calorimetric 
experiments of any kind. The Rowland-Hoare 
graphical scheme is designed to be used especially 
for specific heat experiments by the method of 
mixtures. The involved construction suggested 
by these men is based upon an analysis which we 
have already shown to be faulty (Part II); and, 
therefore, we propose a somewhat different 
method based upon our analysis of the Rowland 
problem. 

Suppose the initial temperature of the calo- 
rimeter 6) just equals @. On a temperature- 
time graph of the observed data (Fig. 5) draw 
lines parallel to the time axis through points 
G (at 0) and A (at @nx). Drop a perpendicular 
from A to B on the line through G and with 
dividers mark off BB’=GB. Now erect a per- 
pendicular to GB’ at B’ and let it intersect the 
6 vs. t curve in C. Draw a line through C parallel 
to the time axis. Extrapolate the ‘“‘afterperiod”’ 
part of the @ curve back to E on AB. With B’ as 
a center and radius BE strike an are which inter- 
sects the parallel through C in E’ and extend 
B’E’ to D’ on the parallel through A. Now with 
B’ as center and the radius B’D’ strike an arc 
which intersects B’C in D. Then B’D is the 
desired corrected temperature rise. 

To show that this construction is valid con- 
sider the similar right triangles FD’B’ and CE’B’ 
in Fig. 5. Evidently B’D’=B’E’- B’F/B’C. But 
B'E'=BE=(0,—6 9) and B’F=BA = 9) 
by construction. From Eq. (42) it follows 
that Therefore B’D’=(0,— 6) 
X (Omax — 90) 9%); and since B’D=B'D’ it 
follows that B’D =(6,* from Eq. (47). 

We record here our pleasure in receiving 
several suggestions from members of both the 
chemistry and physics staffs. We are particularly 
grateful to Professor O. M. Arnold for his many 
excellent criticisms. In conclusion, we express 
here our appreciation for the kind assistance of 
Miss Harriet Peck, librarian at Rensselaer Poly- 
technic Institute. 


33 See ASTM Standards D271-—33, p. 423 (1936). 
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Kinetics of Degradation and Size Distribution of Long Chain Polymers* 


ROBERT SIMHAT 
Department of Chemistry, Columbia University, New York, New York 
(Received March 7, 1941) 


The rate equations for the depolymerization process of linear chain molecules are formulated 


and a general solution is given. Three special cases are considered: (a) Equal disintegration 
probability for all linkages independently of their position in the chain. The results confirm and 
extend the ones previously found from statistical considerations. (b) Preferred breaking at the 
ends. (c) Equal disintegration probability for all chains independently of their size. Equations 


for the decrease of number and weight average molecular weight with time are given in each 


Case. 


N several publications the problem of the size 

distribution of degradation products, gener- 
ated by the splitting of long chains, has been 
considered.' By a treatment, based on statistical 
considerations, an expression for the distribution 
of molecular weights was found as function of the 
molecular weight of the monomer, of the initial 
chain length or distribution of chain lengths and 
of the average degree of depolymerization, i.e., 
the average number of bonds cut per original 
molecule divided by the total number of bonds in 
the original molecule. 

As pointed out at previous occasions, a com- 
plete treatment of the kinetics of degradation has 
two aspects: first a determination of the distri- 
bution of all possible chain lengths in the reaction 
batch at different stages of the reaction; second, 
an investigation of the change of the properties 
of the system, like average molecular weight, 
with time. Acomplete solution contains implicitly 
the answer to both questions. 

Both problems have been considered by Kuhn 
and Klages, respectively, for some special cases. 
However the results are valid only for random 


* Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 

+ Lalor Foundation Fellow 1940-41. 

1W. Kuhn, Ber. 63, 1503 (1930); Zeits. f. physik. 
Chemie A159, 368 (1932); F. Klages, ibid. A159, 357 
(1932); W. H. Durfee and Z. I. Kertesz, J. Am. Chem. 
Soc. 62, 1196 (1940); H. Mark and R. Simha, Trans. 
Faraday Soc. 36, 611 (1940); E. W. Montroll and R. 
Simha, J. Chem. Phys. 8, 721 (1940). Also I. Sakurada and 
S. Okamura, Zeits. f. physik. Chemie A187, 289 (1940). 
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splitting of infinitely long chains. Montroll and 
Simha have given a complete statistical solution 
for random breaking.? This assumption seems 
reasonable in many cases for large degradation 
products. Nevertheless it appears desirable to 
develop a general theory of the breaking of long 
chains comprising both aspects of the problem. 
In case of non-random splitting a_ statistical 
calculation is difficult. It becomes more promising, 
and of course necessary for a determination of the 
time dependence, to start from the rate equations 
for each degrading species. The solution of these 
equations gives the decrease of the amount of 
each component present and therefore also the 
decrease of average molecular weight with time. 


MATHEMATICAL DEVELOPMENT 


Let us assume that at the time ‘=0 there are 
N,(0) molecules consisting of 7 monomers and 
having j — 1 links, where n =7= 2. These molecules 
are exposed to a depolymerization process, e.g., 
hydrolysis or thermal decomposition, charac- 
terized by a set of rate constants k;, where the 
lower index gives the ordinal number of the 
respective bond and the upper one the degree of 
polymerization of the corresponding molecule. 
Then we have the following set of differential 
equations: 

2 The same formalism has been applied to a treatment of 
the removal of substituents from mixed vinyl polymers; 
compare F. T. Wall, J. Am. Chem. Soc. 62, 803 (1940); R. 
Simha, tbid. 63, 1479 (1941). 
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Fic. 1. Fractions of monomers involved in j-mers for 
various values of the reduced time, \=&t, and for initial 


degree of polymerization, n=350, according to Eqs. 


(2a’, 1) and (3’). 


n—l ) 
ins 
dt i=l 
dN,, 1 (n) (n) m—2 (n—}) 
dt i=1 
“dN, (n) eho: 
= (ky j +k; )Nat 
dt 


n—1l) (n- 


j-1 


ki Nj, 


(1) 


it = (k,, +kr-3 )N,, 1 


+h Ne) 


These equations express the fact that a particular 
chain can be generated by disintegration of any 
of the longer ones and is in turn annihilated by 
decomposition. They are somewhat similar to 
those describing the process of radioactive dis- 
integration treated by Bateman, Ruark, Devol 
and others.* Recently Delbriick* has considered 
fluctuations in autocatalytic reactions and has 
obtained similar equations. System (1) however 
is more complicated, due to the fact that a par- 
ticular species can be created by all the preceding 
*H. Bateman, Proc. Camb. Phil. Soc. 15, 423 (1910). 
A. E. Ruark and L. Devol, Phys. Rev. 49, 355 (1936). 
*M. Delbriick, J. Chem. Phys. 8, 120 (1940). The theory 
of polymerization reactions leads to similar equations if the 
process can be described as an addition of monomers; 


compare H. Dostal and H. Mark, Zeits. f. physik. Chemie 
B29, 299 (1935), 
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ones. If both ends of a chain molecule are equal, 


(7) (2) 
then because of reasons of symmetry k; =k,-;. 


Having solved (1), we find the number of 
monomers from the condition that the sum of 
free and bound monomers is constant throughout 
the reaction process. Finally we calculate the 
number and weight average molecular weight 
defined by Kraemer and Lansing’ from the 
equations: 

M, >;jN; Me >; 7*N; 

M. Me JN, 
where Wy denotes the molecular weight of the 
repeating unit. 

We write the solution of (1) in the form; 


N;= > exp (Ait). 
l=1 
Solving the secular determinant for the \, we find 


1-1 
=0, N= for 1=2, 3---n. 


i=l 


The f; are constants determined by the initial 
conditions. The a;, are solutions of the following 
system of linear equations: 


(-> ki 
i=1 


If the above symmetry condition is not fulfilled 
(n) (m) (n) 
we must write k; +, ; instead of 2k; . We 
see from these equations that a;,=0 for j >/. The 
remaining a,;, obey the simpler relationships : 
(—k, 2ke ay =0, 


(2) 


1) 


| 
ut2k, | 


l—2 
(—2 


In this way we find by successive solution the 
coefficients a;, for j7</ in terms of the rate 
constants and of the quantities a,. The final 


5 E. O. Kraemer and W. D. Lansing, J. Phys. Chem. 39, 
153 (1935). 
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result for the N;(¢) appears in the form: 


N,(t) =hja;; exp 
i=1 


i 
thy 


| 
(2a) 

n n | 

Ni(t)= iN,(O)—> iN,(t) | 
if we denote by N,(0) the number of i-mers 
initially present. Equations (2) and (2a) repre- 
sent the most general solution of the rate equa- 
tions (1), except for the case when the solutions 
of the secular determinant are degenerate. We 
shall encounter such a degeneracy in the following 

discussion of special cases. 


DISCUSSION 


1. Let us assume that all the linkages are 
broken at random. In this case k;‘) =k, inde- 
pendent of 7 and j and the rate equations reduce 
to: 

dN, 
-=2k (1, 1) 


t i=)+1 


200 350 
M, 
M. 


Fic. 2. Variation of number average molecular weight 
with reduced time, \, for various initial chain lengths, 
according to Eq. (4a). 
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300 4 
200 
100 4 


Fic. 3. Variation of weight average molecular weight 
with reduced time, A, for various initial chain lengths, 
according to Eq. (4b). 


From (2) and (2a) we find the solution 
N,(t) =e7* — 2h 
j 42542 ake} (2a, 1) 


for n=j=2. If the material was at t=0 ideally 
homogenous, then N,(0)=0 for j7+n. For the 
integration constants we find 


hja;;=(n—j+1)N,(0). 


Introducing 
(3’) 


we can write the solution as: 


1=j=n-1> (2a’, 1) 


| 


This result is identical with that previously found 
from statistical considerations by Montroll and 
the author. In their work a denoted the average 
degree of depolymerization as defined above. 
Here a@ is expressed as function of the time 
elapsed since the start of the reaction. This time 
dependence has been previously derived by 
Kuhn! from the condition that the rate of de- 
crease of the number of linkages is proportional 
to the number of bonds present. That is, in our 
notation : 


d n n 
—~—¥> (i-1)Ni=k (i-1) Ni. 
dt i 1 i=] 


This equation follows immediately from our rate 


571 


| 
) 
| 
| | 
| | Bite 
| 
\ 
\ 
= 
= 
> 


equations (1a). On the other hand, the degree of 
depolymerization is defined as the ratio: 


i=1 


(i—1)N,(0) 
i= 1 


Ni(t)- N,(0) Ni(t)—N,(0) 
i 
- (3) 
(i—1)N,(0) (n—1)N,(0) 
i=1 

where the second identity holds for an initially 
homogenous substance. Inserting from (2a, 1) 
and considering the boundary conditions we find 
after a simple calculation the above result for a. 
This follows also from the rate equation for the 
breaking of linkages. We have in this way 
rederived the statistical results, and furthermore 
found the time dependence of the inhomogeneity 
of the degrading system.® Figures 1, 2, and 3, 
give the weight fraction distribution, the number 
average molecular weight 1, and the weight 
average molecular weight ./, as function of the 
reduced time \=&t. The problem of the size 
distribution produced by the uniform splitting of 
long chains appears in this way completely solved. 

The experiments of Wolfrom and others’ on the 
hydrolysis of potato starch and of methylated 
cellulose and those of Flory*® on the degradation 
of polyesters by alcohols indicate such a mecha- 
nism to be present. However no fractionation of 
the reaction products and subsequent analysis of 
the distribution curve was undertaken. 

Inserting the expression (3’) into the equation 
previously found by Montroll and the author! for 
the number average molecular weight 7 ,,namely, 


n 
 1+a(n—1)' 


we find the equation previously derived by 
Wolfrom, Sowden and Lassettre’ 


(= ) M, n-1 
= / 

Mo Mo n 

® The statistical result involves the trivial assumption 
that the N; are very large numbers. This condition is here 
introduced implicitly by describing the rate process by 
means of differential equations rather than difference 
equations. 

™M. L. Wolfrom, E. N. Lassettre, D. R. Myers and 
J. C. Sowden, J. Am. Chem. Soc. 61, 1072, 2172 (1939). 

*P. J. Flory, J. Am. Chem. Soc. 62, 2255 (1940). 


M, (4a) 
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The decrease of the weight average molecular 
weight with time is found from (3’) and the 
following equation derived by Montroll and 
Simha :! 


My na?+2(1 —a)[(1—a)"+na—1] 
My 


2. The other extreme case is that of preferred 
breaking at the ends, i.e., the probability for 
breaking end bonds is so much larger than the 
corresponding rates for the interior of the chains 
that one has during the whole cracking process 
only very large components and monomeric 
material in the mixture, without any appreciable 
amount of intermediate sizes. There are indica- 
tions according to experiments of Mark and 
Raff® that the thermal depolymerization of 
polystyrene in liquid phase at about 350° 
proceeds in this way. As a first approximation 


(i) 
we assume that ky; =k;1.=k , while all the 


other rate constants are zero. The differential 
equations then reduce to 
dN; 
—— = N;; 
dt 


(4b) 


na" 


(1,2) 
They are of the general type treated by Bateman* 
and have the solution 
N;=>D exp [—2kt] 
i=j 
with 
c;= N,(0) 
x— —— —, (2a, 2) 


This result holds also for 7=2 if we denote the 
rate constant for the dimer by 2k™. A separate 
treatment of the case of equal disintegration 


probabilities, k =k” =k =k gives the solution : 
on] 
N;=N,(0)e-7-— 
V.=N,(O)e —+N,(0)R 
n (2a, 2’) 
n- 


| +2N,(0)R> 


} 


| 
| 


+(n=2)(e- 


r=n—?2 


*H. Mark and R. Raff, private communication. 
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with z= and 


z\* 1 
R= —2,Qn-2 
£6) i! 


Equations (2a, 2) and (2a, 2’) describe the 
course of the degradation process for an initially 
homogenous material of degree of polymerization 


1 

n | (n—2)! 

M. 


gr? 3” 
Mo n 


n. For large values of n and for the first stages of 
the reaction (small values of s) the number of 
monomers produced, increases linearly with time: 


Ni = N,(0)z. 


For the number and weight average molecular 
weight we find, respectively : 


(4a’) 


2 


(n ! 


2 o | 2R 
+{1- ) +—. (4b’) 
n n—3 n 


Figure 4 shows a graph of these functions for 
n= 300. As to be expected, the weight average 
decreases very slowly because the monomers pro- 


("yn) 


“| \ 
| 

| 


Mean Degree of Polymerization ~ (“yn or 
8 


2 
Fic. 4. Variation of number and weight average molecu- 


lar weight with reduced time, z= 2&t, for n = 300, according 


to Eqs. (4a’) and (4b’). 
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duced*contribute very little to M/,, while they 
influence relatively strongly the decrease of .M/,. 

3. In our first case the disintegration proba- 
bility for a chain was proportional to its size or 
more correctly to the number of linkages by 
which it is formed. Let us now consider a case 
where the chance of decomposing is roughly equal 
for all molecules independently of their size. The 
probability for a particular bond in a particular 
chain being cut is then inversely proportional to 
the number of linkages in that chain. To the 
author’s knowledge no experimental data are 
available as yet, but such an assumption might 
give an approximate description of the catalyzed 
cracking of long chain molecules like poly- 
styrene.'” It seems possible that surface catalysis 
affects each molecule approximately in the same 
way, at least within certain size limits, and thus 
justifies the above assumption. Our rate equa- 
tions then read as follows: 


dt i=j 4 


The first step, the formation of a molecule 7 by 
disintegration of an i-mer has a rate constant 
2k/i—1 because there are two possibilities for 
this step. k is the rate constant for the decompo- 
sition of chains of any length. The previous solu- 
tion (2) and (2a) does not apply here because of 


10T am indebted to Professor H. Mark with whom I 
had the opportunity to discuss this point. 
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Fic. 5. Variation of fractional number of chains, 
N,/N,(0), with reduced time, r=2kt, according to Eq. 
(2a, 3). 


the degeneracy of the determinant. Subtracting 
the respective equations for N; and N;,,; from 
each other we get 


dN; dNiat 2 
+kN;= = 
dt dt j 


Therefrom 


C; is an integration constant. If again for f=0, 
N,=N,(0) and all the other N; are zero, C; 
; vanishes. By successive integration we find with 
t= 2kt 
N,(0) | 
1+ 7S, +—S, 07) 
n—1 | 2! 
1 G+) | 


(n—j—1)! | 


(2a, 3) 


G+) | 


‘| 


The coefficients S;“*+" are defined as follows: 


= 1, 


n—2 mi-l 1 


vi=jtk—1 Vy Vo Vi. 
From the relationship 
2k 
Fj41—-F;= ——— ‘Nj, 
j-1 


we obtain a recurrence formula for the S,@*+ 


The values for the upper and lower indices are 
limited by the condition k+j+1=n. It is pos- 
sible to express these coefficients in terms of the 
¥ function. For integer values of the argument it 
is defined by 


1 1 


n 


v(n+1)=y(n)+ 


n+1 


C is the Euler-Mascheroni constant and has the 
value 0.57722. For sufficiently large values of the 
argument we have the expansion :"" 


1141141 141 1 
¥(n) =In n4+—— + +0(-,). 
2n 12m? 120mn* 252 n° 


Fic. 6. Fraction of monomers involved in j-mers for 
various values of the reduced time, rt, and for n=300, 
according to Eq. (2a, 3). 

"E,. Madelung, Die Mathematishen Hilfsmittel des 
Physikers (|. Springer, Berlin, 1936). 
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We shall make use of this series in the numerical 
calculation of the average molecular weight. 
Now obviously 


From the equations derived in the Appendix we 
conclude 
y(n —2) 


1 
= —2) 4—— 
¥(n —2)¥(n—3) 


j-1 1 n—3 1 
+ 


(5a) 


The use of the infinite series 


1 


6 


facilitates the evaluation of the two sums in (5a) 
or at least of the second one. In the range of n 
values in which we are interested here, the value 
of the second sum is practically equal to that of 
the series. For intermediate values of 7 correction 
terms, starting from the series, can be calculated 
for any desired degree of accuracy. It is possible 
to compute higher terms in the expression (2a, 3) 
in a way analogous to that used in the Appendix 
for S,°*. The computations however soon be- 
come very cumbersome. Therefore, at present, a 
more accurate evaluation of the higher terms has 
been omitted. They decrease quite rapidly with 
their ordinal number j, as can also be seen from 
the following numerical example. Therefore even 
for values of r of the order of magnitude of one 
the accuracy obtainable in this way for N; is 
sufficient for intermediate sizes, e.g., 7=100, if 
n= 300. A rough estimation of the coefficients can 
be made by using the first term of the expansion 


for y(n). If 


k+1 
(j—1)[In (n—k—2)—In (j-1) ] 
and In (n—k—2) =In (n—k—1), then 


1 
(n—k—1)—In (j—1) }*. 
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Degree of Polymerization ~ (MyM,) or (My/m,) 


Fic. 7. Variation of number and weight average molecu- 
lar weight with reduced time, 7, for n= 300, according to 
equations derived in the Appendix. 


This equation can be used when k is small; 7 must 
be large enough to allow the substitution of the 
logarithm for the function. For 7=11 the corre- 
sponding error is about 2} percent. 7 —1, however, 
must be small compared with n—k—2. 

For the computation of N,, M, and M, we 
need the sums 


} N;, IN;j, PN;j. 
j=2 
We get 
n N,(0) 
N;=N,(O)e" eT 2 
| ( D4 rktl 
x TiN — ) S, G+) 
| 2! (k4+1)! 
n—k 1 (j+1) T" 2 (a) | 
Che 
(n—2)! 


and similar results for the other sums. Figure 5 
gives the decrease of the number of various 
species with reduced time. Figure 6 shows the 
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total weight distribution during the first stages 
of the reaction. The final graph, Fig. 7, reproduces 
the decrease of the mean degree of polymerization 
with time. As is to be expected, the curves are 
much steeper for random splitting of bonds than 
for random splitting of chains. 

The above equations have all been formulated 
explicitly for an initially homogenous material. If 
we have a whole distribution then the final result 
for N,(t) is a superposition of partial solutions of 
the type written down, each containing a factor 
N,(0) with n=k=j. If we represent the initial 
conditions by a smooth distribution curve we can 
write an integral between the limits m and 7 
instead. 

Theoretical attempts to obtain the disinte- 
gration constant of a chain molecule as function 
of the position of the bond in the chain, its 
strength and the size of the chain hitherto have 
not been very successful. Pelzer'* has examined 
the way in which an isolated chain molecule can 
distribute excess energy among its linkages. If all 
the masses and bonds are identical, equal proba- 
bilities for end and middle bonds and no depend- 


ence upon the length of the chain is found. 
According to Schmidt,'* a different mechanism is 
operative in the cracking, e.g., of hydrocarbons. 
The course of the reaction in the hydrolysis of 
cellulose can be best described, according to 
Freudenberg and Blomquist,'* by the assump- 
tion that end bonds react twice as fast as those in 
the interior of the molecule. 

It is very desirable to have more experimental 
data on depolymerization reactions. The change 
of the mean molecular weight with time as well as 
the whole distribution curve at different stages of 
the reaction should be determined. From the 
point of view of the theory it appears most 
suitable to study thermal cracking processes, 
because in this way the conditions imposed by 
the theory could be best verified. Such investiga- 
tions might also represent a first step to a study 
of equilibria between polymerization and de- 
polymerization processes." 

The author wishes to express his thanks to the 
Lalor Foundation whose grant of a Fellowship 
made this work possible. 


APPENDIX 


(a) Calculation of 


n—2 
It can be shown that 
k 
1 Pr 1 


Proof: 
1 


1 1 


5) 


On the other hand, using the recurrence formula for pad 


v(k)_ 


2 


2H. Pelzer, Zeits. f. Elektrochemie 39, 608 (1933). 
“ O. Schmidt, Zeits. f. physik. Chemie ‘A159, 337 (1932). 


¥(k) 1 ¥(n+1) 


“4K. Freude nberg, and G. Blomquist, Ber. 68, 2070 (1935). See also Alf af Ekenstam. ibid. 69, 549, 553 (1936). 
% The resulting size distribution has been derived in terms of an equilibrium constant under simplifying assumptions 


by G. V. Schulz, Zeits. f. physik. Chemie B43, 25 (1939). 
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Making use of these expressions and taking the proper limits we find Eq. (Sa) given in the paper. 


(b) Calculation of ~ 1, 2. 


n—2 n 


n—3 ; 1 1 n—3 
¥(j) = — +: = 3C (n— = = C+ (n-3); 
n—3 
The calculation of the sum 
n—3 
373 
requires an evaluation of the following partial sums: 


as can be seen from the results derived under (a), and 


n—3 j—1 


~ ¥(n—3)+y(1) 


obtained by direct computation. The final expression is: 


n—3 1 n—3 


j=2 


vin—3)—— 
1 


In the numerical examples the sum 

n—3 

has been evaluated with the aid of Euler’s summation formula," valid if the function and its derivatives appearing, are 
continuous: 


Bu 


(2k-1) (m (2k-) (() m—— f (2042) (9m 
with m, n=1, 2---;0<8<1. The By are Bernoulli’s numbers, Bo=1, B,=1/6, Bg= —1/30, Bp =1/42. In our applica- 
tions it is sufficient to take two terms in the last sum of Euler’s formula. In the calculation of this and the following 
integrals the expansion for the y function up to terms 1/f* was used. The numerical values of ¥(¢) and its derivatives are 
given in Davies, Tables of Functions. We obtain 


and 
We consider the following sums: 
k~1 
Voy 3 
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and the corresponding integrals 


fy2(t)dt. 


ty-(t)dt and 


We finally have | 
n—2 


For the sake of brevity, the sums over S,’*) have not been written down explicitly, but they can easily be found with 
the aid of the previous expressions. 
Taking in our example n= 300 we find for instance 


n 


and similar equations for the other sums showing their relatively good convergence. 


Consolidation Settlement of a Soil with an Impervious Top Surface 


M. A. Biot AND F, M. CLINGAN 
California Institute of Technology, Pasadena, California 


(Received April 8, 1941) 


An investigation is made of the consolidation settlement for a soil whose top surface is 
impervious and loaded uniformly on an infinite strip of constant width. The settlement under 
the load is accompanied by considerable swelling of the unloaded area on both sides of the load. 


INTRODUCTION 


N a previous publication! the settlement of a 

consolidating soil under an area of constant 
width was obtained when the water is allowed to 
flow freely through the top surface both outside 
and under the load. We are going to investigate 
the same problem when the top surface is com- 
pletely impervious both outside and under the 
load. It is clear that in this case the settlement is 
due to the fact that the water contained in the 
soil flows from under the load to unloaded 
regions. Since it is assumed that the water 
cannot escape through the top surface, this will 
produce a swelling of the unloaded area in the 
vicinity of the load. For the basic theory the 
reader is referred to a recent paper.* 


I. SETTLEMENT UNDER A LOAD WITH 
SINUSOIDAL DISTRIBUTION 


In the previous paper? have been established 
the following equations for the two-dimensional 
problem of the consolidation of completely satu- 

'M. A. Biot, ‘Consolidation settlement under a rec- 
tangular load distribution,” J. App. Phys. 12, 426 (1941). 


2M. A. Biot, ‘‘General theory of three-dimensional con- 
solidation,” J. App. Phys. 12, 155 (1941). 
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rated clay with zero Poisson ratio 


de Oa 
GV*u+G— ——=0, 
Ox Ox 
(1.1) 
de Oo 
GV*w+G——— =0, 
02 


In these equations u and w are the components 
of the displacement of the soil, e= du /dx+ dw ‘dz; 
G is the shear modulus for the completely con- 
solidated clay ; c= k/a coefficient of consolidation ; 
k is the coefficient of permeability; a=1/(2G) 
the final compressibility ; ¢ water pressure incre- 
ment in the pores; p=0/dt symbolic operator. 
These equations must be solved with the bound- 
ary conditions (1) that all variables vanish at 
infinite depth s=«; (2) d0/dz=0 at 
(3) —0o+2Gdw/dz= —A sin Ax at z=0; (4) Ou dz 
+dw dx=0 at s=0. The second condition ex- 
presses that no water flows from the top surface. 

The last two conditions express that at the 
surface the normal stress is equal to the load and 
that the shearing stress is zero. 
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It may be verified that a solution of Eqs. (1.1) 
satisfying the boundary condition (1) is 


u=[( + Cod exp [—(A2+ pc) 
— C3(1—Az)e—** | cos Ax, 
p/c)} 
Xexp [—(°+p/c)'s] 
— C3\ze~** sin Ax, 
=2G(C2(p/c) exp [— (A? +p/c)'z] 
— ] sin Ax, 


(1.2) 


in which the arbitrary constants C;, C2. and C3 


are to be determined so that the remaining 
boundary conditions (2), (3) and (4) are satisfied. 
We find 


(p/c)C2(d?+p/c)!— =0, 
2G7(Ci+ C2) = —A, 
ACi + (42+ p/c)'C2— C3 =0. 


(1.3) 


We are interested in the vertical deflection w at . 
the surface, :=0. This value is from (1.2) and a 
(1.3) 


sin Ax 


(1.4) 


Solving (1.3) for Cs we obtain 


A 
G 


The soil deflection at the surface as a function 
of the time for a sudden sinusoidal load distribu- 
tion is written in operational form 


w= . (1.6) 
2Gd 


The function represented symbolically by this 
equation may be calculated from the following 
operational expressions. 


p?/c*4 exp [A*ct/ 2] 5} 5} 
-1(t)= sinh ( niet) +5 cosh ( | (1.7) 
p?/c*h4 — p/cr?—1 5} 2 2 
1 exp [A*ct/ 2] 5} 5} 
1(t)= —sinh ( +5) cosh ( (1.8) 
p?/c?\4— 5} 2 2 
(1+ p/cr*)} exp ] 
—-1(t)= - PLA(ct)*] 
p/crh?-1 5} 
exp 
(1.9) 
5 
The first two operators are tabulated* while A sin \x ’ 
the third is derived by separation into partial “> Gr 1+ P[A(ct)?] 


fractions and applying the shifting formula‘ to 
the well-known operator 


1 
p'1(t) =——. (1.10) 


(xt)? 


With these results the deflection at the soil 
surface is written 

3E. Berg, Heaviside’s Operational Calculus (McGraw- 
Hill Book Company, Inc., New York, 1929), p. 165. 


*V. Bush, Operational Circuit Analysis (John Wiley and 
Sons, New York, 1929), pp. 130, 191. 
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exp [(1+5')d*ct/2 

52 
exp 

5 

(1.11) 
In these formulas we have 


zr 
P(x) = exp (— &)dé 


| 
= 
f 
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Fic. 1. Settlement of soil surface at various instants of 
time for a uniform load po extending from x=0 to 
Xm 


a tabulated function called the probability in- 
tegral. 

The settlement is composed of two parts: a 
purely elastic deflection which occurs at the 
instant of application of the load and a settle- 
ment due to consolidation which occurs gradually 
thereafter. Since we are interested only in the 
latter we proceed as follows: 

For ¢=0 the deflection is purely elastic, its 
value is 

w,=A sin Ax, (1.12) 
The settlement due to consolidation is 
W,=W— Wj. (1.13) 


For convenience we will write Eq. (1.13) in the 
form 


sin Ax. (1.14) 


Il. SETTLEMENT UNDER A UNIFORM Loap 
WITH DISCONTINUITY 


For practical purposes we are interested in the 
differential settlement occurring when the load 
distribution has a discontinuity. Such a load 
distribution will be represented by the following 
function of x 


Po sin Ax 
g(x) = ———dn (2.1) 
0 


in which case ¢(x)=—po/2 for x <0, po/2 for 
x >0 with a discontinuity at x=0. 

The settlement due to this load is obtained 
from (1.14) by the principle of superposition ; 
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Fic. 2. Settlement of soil surface at various instants of 
time for a uniform load po extending from x=—//2 to 
x=+/1/2. 


2 — 


Fic. 3. Rates of the settlement at the center of the 
loaded area in Fig. 2 when the top surface is pervious to 
the settlement, when the top surface is impervious. 


ol LG 


we find 


w,=— —— sin AxdX, (2.2) 


apo f FU X(ct) 
0 


which can be written in the form 


w, 1 * F(y) 
-—f sin (2.3) 
2apo\et 


with &=x/(ct)!. 
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This integral can be evaluated numerically by 
the use of an approximate analvtical expression 
for the function F(y). It may be verified that 
this function is represented within 2 percent by 


F(y) =1—e-7— ye~*7 — 1.06y*e-*7 


(2.4) 


If (2.3) is integrated with the substitution 
(2.4) we find 


w, 1 1 1 
(~) =—£ log tan-'& 
2apo\ct 
1 


0.28€ 


(2.5) 
By adding a load equal to po/2 and extending 
from x=—*x to x=+2 we obtain a total 
uniform load extending from x=0 to x= as 
shown in Fig. 1. In doing so we do not change 
the settlement because the water is imprisoned 
in the soil by the impervious layer at the surface 
and therefore a load extending from x= — ~ to 
x=-+ cannot produce any settlement. Hence 
the settlement for the load po extending over 
the positive x axis is given by the same expression 
(2.5) which we write 


W, = 


where f() denotes the right-hand side of Eq. 
(2.5). It is clear that the volume of the water in 
the soil being constant, the total area under the 
curve f(£) is zero. 

In order to represent the settlement as a 
function of time it is convenient to introduce for 
the abscissa a characteristic length / which can 
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be chosen arbitrarily. We then have 


2 (ct)itx 1 
Ul (ct) 


Then settlement curves are plotted as function 


of x in Fig. 1 at time intervals corresponding to 
(ct)!/l=1/8, 2/8, 3/8, 4/8, 5/8. 


(2.6) 


Ill. SETTLEMENT UNDER A RECTANGULAR 
Loap DISTRIBUTION 


By superposition we may easily derive from 
the previous solution the settlement due to a 
constant load po extending from x=-—/ 2 to 
x=1/2. Using (2.6) we get 


2 x+1/2 (—) 


The settlement curves are represented in Fig. 2 
as a function of x at time intervals corresponding 
to (ct)!/1=1/8, 2/8, 3/8, 4/8, 5/8. 

It will be noted that immediately after loading 
the settlement at the edge is greatly affected by 
the unloaded section adjacent to it. This un- 
loaded section provides an escape for the water 
which is not as readily available for the water 
under the center of the load. Also since this load 
and deflection are symmetrical to the line x =0 
no water will flow across this line and therefore 
the area under the curve from 0 to * must be 
equal to zero. 

If the ratio of the deflection at x=0 is com- 
pared with the case in a previous paper? in 
which the loading is the same but the water was 
allowed to flow freely from the top surface the 
effect of an impervious top surface will be seen. 
This comparison is made in Fig. 3. 


i 
‘ 
| 
q 
581 


Innovations in Instruments 


New Temperature Indicator 


Increased precision in measuring thermocouple tempera- 
tures, in checking thermocouples, and in measuring small 
e.m.f.’s, may be obtained without increase in instrument 
cost by the use of the new portable temperature indicator 
recently put on the market by the Leeds and Northrup 
Company, 4934 Stenton Avenue, Philadelphia. In the new 
indicator 8667, only a small portion of the range is calibrated 
on a continuously adjustable slidewire ; the rest, adjustable 
in fixed steps, is on a dial switch of ten highly accurate 
resistors. Range is from 0 to 111 millivolts with a limit of 


error of +0.1 millivolt. The increased accuracy without 
increase in cost is obtained by sacrificing the convenience 
of a reference-junction compensator; the only accessories 
required for this self-contained potentiometer, illustrated 
in the accompanying figure, are the thermocouple and an 
ice bath for the reference-junction. 


Insulation and Insulating Varnish 


The Irvington Varnish and Insulator Company, 24 
Argyle Terrace, Irvington, New Jersey, announced the 
following new products: Irvington Varnished Fiberglas, 
glass cloth impregnated and coated with heat resisting var- 
nishes, and used as insulation in the manufacture of electric 
motors, generators and similar equipment; and Harvel 
612-C insulating varnish, a phenol-formaldehyde synthetic 
resin-base varnish made from Cashew Nut Shell Liquid 
which polymerizes upon baking to an infusible, insoluble 
state. Literature concerning these products may be ob- 
tained by writing the manufacturer. 


Variety in Decade Resistance Boxes 


Shallcross Manufacturing Company, 10 Jackson Avenue, 
Collingdale, Pennsylvania, announces a wide range of 
decade resistance boxes, manufactured in increments from 
0.1-ohm steps up to 1,000,000-ohm steps. Some have as 
many as six decades in a single box. Details are given in 
Bulletin No. 825. 
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New Vacuum Tubes 


RCA Manufacturing Company, Inc., Harrison, New 
Jersey, announces that it is making available the following 
new tubes: 

RCA-6SG7 and RCA-12SG7, r-f amplifier pentedes of 
the metal type particularly recommended for use in high 
frequency receivers. They feature high transconductance, 
very low grid-plate capacitance, and two separate cathode 
terminals. The 6SG7 and 12SG7 are alike except for 
heater rating. The heater of the 6SG7 is designed so that 
it can be operated in series with other 6.3-volt, 0.3-ampere 
types; likewise, the heater of the 12SG7 can be operated 
in series with other 12.6-volt, 0.15-ampere types. 

RCA-45Z3, a miniature, half-wave high vacuum rectifier 
of the heater-cathode type. Designed specifically for 
service in a.c./d.c./battery-operated portable receivers, 
it offers the advantages of small size and low heat dissipa- 
tion. Its heater requires only 0.075 ampere at 45 volts and 
its output rating makes it capable of supplying rectified 
power for both filament and plate circuits in light-weight 
receivers utilizing the miniature type tubes. 

RCA-304, a new, miniature type, power amplifier 
pentode suitable for operation with 90 volts on plate and 
screen. It has a power output of 270 milliwatts with 7 
percent distortion. The filament of the RCA-3Q4 has a 
center tap which permits economical usage of this tube in 
series-filament arrangements with other miniature tubes. 

The 6SS7, a remote cut-off, r-f amplifier pentode of the 
single-ended metal type having a 6.3-volt, 0.15-ampere 
heater. This new tube provides for a further degree of 
flexibility in the design of a.c./d.c. receivers utilizing single- 
ended metal types, where the total heater voltage of a 
complement of 0.15-ampere types heretofore available 
would exceed 117 volts. 

The 12SN7-GT, a single-ended, twin-triode amplifier 
having separate cathode terminals for each triode unit. 
It is recommended for use in resistance-coupled circuits 
as a voltage amplifier or phase inverter. 

RCA-930, new gas photo-tube recommended for use in 
sound reproduction and relay applications. Electrically, 
the 930 is like the type 923 with its high sensitivity and 
large response to red and near infra-red radiation. 

The RCA-931, a new type of photo-tube in which the 
photo-current produced at a light-sensitive cathode is 
multiplied many times by secondary emission occurring 
between nine successive dynodes within the tube. It is 
capable of multiplying feeble currents produced by weak 
illumination as much as 230,000 times. Focusing of the 
electron stream is accomplished electrostatically within 
the tube. The 931 employs the S4 photo-surface which has 
higher sensitivity to blue-rich light than to blue-deficient 
light. 

The 8001, a transmitting beam pentode suitable for 
power amplifier, modulator, and oscillator service. The 
8001 is particularly useful in all band transmitter designs. 
In class C telegraph service, the 8001 will provide a power 
output of approximately 230 watts at frequencies as high 
as 75 megacycles. Its plate dissipation is 75 watts. 
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PHOTOGRAPHIC PLATES 
for Scientific Purposes 


—_— sixty kinds of special plates for scientific work are made in Kodak 
Research Laboratories. Eight basic types of emulsion are used, each 
carefully prepared to provide desirable combinations of speed, contrast, 
granularity, and resolving power. They are sensitized in various ways to 
permit photography in spectral regions ranging from the short-wave-length 
ultraviolet to the infrared at 12,000 A. 

The characteristics of these plates are described in the booklet, Photo- 
graphic Plates for Use in Spectroscopy and Astronomy, 4th Edition; a copy 


of it will be sent, free, upon request. 


EASTMAN KODAK COMPANY 


Research Laboratories Rochester, N. Y. 


A New Text of Unusual Distinction 


FUNDAMENTALS OF PHYSICAL SCIENCE 


By KonrAD BATES KRAUSKOPF 


Assistant Professor of Geology, Stanford University 
5 


653 pages, 6 x 9, illustrated. $3.50 


In this outstanding new text the sciences of physics, chemistry, astronomy, and geology are pre- 
sented for the beginning student as a unified field of knowledge, rather than as disconnected 
summaries of the separate sciences. The ideas and materials in the text have been developed in 
the well-known pioneering course at Stanford University. The author’s unusually broad training 
in science and his remarkable literary style have enabled him to produce a well balanced and 
scholarly text that will capture and hold the interest of the beginning student. 


Send for a copy on approval 


McGRAW-HILL BOOK COMPANY, Ine. 
330 West 42nd Street New York, N. Y. 
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New Instrument Booklets 


Recent booklets of interest particularly to the physicist 
include: 


The General Radio Experimenter, May, 1941, issue, con- 
taining articles on ‘‘Radio-Frequency Characteristics of 
the Type 726-A Vacuum-Tube Voltmeter,” and “The 
General Radio Standardizing Laboratory.” Published by 
General Radio Company, 30 State Street, Cambridge, 
Massachusetts. 


Fisher Scientific Company’s The Laboratory, Vol. 12, 
No. 4, featuring a description and step-by-step photo- 
graphs of the procedure for making burettes. Fisher 
Scientific Company is located at 711-723 Forbes Street, 
Pittsburgh, Pennsylvania. 


Science and Appliance, March, 1941, timely items in 
science gathered by The Ohio State University Research 
Foundation, Columbus, Ohio. 


The spring issue of The Educational Focus, published by 
Bausch & Lomb Optical Company. The cover of this 
issue is attractively printed in color. 


“Varnished Fiberglas,"’ a folder with samples of black 
and yellow varnished fiberglas, insulating material manu- 
factured by Irvington Varnish and Insulator Company, 
40 Argyle Terrace, Irvington, New Jersey. 


The Laboratory, Vol. 12, No. 3, featuring a story on 
“The Invention of the Mercury Barometer,” distributed 
by Fisher Scientific Company, 711-723 Forbes Street, 
Pittsburgh, Pennsylvania. 


“University Instruments,” planned especially for use in 
the training of engineering students, are described in the 
leaflet, Electrical Measurements, Vol. 6, No. 5, publication 
of the Sensitive Research Instrument Corporation, 4545 


Bronx Boulevard, New York, New York. 


Three catalogs recently received from the Roller-Smith 
Company, 1766 West Market Street, Bethlehem, Pennsyl- 
vania, bear the following numbers and titles: 3650, 
Outdoor Oil Circuit Breakers; 3940, Subway Oil Circuit 
Breakers; and 4540, Precision Balances. 


Brought together for the first time in a single catalog so 
that it might be more generally available to laboratories, 
technical schools, universities, and manufacturers, ‘““RCA 
Measuring Equipment” is described in a 22-page publi- 
cation recently released by the RCA Manufacturing Com- 
pany, Inc., Camden, New Jersey. An RCA leaflet has also 
been recently received discussing “RCA Dynamic Demon- 
strator, ‘The Working Schematic Circuit Diagram.’”’ 


A news-sheet designed to ‘‘keep you up to date on relays 
and relaying methods” is published by the General 
Electric Company, Schenectady, New York. It is entitled 
“Relaying News.” 


Colloidal Graphite as a High Temperature Lubricant and 
Colloidal Graphite—Focusing Anode Material for Cathode 
Ray Tubes, technical bulletins issued by Acheson Colloids 
Corporation, Port Huron, Michigan. 


Du Mont Oscillographer, February-March, 1941, organ 
of the Allen B. Du Mont Laboratories, Inc., 2 Main 
Avenue, Passaic, New Jersey. Contents, “‘A Single-Sweep 
Stopwatch and Its Use in a Biological Problem.” 


“The Successor to Colorimeters. Coleman Universal 
Spectrophotometer,”’ Bulletin 204 of the Burrell Technical 
Supply Company, 1936 Fifth Avenue, Pittsburgh, Penn- 
sylvania. 


Two new technical bulletins of the Acheson Colloids 
Corporation, Port Huron, Michigan, entitled: ‘Colloidal 
Graphite Lubricants,” and “Colloidal Graphite Dis- 
persions.”’ 


Wheelco Condensed Catalog, a convenient listing of the 
principal items of equipment manufactured by the Wheelco 
Instruments Company, 1929-33 South Halsted Street, 
Chicago, Illinois. 


New laboratory apparatus catalog, XV 41, 174 pages 
describing instruments and apparatus for science labora- 
tories, available upon application on business or profes- 
sional stationery, from the American Instrument Company 
of Silver Spring, Maryland. 


“Technical Books on all Subjects, 1941,’’ catalog of the 
Chemical Publishing Company, Inc., 234 King Street, 
Brooklyn, New York. Subjects range from ‘‘aeronautics” 
to “foods.” 


“Vitreosil Pipes and Fittings,’ Bulletin No. 7, published 
by The Thermal Syndicate, Ltd., 12 East 46th Street, 
New York, New York. The vitreous silica pipes and 
fittings described are “invaluable for extreme chemical, 
thermal and electrical conditions.” 


How to apply radiant heat lamps to industrial drying, 
baking, and heating tasks is told in a new 12-page illus- 
trated booklet A-3817 published by the Westinghouse 
Lamp Division, Bloomfield, New Jersey. The fundamentals 
of radiant heat are presented in the first section of this 
publication, followed by four pages of photographs showing 
typical drying lamp installations. 


The Tuckerman Optical Strain Gauge for measuring 
tension and compression strains as small as 0.00002 inch 
in various materials, structural parts and structures, is 
the subject of the 8-page Bulletin No. 2084 of the American 
Instrument Company, 8010-8020 Georgia Avenue, Silver 
Spring, Maryland. The instrument has been on the market 
for a number of years, but the present bulletin is the most 
comprehensive yet published. 
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GAERTNER 


Range .5461 to 15 microns 


with 


Bilateral Slits of 40 mm clear aperture and pro- 
vided with micrometer, reading from 0.01 to 3 
mm width of opening. 


Rocksalt or Potassium Bromide Prisms of best 
optical quality. 


Mirrors are parabolized 10° off axis, focal length 
1 m, figured to high accuracy. 


Baffled to eliminate stray light. 


Cover is air tight except for the breather tubes. All adjust- 
ments of slits, wavelength, etc., are made from outside of case. 


THE GAERTNER SCIENTIFIC CORPORATION 
1212 Wrightwood Ave. 


INFRA-RED SPECTROGRAPH (Littrow Type) 


Chicago, U.S.A. 


TACHOMETERS 


—centrifugal, chronometric and vibrating- 


work. Guaranteed accuracy to .1%. 


PORTABLE POTENTIOMETERS 


Rubicon portable potentiometers are available in a 
variety of single and double range models suitable 
for thermocouple, pH and instrument calibration 


These sturdy, self-contained instruments are built 
to give years of satisfactory service under all condi- 


reed s. Jagabi Hand Tachoscope and J tions encountered in the laboratory, shop or field. 


Hand Tachograph. Jagabi Tachonorm. : They are outstanding values in their class. 
Write for Bulletin 1505-R Ask for Bulletin 270. 


JAMES G. BIDDLE CO. RUBICON COMPANY 


rical nt Makers 
ELECTRICAL AND SCIENTIFIC INSTRUMENTS Electrical Instrument e 


1211-13 ARCH STREET ~ PHILADELPHIA, PA. Ridge Ave. at 35th Street Philadelphia, Pa. 
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INDEX TO ADVERTISERS 


The following firms have shown their interest 
in the advancement of physics through their 
support of the journals and other services of the 
American Institute of Physics and its Founder 
Societies. This should entitle them to consid- 
eration by physicists whenever possible. 


Name 


AMERICAN INSTRUMENT COMPANY 
Ion Type X-Ray Tube. 
Cameras. 


Powder Spectrum Diffraction 


Scientific Instruments and Supplies. 


Bauscu & Loms Optica CompANY 


Makers of 
fractometers, 


Microtomes, Colorimeters, Re- 
Spectrometers, Balopticons, Photomicro- 
graphic and Microprojection Apparatus and related in- 
struments. Also makers of Orthogon Eyeglass 
for Better Vision. 


Microscopes, 


Lenses 


JAMEs G, Bippte ComPANY 


“Jagabi’’ Rheostats; Adam Hilger and Kipp & Zonen 
Optical Instruments; ‘Pointolite’’ Lamps; Electrical Test- 
ing and Speed-measuring Instruments. 


CENTRAL Screntiric COMPANY . 
Manufacturers of Cenco Physical Apparatus and Instru- 
ments to meet all requirements of University, 

High School Physics Laboratories. Specializing in 

high vacuum pumps and development of instruments and 

apparatus for various sciences. 


( lege 
and 


EASTMAN KopAK COMPANY . 


Purified Organic Chemicals for research purposes; Plates 
for Photography, Pho- 
tometry, Astronomy;.Wratten Light Filters; Cameras and 
Films. 


Photomicrography, Spectroscopy, 


Tue Laporatory, 


Precision electrical 
instruments; potentiometers, bridges, temperature bridges, 
volt’ boxes, 


Standards of e.m.f. (standard cells). 


Thermopiles and pyrheliometers. 


GAERTNER SCIENTIFIC CORPORATION 


Spectroscopes, Spectrometers, Spectrographs, Spectropho- 
Heliostats, 
tors, Cathetometers, Reading Telescopes, Interferometers, 
Chronographs, Dividing Machines, etc. 


tometers, Measuring Microscopes, Compara- 


GENERAL RapIo COMPANY 


Manufacturers of electronic measuring instruments; 
vacuum-tube voltmeters, amplifiers, and oscillators; wave 
analyzers, noise meters and analyzers, stroboscopes; lab- 
inductance and fre- 
and 


rheostats, 


oratory standards of capacitance, 


quency; impedance bridges, decade resistors con- 


densers; air condensers and variable inductors; 


Variacs, transformers; other laboratory accessories. 


iil 


INDEX — Continued 


Name 


INTERNATIONAL RESISTANCE COMPANY 


Manufacturers of 
variable 


metallized and wirewound fixed and 
Resistors including high voltage and high fre- 


quency types. 


Leeps & NortHruep CoMPANY 


Manufacturers of Galvanometers, Resistors, Bridges, Con- 
densers, Inductances, Potentiometers, Testing Sets; 
perature Measuring, Recording and Controlling Appa- 
ratus; Instruments for Measuring and Controlling 
Conductivity of Electrolytes and Hydrogen Ion Concen- 
trations. 


Tem- 


McGraw-Hitt Book Company, INc. 


NATIONAL RESEARCH CORPORATION 


Specialists in industrial applications of high vacuum in 
the fields of physics and 
consulting; 


chemistry; research, design, 
low reflecting glass surfaces, special products 


produced under low pressure. 


Galvanometers, electrometers, potentiometers, Wheatstone 
and Kelvin bridges, resistance boxes, hydrogen ion and 
conductivity apparatus. 

Tue Society oF EXpLorAtioN GEOPHYSICISTS .... il 
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Windows Flat Over Entire Area 
to Within 6 Wavelengths 


* Parallelism of liquid-glass inter- 
face is +0.01 mm. in high-preci- 
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standard cells. 


*Parallelism of the faces of each 
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*Actual mean length marked on 
each cell to nearest 0.001 mm. 


» Actual mean length between faces 
of windows is + 1% or +0.005 mm. 
of the nominal length. 


«Cells matched to +0.01 mm. are 
obtainable. 


* Pyrex Bodies with Pyrex or Corex 
Windows. Fused Quartz 
Bodies and Windows. 


Strain-free Cells in All Pyrex or All Quartz 
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Fic. 1. Settlement of soil surface at various instants of 
time for a uniform load po» extending from x=0 to 
x= 


a tabulated function called the probability in- 
tegral. 

The settlement is composed of two parts: a 
purely elastic deflection which occurs at the 
instant of application of the load and a settle- 
ment due to consolidation which occurs gradually 
thereafter. Since we are interested only in the 
latter we proceed as follows: 

For ¢=0 the deflection is purely elastic, its 
value is 


sin \x/2G). (1.12) 
The settlement due to consolidation is 
W, = W— Wj. (1.13) 
For convenience we will write Eq. (1.13) in the 
form 
sin Ax. (1.14) 


II. SETTLEMENT UNDER A UNIFORM Loap 
WITH DISCONTINUITY 


For practical purposes we are interested in the 
differential settlement occurring when the load 
distribution has a discontinuity. Such a load 
. distribution will be represented by the following 
function of x 


f (2.1) 
WT 


in which case ¢(x)=—po0/2 for x<0, po/2 for 
x>0 with a discontinuity at x=0. 

The settlement due to this load is obtained 
from (1.14) by the principle of superposition ; 
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Fic. 2. Settlement of soil surface at various instants of 
time 4 a uniform load po extending from x=—//2 to 
x=+//2. 


Fic. 3. Rates of the settlement at the center of the 
loaded area in Fig. 2 when the top surface is pervious to 
the settlement, when the top surface is impervious. 


4: 


we find 


Ws 


f FING" Axdd, (2.2) 


which can be written in the form 
w, 1. Fly) 
=— ——-sin (yi)dy (2.3) 
with §£=x/(ct)}. 
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This integral can be evaluated numerically by 
the use of an approximate analytical expression 
for the function F(y). It may be verified that 
this function is represented within 2 percent by 


F(y) =1—e-7— 1.062e-*7 


(2.4) 


If (2.3) is integrated with the substitution 
(2.4) we find 


— ta —1 


1.06€ 0.28& 


2 + 


By adding a load equal to po/2 and extending 
from to x=+o we obtain a total 
uniform load extending from x=0 to x= © as 
shown in Fig. 1. In doing so we do not change 
the settlement because the water is imprisoned 
in the soil by the impervious layer at the surface 
and therefore a load extending from x=— ~ to 
x=-++© cannot produce any settlement. Hence 
the settlement for the load po extending over 
the positive x axis is given by the same expression 
(2.5) which we write 


= 2apo(ct/m)*f(é), 


where f(¢) denotes the right-hand side of Eq. 
(2.5). It is clear that the volume of the water in 
the soil being constant, the total area under the 
curve f(£) is zero. 

In order to represent the settlement as a 
function of time it is convenient to introduce for 
the abscissa a characteristic length / which can 
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be chosen arbitrarily. We then have 


(ct)itx 1 (2.6) 


Then settlement curves are plotted as function 


of x in Fig. 1 at time intervals corresponding to 
(ct)3/1=1/8, 2/8, 3/8, 4/8, 5/8. 


III. SETTLEMENT UNDER A RECTANGULAR 
Loap DISTRIBUTION 


By superposition we may easily derive from 
the previous solution the settlement due to a 
constant load po extending from x=-—I//2 to 


x=1/2. Using (2.6) we get 
—1/2 
)}} (3.1) 


2 (ct)? x+l/2\ x 
(ct)! 
The settlement curves are represented in Fig. 2 
as a function of x at time intervals corresponding 
to (ct)!/l=1/8, 2/8, 3/8, 4/8, 5/8. 

It will be noted that immediately after loading 
the settlement at the edge is greatly affected by 
the unloaded section adjacent to it. This un- 
loaded section provides an escape for the water 
which is not as readily available for the water 
under the center of the load. Also since this load 
and deflection are symmetrical to the line x=0 
no water will flow across this line and therefore 
the area under the curve from 0 to © must be 
equal to zero. 

If the ratio of the deflection at x=0 is com- 
pared with the case in a previous paper? in 
which the loading is the same but the water was 
allowed to flow freely from the top surface the 
effect of an impervious top surface will be seen. 
This comparison is made in Fig. 3. 
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Innovations in Instruments 


New Temperature Indicator 


Increased precision in measuring thermocouple tempera- 
tures, in checking thermocouples, and in measuring small 
e.m.f.’s, may be obtained without increase in instrument 
cost by the use of the new portable temperature indicator 
recently put on the market by the Leeds and Northrup 
Company, 4934 Stenton Avenue, Philadelphia. In the new 
indicator 8667, only a small portion of the range is calibrated 
on a continuously adjustable slidewire ; the rest, adjustable 
in fixed steps, is on a dial switch of ten highly accurate 
resistors. Range is from 0 to 111 millivolts with a limit of 


error of +0.1 millivolt. The increased accuracy without 
increase in cost is obtained by sacrificing the convenience 
of a reference-junction compensator; the only accessories 
required for this self-contained potentiometer, illustrated 
in the accompanying figure, are the thermocouple and an 
ice bath for the reference-junction. 


Insulation and Insulating Varnish 


The Irvington Varnish and Insulator Company, 24 
Argyle Terrace, Irvington, New Jersey, announced the 
following new products: Irvington Varnished Fiberglas, 
glass cloth impregnated and coated with heat resisting var- 
nishes, and used as insulation in the manufacture of electric 
motors, generators and similar equipment; and Harvel 
612-C insulating varnish, a phenol-formaldehyde synthetic 
resin-base varnish made from Cashew Nut Shell Liquid 
which polymerizes upon baking to an infusible, insoluble 
state. Literature concerning these products may be ob- 
tained by writing the manufacturer. 


Variety in Decade Resistance Boxes 


Shallcross Manufacturing Company, 10 Jackson Avenue, 
Collingdale, Pennsylvania, announces a wide range of 
decade resistance boxes, manufactured in increments from 
0.1-ohm steps up to 1,000,000-ohm steps. Some have as 
many as six decades in a single box. Details are given in 
Bulletin No. 825. 
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New Vacuum Tubes 


RCA Manufacturing Company, Inc., Harrison, New 
Jersey, announces that it is making available the following 
new tubes: 

RCA-6SG7 and RCA-12SG7, r-f amplifier pentodes of 
the metal type particularly recommended for use in high 
frequency receivers. They feature high transconductance, 
very low grid-plate capacitance, and two separate cathode 
terminals. The 6SG7 and 12SG7 are alike except for 
heater rating. The heater of the 6SG7 is designed so that 
it can be operated in series with other 6.3-volt, 0.3-ampere 
types; likewise, the heater of the 12SG7 can be operated 
in series with other 12.6-volt, 0.15-ampere types. 

RCA-45Z3, a miniature, half-wave high vacuum rectifier 
of the heater-cathode type. Designed specifically for 
service in a.c./d.c./battery-operated portable receivers, 
it offers the advantages of small size and low heat dissipa- 
tion. Its heater requires only 0.075 ampere at 45 volts and 
its output rating makes it capable of supplying rectified 
power for both filament and plate circuits in light-weight 
receivers utilizing the miniature type tubes. 

RCA-3Q4, a new, miniature type, power amplifier 
pentode suitable for operation with 90 volts on plate and 
screen. It has a power output of 270 milliwatts with 7 
percent distortion. The filament of the RCA-3Q4 has a 
center tap which permits economical usage of this tube in 
series-filament arrangements with other miniature tubes. 

The 6SS7, a remote cut-off, r-f amplifier pentode of the 
single-ended metal type having a 6.3-volt, 0.15-ampere 
heater. This new tube provides for a further degree of 
flexibility in the design of a.c./d.c. receivers utilizing single- 
ended metal types, where the total heater voltage of a 
complement of 0.15-ampere types heretofore available 
would exceed 117 volts. 

The 12SN7-GT, a single-ended, twin-triode amplifier 
having separate cathode terminals for each triode unit. 
It is recommended for use in resistance-coupled circuits 
as a voltage amplifier or phase inverter. 

RCA-930, new gas photo-tube recommended for use in 
sound reproduction and relay applications. Electrically, 
the $30 is like the type 923 with its high sensitivity and 
large response to red and near infra-red radiation. 

The RCA-931, a new type of photo-tube in which the 
photo-current produced at a light-sensitive cathode is 
multiplied many times by secondary emission occurring 
between nine successive dynodes within the tube. It is 
capable of multiplying feeble currents produced by weak 
illumination as much as 230,000 times. Focusing of the 
electron stream is accomplished electrostatically within 
the tube. The 931 employs the S4 photo-surface which has 
higher sensitivity to blue-rich light than to blue-deficient 
light. 

The 8001, a transmitting beam pentode suitable for 
power amplifier, modulator, and oscillator service. The 
8001 is particularly useful in all band transmitter designs. 
In class C telegraph service, the 8001 will provide a power 
output of approximately 230 watts at frequencies as high 
as 75 megacycles. Its plate dissipation is 75 watts. 
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